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Diffuse Scattering of X-Rays by Conduction Electrons 


CLARENCE ZENER,* Wayman Crow Hall of Physics, Washington University 
(Received July 8, 1935) 


A theoretical expression is obtained for the scattering 
of x-rays from a completely degenerate free electron gas. 
Both the classical and the Compton methods are used. 
The modifications introduced by the periodic potential of 
the lattice are examined. The scattering by the conduction 
electrons is found to decrease with the first power of the 





scattering angle. It is shown that the diffuse scattering by 
the core electrons decreases as the second power of the 
angle of scattering. Hence the scattering by the conduction 
electrons will dominate all other diffuse scattering at 
sufficiently small angles. 








$1. 

HE QUANTUM-MECHANICAL theory 

of the scattering of x-rays has been thor- 
oughly investigated.' In practice, however, we 
may wish to use classical concepts. Again, for 
incoherent scattering, we may find it convenient 
to think in terms of photons being scattered by 
single electrons, as in the Compton effect. In the 
first, the classical method, the interference be- 
tween two electrons arises from the addition of 
the amplitudes of the scattered x-rays. In the 
second, the Compton method, the interference 
between two electrons arises through the opera- 
tion of Pauli’s exclusion principle. Thus an elec- 
tron cannot jump into a state already occupied 
by another electron. 

The first method gives the correct total in- 
tensity of scattered x-rays (relativity and reso- 
nance effects neglected). The formula developed 
by Raman for the ratio S of the actual scattering 
to that which would occur if the intensities of the 
scattering from each electron were added, is* 

*The author was aided in part by a grant from the 
Rockefeller Foundation to Washington University for 
research in science. 

1See A. H. Compton, Phys. Rev. 47, 367 (1935) for a 


comprehensive review of the literature. 
?C. V. Raman, Ind. J. Phys. 3, 359 (1928), Eq. (3). 


573 






S=1+time average of N~' > cos k(z,—z,). (1) 


Tx 8 


Here N is the total number of electrons, and 
k=(42/X) sin ¢/2, 


where J is the wave-length of the incident x-rays, 
and ¢ is the angle of scattering. Wave mechanics 
gives the formula’ (relativity and resonance 
effects neglected) 


S= 1+ fy > cos k(z,—2,)~*dvy---dvy. (2) 

rs 
These two expressions for S are identical if we 
identify the classical time average with the 
average over y*. 

If the probability that one electron were at a 
given position did not depend upon the positions 
of other electrons, ¥ could be written as a product 
of eigenfunctions of single electrons 


y= Ily.(r., 4,), 


where 6, denotes the spin variable. In this case 
(2) reduces to 

37. Waller and D. R. Hartree, Proc. Roy. Soc. Al24, 124 
(1929), Eq. (14). 
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We shall use the notation 
fan = Sn Cos key, *dv 


with the convention that the integral is zero when 
n and n’ refer to states with opposite spins. The 
second term usually represents constructive 
interference. Thus two electrons bound to the 
same atom may move independently of one an- 
other, but their mere confinement to the same 
atom gives constructive interference, at least at 
small angles. 

The probability of one electron being at a 
definite position is dependent upon the presence 
of the other electrons by virtue of Pauli’s exclu- 
sion principle. This principle is taken explicitly 
into account by writing the wave function in the 
form 


y=>d(- 1)*rPIYys(s, 5s). 


Here P is an operator which permutes the coordi- 
nates (including the spin coordinate) of the elec- 
trons, and gp is the order of the permutation. 
With this wave function (2) becomes* 


S=1+N" > (fant nin’ —Sinn’)- (3) 


nn’ 

The last term is intrinsically negative. Its 
interpretation depends upon our picture of the 
scattering process. According to the Raman 
viewpoint, Pauli’s exclusion principle lessens the 
constructive interference between electrons of 
like spin. According to the Compton viewpoint, 
Pauli’s exclusion principle prohibits those colli- 
sions between photons and electrons in which the 
recoil electron would find itself in a state already 
occupied by another electron of like spin. 

The first purpose of this paper is to discuss in 
detail the scattering of x-rays in a case where both 
methods may be carried out exactly (§2 and §3). 
This case is the scattering by a free electron gas. 
The second purpose is to examine the scattering 
by the conduction electrons of actual metals 
(§4). Finally, we investigate the conditions 
where this scattering may be observed (§5). 


§2. 
From the classical standpoint the scattering of 
x-rays from a free electron gas presents the same 
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problem as the scattering from a monatomic 
liquid. In the first problem we must find the 
interference between the individual electrons, 
while in the second we must find the interference 
between the individual atoms. The following 
analysis for a free electron gas follows step by 
step the well-known analysis‘ for a monatomic 
liquid. 

The probability that a particular electron lies 
within the volume element dz,, at the same time 
that another particular electron lies within the 
volume element d?,, is 


P(dv»,/ V)(dv,/V), 


where V is the total volume of the electron gas. 
Since the medium is isotropic, P is a function 
only of the distance r between the two volume 
elements, i.e., between the two electrons. Eq. (1) 
may now be transformed to 


oo sin kr 
S= tnt (1 —P)———r*dr, (4) 


0 kr 


where 1 is the density of electrons. 

If the electrons moved independently of one 
another, then P=1, and the interference would 
average out to zero. Actually, Pauli’s exclusion 
principle tends to keep electrons of like spin 
apart. Since the configurations favorable to con- 
structive interference are thereby reduced, the 
average interference will be destructive, at least 
for small angle scattering. 

The function 1—P for a free electron gas is® 


(9) 





sin (r/d) —(r/d) cos (7, “) : 


1—-P=(9 2( 
(r/d)* 


Here d=) o/27, where Xo is the wave-length of the 
fastest electrons, namely, 2(32/7)~8. 1—P has 
been plotted in Fig. la. It has the important 
property that 


nan fy (1—P)r*dr=1, 


i.e., the integral of the ‘electron hole’’ density is 
equivalent to one positive electron. Hence scat- 
tering in the forward direction will be completely 
destroyed by interference. 

4 Compton and Ellis, X-Rays in Theory and in Experi- 
ment (1935), pp. 177-185. 


5 Reference 4, Eq. (3.80). 
® E. Wigner and F. Seitz, Phys. Rev. 43, 807 (1933). 
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Fic. 1. Classical representation of scattering by free 
electron gas. The ‘‘electron hole’’ density (a) gives rise 
by destructive interference to an S value (b) less than 
unity. 


In order to obtain an explicit expression for S, 
we have merely to substitute (5) into (4). We 
obtain (see Fig. 1b) 


(1/2)(3x—x*), x<l 
| (6) 


1 , ae 


with x=kd/2=(Xo/A) sin (¢/2). 

The method of performing the integration in 
(4) will be briefly sketched. We obtain with ob- 
vious transformations 


3 © /sin Z—2 COS Z\? | 
S=1-{ — { e*2d¢. 
2rix] J_.,, 2* 


Since the integrand has no poles in the complex 
plane, we may choose any path for the integra- 
tion, say the path IT, of Fig. 2. Express the inte- 
grand as the sum of three terms, containing the 
factors exp 2i(x+1)s, exp 2txz, exp 2i(x—1)z, 
respectively. The integral of the first two terms 
vanishes, since no pole lies between the path I; 
and an infinite semicircle in the upper half- 
plane, on which the two terms vanish. If x>1, 
the third term likewise vanishes. Now the inte- 
gral of the third term along I is equal to — 277 
times its residue at z=0, plus its integral taken 
along the path ls. If x <1, this integral along Tl. 
vanishes, since the integrand is zero along an 
infinite semicircle in the lower half-plane. We are 











* 
a Ft ews 
@eeee —z Z-0 aeee 
. 


Fic. 2. Paths of integration used in evaluating S. 
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thus left with — 277 times the residue of the third 
term at z=0. This leads us directly to (6). 


$3. 


In the case of a free electron gas, Eq. (3) re- 
duces to 
S=1-N—“T fan’. 
ngtn’ 
A simple phenomenological interpretation of this 
equation is the following. 

We regard the scattering of x-rays as the simple 
Compton scattering of photons by individual 
electrons. The intensity of scattering in any 
direction is given by the Thomson formula (apart 
from relativity corrections). During the collision 
the scattering electron suffers a change of momen- 
tum which is calculable from the principles of 
conservation of energy and of momentum. Since 
the wave properties of the photon are not con- 
sidered, constructive and destructive interference 
between two electrons does not occur. However, 
we must recognize that the electrons obey Pauli’s 
exclusion principle. We thus exclude those colli- 
sions for which the final quantum state of the 
electron is already occupied. In the case of a 
(completely degenerate) Fermi gas, the final 
momentum must thus lie outside the Fermi sur- 
face in momentum space. 

The scattering factor S may thus be calculated 
as follows. We calculate the momentum Ap, 


| Ap| =(2hv/c) sin (¢/2), 


which an electron must absorb to scatter a 
photon through the angle ¢. Then S(¢@) is equal 
to that fraction of the volume inside the Fermi 





(b) «) 


Fic. 3. Compton representation of scattering by a free 
electron gas. The momentum change of a photon (a) during 
a collision is accompanied by an opposite change in mo- 
mentum of the scattering electron (b). Electrons at such 
positions as A in momentum space are effective, while 
those in such positions as B are not effective. (c) is a 
schematic diagram analogous to (b) for the case of con- 
duction electrons. The full lines represent the Fermi 
surface. The dotted square represents the boundary of the 
(multiply connected) first zone. 
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surface for which Ap+p lies outside the Fermi 
surface (see Fig. 3). A simple calculation shows 
that this fraction is given by Eq. (6). 


wre 


4. 


In an actual metal the conduction electrons are 
not entirely free. A quantum state may no longer 
be specified by its momentum. Instead, each 
quantum state is specified by a vector k which, 
however, has many properties analogous to that 
of momentum.’ In particular, if a photon changes 
its momentum by Ap in colliding with an electron, 
the electron must change its k vector by an 
amount Ak= —Ap. Thus only those states k are 
effective in this particular type of scattering for 
which the state kK—Ap is unoccupied. 

The only difference between conduction and 
free electrons which we need consider is that in 
the former the Fermi surface in k space is not 
necessarily a sphere. In fact, k space is divided 
into an infinite number of zones, each zone being 
multiply connected.’ At small angle scattering 
only transitions between states in the same zone 
are important. 

As the angle of scattering is decreased it is 
thus evident that S becomes proportional to the 
area of the Fermi surface. In divalent metals 
there are just enough conduction electrons to fill 
the first zone. Hence as the binding in these 
metals becomes stronger (case of poor conduc- 
tors) the area of the Fermi surface, together with 
S, becomes smaller. For extremely tight binding 
(case of insulators) the area of this surface 
vanishes. 

7 An excellent review of the quantum theory of conduc- 
tion electrons is given by H. Bethe, Handbuch der Physik 


(1933), Vol. 24, No. 2. 
’ For example, see reference 7, Fig. 19. 
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5. 


It remains to investigate the conditions under 
which diffuse scattering from conduction elec- 
trons may be experimentally observed. It is 
necessary of course that the scattering due to the 
conduction electrons be comparable to the scat- 


Se 


tering due to the core electrons. 
The diffuse scattering due to each core is given 
by® 


I=1,{(2'-f?/Z’)+(1-e™) f}. 


Here J, is the Thomson scattering of a single 
electron, Z’ is the number of electrons per core. 
The structure factor f is given approximately by 


(4r/k)fo°u(r) sin (kr)rdr, 


where u(r) is the density of the core electrons, 
and k= (47/X) sin (¢/2). For our present purpose 
we need only know that the temperature de- 
pendent Debye-Waller term . contains k® as 
a factor. 

Inspection shows that a Taylor expansion of J 
starts with k®. The expansion of the first term 
starts with k? since f(k) is an even function of k, 
and f(o)=Z’. That the expansion of the second 
term starts with k? is evident from the dependence 
of MW upon k. 

Hence the diffuse scattering from the core 
electrons at small angles varies as ¢*. From Eq. 
(6) we see that the diffuse scattering by the con- 
duction electrons varies as ¢. Hence at sufficiently 
small angles the latter will become greater than 
the former. Calculation shows that for the par- 
ticular case of lithium at room temperature the 
two types of scattering become equal at scattering 
angles such that 


(A/d) sin (¢/2)=0.04. 


°G. E. M. Jauncey, Phys. Rev. 37, 1193 (1931); Y. H- 
Woo, Phys. Rev. 38, 6 (1931). 
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Extra Negative Term in the Incoherent Part of the Diffuse Scattering of X-Rays 
from Neon-Like Crystals 


G. E. M. JAuncey* anp J. H. DeminG, Washington University, St. Louis 
(Received July 8, 1935) 


Following the suggestion of Harvey, Williams and 
Jauncey that experimental evidence for the need of the 
Waller-Hartree extra negative term which arises from the 
Pauli exclusion principle can best be obtained in the 
diffuse scattering of x-rays from crystals, careful measure- 
ments of the Sexp values for NaF, MgO and SiC have 
been made at 295°K using the method of Jauncey and 
Claus. Various precautions are described. Since the the- 
oretical St, values depend somewhat critically upon the 
characteristic temperature © for a given crystal, a dis- 
cussion of the determination of © from specific heat data 
is given. It is shown that for x-ray purposes the proper 
value for © is that obtained from specific heat data for 


fairly low temperatures even though the temperature at 
which x-ray scattering occurs may be fairly high. Tables 
of quantum-mechanical f and SZ, values extend only to 
(sin ¢/2)\=1.1. Beyond this the crystals approach neon 
gas in their scattering and it was found that Brown’s 
analytical formula for neon gives better agreement with 
experiment than does Froman’s extrapolation formula. 
The experimental results showed excellent agreement with 
the Waller-Hartree theory and thus give definite evidence 
for the operation of the Pauli exclusion principle on the 
outer shell electrons. Some remarks on the classical 
analog of the extra term are included. A special discussion 
of the results for SiC is given. 





1. INTRODUCTION 


HE theory of the diffuse scattering of x-rays 
from crystals leads to the formula! 


= . + Sine (1 +a vers >)’, (1) 


in which S,,, and Sj,- refer respectively to the 
coherent and incoherent parts of the scattered 
rays. For crystals consisting of two kinds of 
atoms 


Scoh = fri —_ M1) + f-?(1 —e~ 2)! (Z,:+Z2) (2a) 


and Waller and Hartree’s theory of x-ray scat- 
tering’ requires that 


Sine= 1— {D0(E;;2) 1 +d (K;j;7)2} (Z:+Z2) 
—{> | Eel P2+D | Bjx|2?}/(ZitZe2), (2b) 


j#k jek 


in which the symbols have previously been 
defined* and the numerical subscripts refer to 
the two kinds of atoms. Because of its algebraic 
sign, the term containing the subscripts jk 
has been called the extra negative term.‘ Harvey, 
Williams and Jauncey*® have pointed out that 
the most favorable situation for showing the 


* The senior author was aided in part by a grant from 
the Rockefeller Foundation to Washington University for 
research in science. 

1G. E. M. Jauncey, Phys. Rev. 42, 453 (1932). 

2 1. Waller and D. R. Hartree, Proc. Roy. Soc. A124, 119 
(1929), 

3G. G. Harvey, P. S. Williams and G. E. M. Jauncey, 
Phys. Rev. 46, 365 (1934). 

4E. O. Wollan, Phys. Rev. 43, 955 (1933). 


a) 


need for the extra negative term is in the diffuse 
scattering of x-rays from crystals. In addition, 
they have calculated the values of =| #;,/|* for 
neon. Values of f, the atom form factor, and of 
YE;* have been calculated by James and 
Brindley’ and by Heisenberg,® respectively. 
Tables of these values are now given by Compton 
and Allison.’ We have used the values for Mg*?, 
Nat, F-, O- in calculating the theoretical 
values for our crystals of NaF and MgO. 

Compton* remarks that there seems to be no 
exact classical analog of the extra negative 
term, although he suggests that it represents the 
constraints upon each electron’s motion due to 
the presence of the other electrons in the atom. 
At a recent symposium on x-ray scattering, 
Jauncey® also suggested this idea of constraints 
between electrons. In his extension of the 
Raman-Compton!’:" classical theory of x-ray 
scattering, Jauncey” introduced the. probability 
distribution function p(2), 22, - ++ 2,)d2,d2q: - «dz, 
for the probability of electrons 1, 2, --- being 
in dz,, dz_ +--+ dz, at 21, 22, Zn, respectively, 
from the center of the atom. Jauncey then made 
the simplifying assumption 


> R. W. James and G. W. Brindley, Phil. Mag. 12, 81 
(1931). 

®W. Heisenberg, Physik. Zeits. 32, 737 (1931). 

7 A. H. Compton and S. K. Allison, X-Rays in Theory and 
Experiment, pp. 781, 782. 

* A. H. Compton, Phys. Rev. 47, 367 (1935). 

G. E. M. Jauncey, Phys. Rev. 47, 196 (1935). 

”C.V. Raman, Ind. J. Phys. 3, 357 (1928). 

tA. H. Compton, Phys. Rev. 35, 925 (1930). 

2G. E. M. Jauncey, Phys. Rev. 37, 1193 (1931). 
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p21, Z2, ++ Sn) = Pili): po(Ze): ++ PalSn). (3) 
This means that the probability distribution 
function for any one electron is independent of 
the position of any other electron in the atom. 
This assumption leads in the case of crystals 
consisting of atoms of one kind or of monatomic 


gases to 


Z 
Sine= 1 = (> E,*)/Z, (4) 
r=l1 
where for the rth electron 
E,= JS p, cos kz,dz, (5) 


and k=4n(sin ¢/2)/d. If, however, the prob- 
ability distribution function for one electron is 
not independent of the positions of the other 
electrons in the atom, the relation (3) no longer 
holds and Sj, is no longer given by (4). As an 
example of the change in Sj, caused by intro- 
ducing constraints, let us consider the two 
models of the helium gas atom discussed by 
Compton." In the first model each electron is at 
random on a sphere of radius 0b. In this case 


(21, 22) = pr(21) - p2(Z2), (6) 
where 
pi(ti) = po(ze)=1/26, —b<s<b a 
=0, |s| >b 
and we arrive at 
z 
f= LE= 2(sin kb)/kb, (8) 
Scon= f?/Z= 2(sin® kb) /k*b? (9a) 
and 
Sinc= 1— (sin? kb) /k?b?. (9b) 


In the second model, each electron is again on a 
sphere of radius ) but the two electrons are 
always at opposite ends of a diameter. Now, 
although (7) is true, the relation (6) is not true. 


For the second model, Debye" gives 
S=1-+(sin 2kb)/2kb. (10) 


'’ P. Debye, Ann. d. Physik 46, 8093(1915). 


AND J. H. 


DEMING 


According to the classical theory the coherent 
part of the rays scattered from an atom can be 
eliminated in all directions, except the Laue 
directions, by assembling many like atoms into a 
perfect lattice.'* The scattering which does occur 
in the non-Laue directions is then incoherent. 
Using this criterion, we find that 


Sinc= 1+ (sin 2kb) /2kb—2(sin? kb) /k*d?. (11) 
This expression differs from (9b) because of the 
electrons being constrained to move at the 
opposite ends of a diameter. Each model has 
the same atom form factor (8) and so for either 
model of helium gas S.on is given by (9a). 

A comparison of Williams’ experimental S 
values’ for NaF with theory indicates a better 
agreement when the extra negative term is 
included than when not.* However, the dis- 
crepancy between theory and experiment is too 
great for comfort ; and, in view of the importance 
of the theory, a more accurate test of the formula 
was considered necessary, and so this present 
research was undertaken. 


2. EXPERIMENTAL METHOD 


Following the method described by Jauncey 
and Claus,'® we used the continuous spectrum 
from a tungsten target tube operated at 34.3 kv 
peak and 8 ma. As always, great care was taken 
in evading Laue spots. The distribution of in- 
tensity in the continuous spectrum was obtained 
by reflection from calcite and is shown in Fig. 1. 
This curve is uncorrected for change in the 
reflecting power of the calcite for different wave- 
lengths. Thicknesses of aluminum sufficient to 
reduce the intensity of the primary rays from 
1.00 to 0.86, 0.61, 0.42, and 0.22, respectively, 
were determined experimentally. The reduction 
of total intensity in the continuous spectrum due 
to each of these thicknesses was calculated as 
previously described'® and found to agree very 
well with each of the respective values just 


4 This is evident from the fact that the double summation 
in Jauncey and Harvey’s paper (Phys. Rev. 37, 1203 
(1931)) becomes zero when the atom centers are always 
exactly at lattice points and there is thus no thermal 


motion. 
1% P. S. Williams, Phys. Rev. 46, 83 (1934). 
(4. E. M. Jauncey and W. D. Claus, Phys. Rev. 46, 941 


(1934). 
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Fic. 1. Spectral distribution curve for the x-rays falling 
upon, but not penetrating through, the crystals. 





given. It seemed better not to attempt any 
correction of the curve of Fig. 1. 

Instead of comparing the intensity of the rays 
scattered from the crystal under observation 
with that of the rays scattered from paraffin as 
was usually done by Jauncey and Claus,'® we 
always compared the intensity of the diffusely 
scattered rays with that of the primary rays 
penetrating the crystal. To do this we cut off a 
flash of primary rays of duraticn 1,60 sec. by 
means of a pendulum. We then calculated 


Sexp = (R?W/A NZ pt) -(2 cos ¢/2)/(1+ cos? ¢) 
-(m®ct/et)-D,/D, (12) 


where D, is the ion current due to the rays 
scattered in a direction ¢, D is that due to the 
primary rays transmitted through the crystal 
slab, and the remaining quantities have pre- 
viously been defined.’ It is to be noted that the 
ratio D,/D is of the order 1/10,000. The dis- 
parity in the concentration of the ions in the 
ionization chamber is even greater than this 
because of the fact that the primary beam is 
fairly narrow. It was feared that there might be 
some unsuspected ion-recombination effect. We 
therefore reduced D from 1.00 to 0.66, 0.25, and 
0.11 and obtained no certain change in the ratio 
D,/D, as compared with the corresponding 
theoretical value of D,/D, for a scattering angle 
of 30°. Since D,/D is not a function of D, we 


17. E. M. Jauncey and Ford Pennell, Phys. Rev. 43, 
505 (1933). 














0 1 rl i | l i 
sin 5 
Fic. 2. Comparison of theory and experiment. Full 


curves include the Waller-Hartree extra negative term; 
broken curves do not. 





conclude that the experimental D, and D are 
saturation values. Of course, we had, in addition, 
applied the usual test for saturation by varying 
the voltage across the ionization chamber. We 
estimate our experimental error at about 0.02 
in the S.x, values. 

The values of S.., for single crystal slabs of 
NaF, MgO and SiC (carborundum) are given in 
Table I and are shown as circles in Fig. 2. The 


TABLE I. Sexp values. 





sin ¢/2 Nak MgO SiC 





0.087 0.56 0.46 0.29 
.104 0.68 0.55 — 
130 0.86 — 51 
.148 0.93 72 — 
174 1.06 .76 67 
.216 1.11 87 .76 
.259 1.12 93 .84 
342 1.08 .93 .89 
462 1.03 .94 88 


707 0.93 .88 85 
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values of pf for the crystal slabs were: NaF, 
0.349; MgO, 0.681; SiC, 0.326 g cm®. It is to be 
noted that our S,,, values are all absolute and 
that there has been no fitting of the experi- 
mental and theoretical values at any point. 


3. COMPARISON WITH THEORY 


The theoretical quantity S,), is defined by!” 


Sinchk ol ¢ 
Sax f {Son — Sr eieaetiae “ tan / { 1n 
K.(1+e vers o)3) (13) 


where J is a function of \ representing the 
distribution of intensity in the spectrum of the 
X-rays penetrating the crystal slab, so set that 
the normal to the slab bisects the angle of 
scattering. The spectral distribution curve shown 
in Fig. 1 is not the curve of J vs. \ used in (13) 
but is the distribution curve for the x-rays 
falling on the slab. Knowing the absorption 
coefficients for various wave-lengths in a given 
crystal slab and knowing its effective thickness 
t/(cos ¢/2), the curve of J vs. X for the rays 
penetrating the crystal can easily be calculated 
from Fig. 1. There is a separate curve for each 
angle of scattering and for each crystal. 

From (2a) we see that it is necessary to know 
the Debye-Waller exponents 7; and /.. In the 
case of cubic crystals consisting of atoms of 
nearly equal mass, J; may be taken as equal to 
Mz. This is the case for NaF. Also 1/,;~ 1/2 when 
the experimental temperature, which for our 
experiments was 295°K, is considerably lower 
than the characteristic temperature of the 
crystal. This is the case for MgO and SiC. 
Waller’s theory'® shows that at temperatures 
below the characteristic temperature it is the 
mass per unit volume!’ which is involved in 
the formula for 1/7, and so we have used the 
average atomic weight in this formula. But, 
before M can be calculated, the characteristic 
temperature © must be known. Shonka®’ has 
obtained © = 442°K for NaF by comparing the 
intensities of the various Debye-Scherrer circles 
for the powdered crystal at 295°K with the in- 
tensities of the corresponding circles for 86°K. 


18], Waller, Diss. Uppsala (1925). 

19 Dr. C. Zener of this laboratory has informed us of this 
interesting theoretical fact. 

20 J. J. Shonka, Phys. Rev. 43, 947 (1933). 
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Unfortunately, we are unaware of any x-ray 
determination of © for MgO or SiC. We find that 
the uncertainties in the values of © as calculated 
from specific heat data lead to variations in S,, 
which are greater than our experimental error in 
Sexp and are a considerable fraction of the 
change in S;, expected as the result of the in- 
clusion or otherwise of the Waller-Hartree extra 
negative term. 

In order to obtain a method of calculating © 
from specific heat data which will be valid for 
X-ray purposes, we have taken the tabulated 
specific heat values for NaCl and KCI! and for 
Nak .2 We have changed these values from C, 
to C, by means of the formulas” 


Cp—Cy= C,?2T X.0214/Tn, (14) 
where 7, is the melting point, and 
Crp—C.=8°MT/ Kop, (15) 


where 8 is the coefficient of volume expansion, 
M the molecular weight, K the compressibility 
and p the density of the crystal. The values of C, 
for our crystals determined by either formula are 
in good agreement, but lack of extensive thermal 
data makes (14) the more useful. We then com- 
pared the C, values with a table of Debye’s 
specific heat function** and found the corre- 
sponding 0/7 values. Knowing the temperature 
for each of these values, we have calculated 
values of © and have then plotted curves of 
© vs. O/T as shown in Fig. 3. It is seen that the 
points give evidence of falling on curves which 
approach a constant value. The dotted lines 
show the values of © which are satisfactory for 
X-ray purposes as found by Shonka for NaF, by 
James and Brindley® for KCI, and by James and 
Firth®® for NaCl. It is seen that the x-ray value 
of © for each of these crystals corresponds to the 
constant value which is obtained from specific 
heat data at fairly low temperatures. The 
specific heat values of 0 for MgO*! are also shown 


21 Int. Crit. Tables, V, 91. 

2 Int. Crit. Tables, V, 100. 

*3 See F. K. Richtmyer, Jntroduction to Modern Physics, 
(first edition) pp. 244, 245. 

*4See J. K. Roberts, /Jeat and Thermodynamics, p. 404. 

2° R. W. James and G. W. Brindley, Proc. Roy. Soc. 
A121, 155 (1928). 

ant W. James and E. Firth, Proc. Roy. Soc. A117, 62 
(1927). 
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Fic. 3. Characteristic temperature curves obtained from 
specific heat data. The value of © shown in the body of 
each graph is that most suitable for x-ray purposes. Errata: 
The ordinates should be labeled © instead of T. 


in Fig. 3. It seems, therefore, that the x-ray 
value of © for MgO should be 734°K. 

The rest of the calculation of the S,;, values 
was done as described by Jauncey and Claus,'® 
these values being averaged for ionization 
chamber window width and height and for cross- 
sectional area of the primary beam of x-rays 
at small angles of scattering. The full curves of 
Fig. 2 show the S,;, values in which the Waller- 
Hartree extra negative term has been included, 
while the broken curves are for the omission of 
this term. The agreement of experiment with 
the Waller-Hartree theory is remarkably good. 
It should be noted that it is necessary in the case 
of MgO to use 0 = 734°K and not 674°K (shown 
as the lowest value in Fig. 3) in order to get an 
agreement between theory and experiment. If 
one assumes the Waller-Hartree theory to be 
correct, this means that the x-ray value of © for 
MgO is about 734°K. Our results also imply the 
excellence of the quantum-mechanical f and SE, 
values as calculated approximately by James and 
Brindley and by Heisenberg.*? 

*7 These values do not extend bevond (sin ¢/2)/A=1.1 
and some kind of extrapolation becomes necessary at the 
larger angles of scattering. We have fcund that Brown's 
(Phys. Rev. 44, 214 (1933)) analytical formula for neon 
gives much better agreement with experiment for all three 


crystals than does Fkroman’s extrapolation formula (Phys. 
Rev. 36, 1339 (1930)). 


The experimental points and theoretical curves 
for SiC are worthy of comment. The ions in 
Nak and MgO are neon-like, but there is no 
certainty that the Si and C atoms become 
quadruply ionized and therefore neon-like. Also 
we were hard put to it to determine © for SiC. 
The structure of this crystal is more complicated 
than that of NaF or MgO. Perhaps for this 
reason, when we tried to determine © by the 
method of Fig. 3, we found no tendency for the 
points to fall on a curve similar to those of 
Fig. 3. We therefore resorted to an average 
value which we took as 1023°K. Using this 
value and assuming the ions C~* and Si*‘, we 
obtained the theoretical curves shown.** The 
numerical values of the Waller-Hartree extra 
negative term as worked out by Harvey, 
Williams and Jauncey’® are for neon, and in this 
paper we have assumed these neon values to be 
valid for neon-like ions. The formulas given by 
Waller and Hartree only apply to cases where 
the outer shell of electrons is complete. We do 
not know how to correct for the operation of the 
Pauli exclusion principle in other cases, and so 
we have made no attempt to calculate values of 
Sin for SiC for possibilities other than Si** and 
C~*, It is probably a coincidence, but it is 
possibly of interest, that the experimental points 
for SiC fall so well on the Waller-Hartree curve 
for the quadruply ionized atoms. 

The authors wish to thank Dr. Preston Harris, 
Chemistry Department, Ohio State University 
for the crystal of MgO and Mr. R. C. Brenner, 
Director of the Research Laboratory, Car- 
borundum Company for the large crystals of SiC. 


28 The following f and SE,* values were necessary for our 
calculation of St, values but were not available in the 
literature. Using the information contained in James and 
Brindley’s paper® we have obtained 
(sin @/2) X 

0.0 0.1 0.2 0.3 04 0.5 0.6 0.7 0.8 0.9 1.0) 1.1 


f(C-4) 
' 909 62 3.7 26 19 1.7 16 14 13 12128 069 
SE?(C-4) 

10.0 5.8 3.1 1.65 1.4 1.3 1.0 0.8 0.7 0.6 0.5 0.4 
TE (Si**) 


10.0 9.4 8.3 68 5.2 4.0 3.1 2.5 2.0 1.8 1.65 1.5 
These values should be useful to other workers in this field. 
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An experimental and theoretical investigation of the 
absorption spectrum of MgO in the near infrared is re- 
ported. Part A is devoted to the experimental results, 
which were obtained by use of a rocksalt spectrometer. 
They show that the fundamental absorption is accom- 
panied by a great deal of secondary structure. In Part B 
an interpretation of the results of A is undertaken on the 
basis of a quantum-mechanical treatment of anharmonic 


INTRODUCTION* 


N the work to be described in the following 

paper, both experimental and theoretical in- 
vestigations were carried on and are published 
here simultaneously. In order to keep them 
somewhat apart, we have decided to present the 
details of the first in Part A and those of the 
second in Part B, and to present a general ac- 
count of their connection in the introduction and 
the conclusion. 

As a result of conclusions drawn from the 
classical theory of crystal structure of Madelung, 
Born, von Karman and other workers in the 
field,! it was generally believed until the past few 
years, that a cubic crystal possessed but one ab- 
sorption frequency, corresponding to a motion in 
which like ions move together and unlike ions 
are 180° out of phase. For the alkali-halide crys- 
tals, this fundamental should be expected to vary 
inversely as the effective mass mme2/(m,+ me) of 
the two ions, and in fact, it was found that those 
corresponding to LiF, NaCl and KBr are at 
wave-lengths 32.6, 61.1 and 83.3u, respectively. 
It is as a result of this fact that these crystals 
show a large transparency in the near infrared, 
KI being transparent as far out as 33u. 

Because of the supposed simplicity of their 
spectra, it was concluded that they should be 


* The experimental section of this paper (Part A) is 
principally due to R. B. Barnes and R. R. Brattain, and 
the theoretical section (Part B) was developed principally 
by F. Seitz while a fellow in theoretical physics at Prince- 
ton University. 

1E. Madelung, Nachr. Gétt. M. P. Klasse, March 1909 
and Jan. 1910; Physik. Zeits. 11, 898 (1910); M. Born and 
Th. von Karman, Physik. Zeits. 13, 297 (1912); M. Born 
and M. Goeppert-Mayer, Geiger-Scheel Handbuch der 
Physik, Vol. 24b. 


potential forces. On this basis, it is shown that a cubic 
crystal of the MgO type should not possess but a single 
near infrared absorption as previous calculations had 
indicated, but should possess an extremely complicated 
secondary structure. In Part C the possibilities of obtaining 
a quantitative correlation between experiment and theory 
are discussed and it is concluded that the former must be 
made the object of additional refinements. 


ideal substances on which to test the classical 
theory of dispersion in the infrared, and in 1930 
Czerny® presented the results of his attempt to 
use NaCl for this purpose. He reported the pres- 
ence of a weak secondary maximum of absorption 
upon the short wave-length side of the funda- 
mental. A similar result was obtained for KCI. 
In a second paper* it was shown that these 
crystals possessed not one but several secondary 
maxima, giving rise to absorption curves which 
differed considerably from those predicted on the 
basis of the previously mentioned theories, and 
only by using a layer of crystal 0.174 thick were 
the authors able to obtain a curve of the theoret- 
ical type. Fig. 1, taken directly from the above 
paper shows the type of curve obtained for NaCl. 
Continuing this work, Barnes* showed that some 
fifteen crystals of this type possessed similar 
secondary maxima. At this time a good deal of 
structure was observed, but since the resolving 
power was very low, the work published was to 
be regarded as being essentially envelopes of 
possible additional structure. 

Following this, Born and Blackman* sought an 
explanation of these results on the basis of an- 
harmonic terms in the potential function of inter- 
action of the lattice atoms, and handled the 
problem of a linear diatomic chain along classical 
mechanical lines. Their results showed that such 
terms would give rise to additional maxima, one 


on the short wave-length side of the fundamental 


2M. Czerny, Zeits. f. Physik 65, 600 (1930). 

* R. Bowling Barnes and M. Czerny, Zeits. f. Physik 72, 
447 (1931). 

+R. Bowling Barnes, Zeits. f. Physik 75, 723 (1932). 

5M. Born and M. Blackman, Zeits. f. Physik 82, 551 
(1933); M. Blackman, Zeits. f. Physik 86, 421 (1933). 
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Fic. 1. Transmission of thin films of NaCl. (Barnes and 
Czerny.) 


and one on the long wave-length side (the former 
agreeing reasonably well with the experimental 
values), so that a qualitative check of the validity 
of this assumption was obtained. 

When good specimens of MgO became avail- 
able, it was decided that a more thorough in- 
vestigation of the structure of the absorption 
spectrum of a cubic crystal could be carried 
through than had been done previously.* This 
crystal is particularly suitable for near infrared 
investigation since its spectrum lies closer to the 
visible than does that of any other cubic crystal, 
and, moreover, it may be obtained in perfect 
crystalline as well as powdered form. Under these 
favorable circumstances an investigation was 
carried through and the principal experimental 
results may be tabulated as follows: 

1. The absorption spectra of (a) powdered, (b) evap- 
orated films and (c) cleavage plates of crystalline MgO 
are identical as far as maxima are concerned as long as 
the particle size of the powder is not such as to give rise 
to anomalous scattering or refracting effects. 

2. The fine structure is reproducible on various thick- 
nesses of material and indicates that there are more than 
forty absorption peaks on the near infrared side of the 
largest absorption peak when the crystal is at room 
temperature. 

3. The fundamental absorption line does not lie at 
14.2u as has been previously supposed, but is further in 
the infrared. 


® R. Bowling Barnes and R. Robert Brattain, Phys. Rev. 
47, 416 (1935). 
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From this work, it was concluded that the re- 
sults of Born and Blackman were not general 
enough to handle the case at hand and another 
approach to the theoretical side of the problem 
was undertaken. For this purpose it was assumed 
that the crystal could be handled as a three-di- 
mensional quantum-mechanical system of inter- 
acting particles arranged at the mesh points of 
the NaCl type lattice and a general potential 
function of interaction was set up using only the 
group theoretical restrictions required by sym- 
metry. This is to be contrasted with the type of 
potential function employed in the Born treat- 
ment of ionic crystals in which central forces are 
assumed between particles. The recent high pre- 
cision work on the elastic coefficients by Quimby, 
Balamuth, Rose and Durand? indicates conclu- 
sively that the Cauchy-Poisson relations are not 
valid, the ci. of Voigt being only about half of the 
corresponding c4; for divalent ionic crystals of the 
type of MgO, while for monovalent crystals of 
the type of NaCl, cy is ten percent larger than C4, 
in the neighborhood of absolute zero. Since 
equality of ci. and ¢44 was essential for the validity 
of the Born theory these facts point to a basic 
lack of generality in it and for this reason it is 
found necessary, for the time being, to employ a 
phenomenological description of potential fields 
in much the same way that is done in handling 
molecular potential functions. 

By employing the anharmonic terms as per- 
turbations, normal coordinates were introduced 
so as to afford a suitable system of quantum- 
mechanical variables and from these the first- 
order energy spectrum arising from the harmonic 
terms was determined. It is of course identical 
with the classical energy spectrum and for the 
diatomic face-centered lattice is included inside 
of two zones, in each of which there are three 
branches corresponding to the three directions of 
polarization of waves. The perturbing terms were 
then regarded as being such as to alter the selec- 
tion rules for transitions without appreciably 
affecting the energy spectrum, and it was found 
that this effect can be regarded as important only 
for transitions in which two waves whose momen- 
tum vectors lie in the neighborhood of the sur- 


7 L. Balamuth, Phys. Rev. 45, 721 (1934); we are deeply 
indebted to these people for the use of results of their work 
that is as yet unpublished. 
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Fic, 2. Rocksalt spectrometer and amplifier. 


faces of the zones and at regions near these 
surfaces where the energy gradient is zero, are 
simultaneously excited, since in this case the inte- 
grated effect of many insignificant absorption 
lines is appreciable in the vicinity of one wave- 
length. In this case both the sum and difference 
energies of the two waves will give rise to ab- 
sorption peaks. The conclusions drawn were as 
follows: 


1. The first-order anharmonic terms give rise to a large 
number of absorption lines of types to be described in a 
later section. 

2. Any attempt at determining correlations between 
the experiment and theory at this stage of the develop- 
ment of the former is premature since the conditions under 
which the theory is presented are such that they will be 
realized best at low temperatures and it is only when such 
experiments are carried out that correlations may be 
sought for. There is no reason for expecting that such 
correlations do not exist, however. 


Part A. EXPERIMENTAL SECTION. THE NEAR 
INFRARED SPECTRUM OF MgO 


In order to be able to obtain better experi- 
mental results, a cubic crystal was sought whose 
fundamental corresponded to a wave-length in 
the near infrared where higher resolution could 
be employed. Consequently, the authors were 


pleased to learn that The Norton Company, 
using an electric furnace which operated at a 
temperature above 2500°C were able to fuse 
MgO, and prepare large flawless single crystals 
from the molten material. The samples obtained 
possessed almost every property which would 
make them suited for this purpose. They were 
clear, hard, easy to cleave (down to 0.075 mm), 
were not affected by water, were obtainable in 
sizes up to 3X41 cm, possessed almost perfect 
cleavage planes so that no polishing or grinding 
was necessary, were very tough and resistant to 
acids and heat, and the material had its funda- 
mental in the near infrared. Both transmission 
and reflection measurements were made through- 
out the region from 1y to 15.6u. Previous x-ray 
investigations have shown that MgO crystals 
are cubic in their structure being of the NaCl or 
face-centered type.® 


$1. Apparatus and procedure 

For this investigation in which it was desired 
to be able to study the absorption of the samples 
throughout the near infrared, the experimental 
set-up shown in Fig. 2 proved very satisfactory. 
The radiation from a hot platinum strip (S), 


’ P. P. Ewald and C. Hermann, Structurbericht (Leipzig, 
1931). 
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Fic. 3. Transmissions of cleaved plates of MgO. (a) 0.12 mm; (b) 0.47 mm; (c) 1.75 mm; 
(d) 3.05 mm. 


which was operated from a 22-volt transformer 
at a current of 20 amperes, was brought to a focus 
upon the sample by ./,, and after diverging was 
focused upon the slit of the spectrometer by the 
mirror >. After emerging from the spectrometer 
which contained a rocksalt prism and a Wads- 
worth mirror, this radiation was then concen- 
trated upon a Moll linear thermopile. The latter 
was connected as shown to an original Moll gal- 
vanometer (55 ohms internal resistance), the 
deflections of which were then amplified by 
means of the photoelectric relay? system. The 
filament F (4 volts, 4 amperes) was focused by 
the two lenses over the galvanometer mirror to 
the photoelectric cell, P, which was a selenium 
cell of the blocking layer type and had an active 
surface 3X3 cm. The use of the grids as shown in- 
creased the sensitivity threefold. An H.S. Leeds 
and Northrup 16.4-ohm galvanometer was then 
used to read the currents produced by the photo- 
electric cell. After taking extreme precautions to 
avoid any unsteadiness due to fluctuations in the 
temperature and pressure of the room, the relay 
system was set and used at a magnification of 
about 150-fold. The thermocouple, although not 
of the vacuum type, was placed in an evacuated 
container in order to avoid the pressure effects 


®*R. Bowling Barnes and F. Matossi, Zeits. f. Physik 76, 
24 (1932); L. Bergmann, Physik. Zeits. 32, 688 (1931). 


referred to just above. The scale distance was 
2.5 meters. 

The calibration of the instrument was deter- 
mined empirically, by mapping the positions on 
the scale of many well-known sharp bands such as 
those of H,O, the emission and absorption bands 
of COs, the reflection bands of quartz, the ab- 
sorption bands of benzene, toluene and cyclo- 
hexane. Points were thus obtained at rather close 
intervals out to 15.24. The transmission curves 
obtained were accurately reproducible as far out 
as the measurements were made. The authors 
realize that the extreme wave-lengths may be 
uncertain, due to errors in the calibration curve, 
by as much as +0.05u. The calibration curve was 
checked frequently, and aside from a shift back 
and forth from time to time of less than +0.01y 
due to changes in the temperature of the prism, 
was found to be constant. 

At the position marked A in Fig. 2, a sliding 
tower was built which enabled the sample to be 
moved in and out of the light beam, several 
measurements being taken in each of these two 
positions. The curves plotted below are in every 
case the actual percentage transmissions as 
measured, no corrections having been made for 
the losses of energy at the two surfaces of the 
crystal. The accuracy of the curves with regard 
to the ordinates is about +2 percent. 
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Fic. 4. (a) Transmission of a thin cleaved plate, £=0.075 mm. (b) Transmission of a 1y layer 
of MgO evaporated onto a nitrocellulose film. (c) Transmission of a layer of MgO fumed onto 


a plate of NaCl. 


§2. Results 

(a) Fig. 3 shows the results plotted to a yu 
scale for different thicknesses of MgO prepared 
by cleaving plates from one of the large single 
crystals. The similarity between the curves will 
be noted. The agreement of the secondary struc- 
ture in the various plates is good considering the 
resolving power of the instrument and the fact 
that as we proceed to greater thicknesses, the 
relative strengths and breadths of the absorption 
lines vary giving rise to absorption peaks of new 
shapes. 


In Fig. 4 curve (a) the composite curve result- 
ing from many separate runs on a cleaved plate 
0.075 mm is given. No experimental points are 
indicated but each of the minima shown has been 
checked in at least two or three of the individual 
sets of measurements. The structure is quite 
complicated, yielding some 40 absorption lines 
between 6u and 15.5y. 

(b) In the work of Czerny and Barnes’ thin 
films of the alkali halides were prepared by heat- 
ing the salts in a vacuum by means of a tungsten 
spiral and allowing them to condense onto a film 
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Fic. 5. (a) Retlection power of a crystal of MgO. (b) Dispersion curve of MgO calculated by 
use of the relation R=(n—1)?/(n+1)?. 


of nitrocellulose. It was shown in their paper that 
the optical constants of a NaCl film prepared in 
this manner were identical with those of a mas- 
sive plate of salt. Hilsch and Pohl'® and their co- 
workers have also obtained results which showed 
that films prepared in this manner were really 
single crystals. Further proof that such films are 
really crystalline has been obtained by Kirchner,"! 
using electron diffraction. 

Curve (b) in Fig. 4 shows the result of a series 
of measurements upon such a thin film of ap- 
proximately 1y thickness. It is immediately ob- 
vious that the wave-lengths of the fine structure 
lines agree in almost every detail with those of 
the composite curve for the thin cleavage plate 
shown above in curve (a). 

(c) In order to obtain still another check upon 
the existence of this secondary structure the 
transmissions of several thin films of finely 
powdered MgO were studied. The powders were 
prepared by two different methods, namely, by 
crinding up pieces of the large single crystals and 
by burning a piece of magnesium ribbon and al- 
lowing oxide to deposit upon a thin plate of NaCl. 
Curve (c) Fig. 4. shows the results of these 
measurements from which it is at once evident 
that the secondary structure is the same as that 


” R, Hilsch and R. W. Pohl, Zeits. f. Physik 57, 145 


(1929). See later papers. 
1 F, Kirchner, Ann. d. Physik 11, 741 (1932); Zeits. f. 


Physik 76, 576 (1932). 


found for the thin cleaved plate and the evapo- 
rated film. 

(d) The positions of the source and J, were 
now changed and an extra plane reflection put in 
at the focus A in order to be able to study the 
reflection power of one of the large crystals. At 
each wave-length the crystal was interchanged 
with the plane silvered mirror and the ratios of 
the deflections from the two taken as the percent 
reflection power of MgO. These results are shown 
by curve (a) of Fig. 5. A double reflection maxi- 
mum is seen at 14.84 and 15.34. From the values 
of R for the region below say 8u, the refractive . 
indices of MgO may be obtained by using the 
simple equation R= (n—1)*?/(m+1)*. The disper- 
sion curve thus arrived at is shown by curve (b). 
The value of n=1.86 at 1y is in fair agreement 
with the value y= 1.747 reported by Brice and 
Strong.” 


§3. Discussion of results 


(a) MgO has been the subject of several in- 
vestigations. Angstrém in 1889 studied its ab- 
sorption, while Coblentz in 1913 reported that 
with increasing wave-lengths its reflection power 
increased. In the latter investigation measure- 
ments were made at only 6 wave-lengths. In 
1928, Fraulein S. Tolksdorf made a study of 


= RF Brice and J. Strong, Phys. Rev. 47, 255 (1935). 
13S, Tolksdorf, Zeits. f. physik. Chemie 132, 161 (1928). 
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MgO from 2 to 154 by means of a prism spectrom- 
eter. Following a method suggested by herself and 
Friulein Laski,'* she studied the transmission of 
a thin film of finely powdered MgO held between 
two plates of NaCl. A transmission curve was 
obtained which ranged from about 40 percent at 
2u to 95 percent at 154 and which showed three 
distinct minima, located at 3.85, 7.65 and 14.2u. 
These she mistook for 4v, 2v and », respectively. 
In view of the predictions of the classical theories 
mentioned in the introduction, she called these 
vibrations. 


’ 


supposed overtones “‘ anharmonic’ 

Comparison of her results with those shown in 
Fig. 4, shows that Friiulein Tolksdorf had mis- 
taken the 14.24 secondary maximum for the 
fundamental, and had not measured any further 
out in the infrared. 

In 1931, Strong,'® using reststrahlen, showed 
that a layer of MgO fumed onto a NaCl plate 
showed a very strong absorption in the neighbor- 
hood of 234. A crystal of B-MgO showed 80 per- 
cent reflection at this same wave-length. He 
offered these results to show that the maximum 
found by Tolksdorf at 14.24 could not be the 
fundamental of MgO. 

In a short note, two of us reported the pre- 
liminary results of this investigation.'® At that 
time we stated that we had found a strong ab- 
sorption around 14u which we identified with the 
fundamental reported by Tolksdorf, the results of 
Strong not being known to us. In view, however, 
of the absorption band which we later observed 
at 15.44 and the work of Strong the identi- 
fication of any one absorption band as the funda- 
mental is quite uncertain. At present we are not 
in a position to say, for example, whether the 
reflection maxima of Fig. 5 are due to strong 
secondary maxima or to the fundamental. This 
question must be decided at a later date. From 
our present results one can only say that the 
fundamental absorption of MgO is accompanied 
by many secondary maxima. 

(b) The curves shown in Fig. 4 show very con- 
clusively that the infrared spectrum of a solid 
may be studied successfully in either of the three 


forms, namely, cleavage plates, powder or 


144 (5. Laski and S. Tolksdorf, Naturwiss. 14, 21 (1926). 

' J. Strong, Phys. Rev. 37, 1565 (1931). 

” R. Bowling Barnes and R. Robert Brattain, Phys. Rev. 
47, 416 (1935). 
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evaporated films. The last two methods require 
great care and sometimes lead to spurious results. 

Before using the powder method the investi- 
gator should be familiar with the authors’ paper 
and those of Pfund on the infrared of powders."” 
In the light of these papers, the increasing trans- 
parency of the powdered MgO found by Tolks- 
dorf as longer wave-lengths were reached, is 
readily understood. In all cases where powders 
are used the particle size of the powder must be 
smaller than the wave-lengths of the regions being 
studied in order to avoid spurious effects, and 
shifts in the wave-lengths. 

Hirsekorn'’ measured the transmission of films 
of NaCl evaporated at atmospheric pressure, and 
found spurious effects due to the particles. These 
evaporated films must be prepared in the very 
highest vacuum obtainable. 

These last-named methods have been used to 
advantage by many people, and it is therefore 
very gratifying to see that all three methods when 
applied to the same substance yield results 
strictly comparable in every detail. 

(c) In conclusion it must be mentioned that, 
since these secondary maxima are probably due 
to anharmonic terms in the force equation, they 
should be sensitive to changes of temperature, as 
will be discussed in Part B. Accordingly, prepara- 
tions are under way to repeat these measurements 
first at liquid air then at liquid hydrogen tem- 
peratures. Then and only then could a complete 
correlation with theory be expected. 


$4. Experiments on other crystals 


From time to time the complex structure of the 
infrared absorption bands of various crystals, 
as evidenced by reflection as well as transmission 
measurements has surprised experimenters. As 
increasingly higher resolving power has been em- 
ployed, it has become more and more evident 
that these bands possessed a real fine structure 
the exact interpretation of which is unknown. 

No attempt will be made at this point to give 
a complete list of papers bearing upon the above 
statements, however it is felt that several in- 
vestigations should be cited. 

7 R. Bowling Barnes, R. Robert Brattain and R. S. 
Firestone, Phys. Rev. 47, 792A (1935); A. H. Pfund, Phys. 
Rev. 36, 71 (1930); J. Opt. Soc. Am. 24, 143 (1934); J. Opt. 


Soc. Am. 23, 375 (1933). 
18H. G. Hirsekorn, Ann. d. Physik 6, 985 (1930). 
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In addition to MgO Tolksdorf studied CaO, 
BeO, ZnO, CdO and HgO. For the first three she 
found the spectra to consist of several distinct 
maxima of absorption where only one or two 
were expected. For CdO and HgO, she was un- 
able to locate the fundamentals because they 
were beyond 22u, the limit of her spectrometer; 
however, she did find in the near infrared, several 
weak absorptions which are probably analogous 
to the secondary maxima reported above for 
MgO. A complete account of this work, and an 
attempt to identify these secondary maxima as 
harmonics and combinations will be found in the 
book of Schaefer and Matossi, Das Ultrarote 
Spektrum. 

Silverman!® measured, in the neighborhood of 
the 7u band, the reflection and transmission of 
calcite, and found between 6.4u and 7.44 some 
33 lines in each case. This unexpected fine struc- 
ture he suggested was possibly due either to the 
disturbing effects of imperfections in the crystal, 
or to combinations with the far infrared bands of 
calcite. 

Silverman and Hardy’ showed that the 9u 
quartz band also possessed an extremely compli- 
cated structure. 

Recent papers from the laboratory of Schaefer 
and Matossi also point to the presence of compli- 
cated structure in the cases of many crystals. 

Besides the original work discussed above, 
other investigators”! have studied the alkali hal- 
ides and verified the existence of the secondary 
maxima in every case. 

Although the absorptions of all types of crys- 
tals seem to possess this fine structure, or sec- 
ondary structure, the most important thing to be 
done right now is undoubtedly to obtain very 
accurate measurements on cubic crystals. Ex- 
periments performed at low temperatures, with 
polarized radiation, or with the radiation travers- 
ing the crystal in directions other than normal to 
the 100 plane, will yield results which will be of 
the greatest theoretical importance with regard 
to the study of the solid state. 


19S. Silverman, Phys. Rev. 39, 72 (1932). 

20 J. D. Hardy and S. Silverman, Phys. Rev. 37, 176 
(1931); S. Silverman, Phys. Rev. 45, 158 (1933). 

21K. Korth, Nachr. d. Ges. d. Wiss. zu Gétt. Math- 
Phys. KI. 1932, p. 576; A. Mentzel, Zeits. f. Physik 88, 178 
(1934). 
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PaRT B. THEORETICAL INTERPRETATION 


$1. Previous developments: a criticism of the 
Born approximation 


A sketch of the general scheme followed in an 
attempt to approach a theoretical interpretation 
of the previous results was presented in the intro- 
duction, and before proceeding with a detailed 
discussion of this, we shall pause to review some 
of the previous developments. 

In the earliest period of observation of the 
reflection spectra of solids, it was found® that 
for cubic crystals the spectra could be accounted 
for by assuming that they possessed one funda- 
mental absorption frequency and that the reflec- 
tion properties were related to this by means of 
the well-known formulae of optics. In noncubic 
crystals, it was necessary to assume the presence 
of several absorption maxima. Following this, a 
number of dimensional equations relating this 
frequency for cubic crystals to other physical 
properties were suggested by Einstein, Linde- 
mann, Braunbeck and others, and were found to 
give rough agreement with the established facts. 

On the basis of some remarkably good assump- 
tions concerning the arrangements of atoms in the 
alkali halides, Madelung*’ succeeded in relating 
the fundamental absorption frequency to the 
interaction forces and connected these with the 
elastic coefficients. This was followed by a paper 
of Born and von Karman” in which the ideas be- 
hind Madelung’s work were placed upon a firmer 
basis for analytical treatment. Although the con- 
clusions which they drew from their development 
are not strictly correct when viewed in the light 
of modern knowledge, their treatment of vibra- 
tional waves has continued to be of service in 
many fields. In addition, this work gave an indi- 
cation of the weaknesses of the Debye theory of 
specific heats in which the crystal was treated as 
a continuous solid—a question which has been 
the object of considerable qualitative and semi- 
quantitative investigation in recent times.*° 

22 A complete sketch of the work discussed in this para- 
graph is to. be found in Schaeffer and Matossi, Das 
Ultrarote Spektrum (Berlin, 1930). 

23 E. Madelung, Nachr. d. Ges. d. Wiss. zu Gétt. M. P. 
Klasse, March 1909, Jan. 1910; Physik. Zeits. 11, 898 
(1910). 

24M. Born and Th. von Karman, Physik. Zeits. 13, 297 
(1912). 


* M. Blackman, Proc. Roy. Soc. A148, 365, 389 (1935); 
A149, 117, 126 (1935). 
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During a subsequent period, Born and co- 
workers succeeded in developing a theory of 
ionic crystals,?* based upon the assumption of 
central force interaction of ions, by means of 
which a number of properties of such crystals 
such as dielectric constant and fundamental 
frequency could be correlated. Such correlations 
are to be regarded as first approximation connec- 
tions which are valid when the assumptions made 
are correct to within factors that are essentially 
refinements. The principal assumption which is 
made is that of central force interaction of ions. 
For the purpose of many property-correlations, 
such as those mentioned previously, considerable 
inaccuracy in these assumptions may not lead to 
wide discrepancies in results, simply because the 
quantities correlated are not sensitive to alter- 
ations in the fundamental assumptions; that is, 
for example, neither the dielectric constant nor 
the fundamental frequency of a crystal will be 
greatly affected by small deviations of the inter- 
action potential from one of the central type. This 
is not true of many other properties of crystals, 
however, particularly those which are the differ- 
ences of relatively insensitive quantities. This 
question is of considerable importance when the 
problem of determining the positions of the 
secondary absorption maxima of crystals relative 
to the fundamental is raised, as will be shown in 
the following sections, and for this reason it is 
important that the extent to which the results of 
the assumptions of Born are valid be known. If 
they are found to be incorrect to a sufficient 
degree, it will be necessary to proceed from an 
entirely phenomenological standpoint, as will be 
shown to be the case. 

One of the principal results of the assumption 
of central forces when dealing with arbitrary 
crystals is the so-called Cauchy-Poisson rela- 
tions’? between the components of the stress- 
strain tensor. For cubic crystals, with which we 
are principally concerned, these relations are 
such that equality is to be expected between the 
coefficients C2. and cq of Voigt, the first of which 
connects the transverse stress to' the strain along 
a given crystallographic axis, while the second 


26 A thorough review of this work is given by M. Born 
and M. Goeppert- Mayer in the Geiger-Scheel Handbuch der 
Phystk, Vol. 24b. 

27 W. Voigt, Lehrbuch der Kristallphysik (Teubner, 1910), 
p. 607, 
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TABLE I. 
e107! 121078 ca10-"8 

°K NaCl MgO NaCl MgO NaCl MgO 
273.1 5.035 28.841 1.277 15.350 1.279 8.705 
210 5.284 29.128 1.298 15.421 1.253 8.585 
150 5.508 29.396 1.315 15.472 1.127 8.535 
80 5.754 —_ 1.332 —_— 1.159 _- 

77.4 29.547 - 15.501 — 8.462 


relates the shearing stress in a plane orthogonal 
to one of the principal axes to a shearing strain in 
the same plane. A number of results concerning 
these relations for valence crystals and metals** 
revealed that they are not nearly satisfied. Previ- 
ous*® to some very recent work by Quimby, 
Balamuth, Rose and Durand’ using oscillating 
circuit methods, the only measurements carried 
on for the determination of the constants for 
ionic crystals dealt with the monovalent alkali 
halides and were performed at ordinary tempera- 
tures. In these earlier cases the relationships were 
found to be valid to within an experimental error 
which we judge to be about 5 percent. The room 
temperature values are not particularly valuable 
for the present discussion, however, since the 
elastic coefficients are known to vary consider- 
ably with temperature and we are really inter- 
ested in the absolute zero values. The work of 
reference 7, however, has furnished a good deal of 
information for both monovalent (NaCl) and di- 
valent ionic crystals (MgO) over a temperature 
range extending down to 80 degrees Kelvin with 
an accuracy exceeding that previously obtained 
by about one power of 10. A sample set of values 
of the adiabatic coefficients of both of these 
crystals are given in Table I. The isothermal 
values may be determined from these by ade- 
quate small corrections which are always such as 
to lower cy, and leave ¢2. unchanged. 

From these results it is clear that in the vicinity 
of absolute zero the cy; coefficient will be about 
15 percent less than cy. for NaCl, while in the 
case of MgO the ratio is almost 2 to 1. This 
presents conclusive evidence that the Born ap- 
proximation is a rather rough one for ionic 
crystals. 

28 A review of the work of Voigt on this subject may be 
found in M. Born’s Atomtheorie des Festenzustands 
(Teubner, 1923), Chapter I. 

29 W. Voigt, reference 26, p. 741; P. W. Bridgman, Proc. 
Am. Acad. 64, 19 (1929). We are indebted to Professor 


Bridgman for correspondence concerning details of this 
work. 











INFRARED ABSORPTION 

With the discovery of secondary maxima, 
Born and Blackman,* as we have stated in the 
introduction, sought an explanation of this phe- 
nomena on the basis of anharmonic potential 
terms and solved the problem of the diatomic 
linear chain, considering such terms to be small 
perturbations, as a classical mechanical one. As 
is well known, the frequency spectrum for this 
problem as a function of the variable o=27/, 
\ being the wave-length, is represented by a curve 
of the type shown in Fig. 6 in which the discon- 
tinuity occurring for «= +7/Na, where a is the 
distance between nearest like atoms, is a charac- 
teristic of the diatomic lattice that disappears 
when the masses of the particles are equal and 2 
is the total number of atoms. The results of Born 
and Blackman show that on the basis of their 
classical reasoning one would expect frequencies 
vi +ve and v;—v2 as well as the fundamental fre- 
quency v to be absorbing. These results were then 
extended to the three-dimensional case again on a 
classical basis, but with neglect of a suitable 
analysis of the discontinuous zone structure that 
must exist by analogy with the results of the 
theory of metals. This last point represents a 
serious omission, even if the use of classical 
mechanics and the use of central forces were 
valid, and will be discussed more thoroughly 
in §3. 

In view of these facts, the problem will be 
handled entirely anew in the following sections 
and with a view toward eventually obtaining 
correlations between experimental and theoret- 
ical results. 
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Fic. 6. A representative energy-wave number distribu- 
tion for the linear chain. In the classical case the energy is 
proportional to the square of the frequency, while in the 
quantum case it is proportional to it. In the three-dimen- 
sional lattice, this is characteristic of the energy spectrum 
for one polarization and a line passing through the origin in 
o-space. The discontinuities occur when this line intersects 
the zones. 


SPECTRA OF CRYSTALS 591 
$2. The symmetry restrictions induced on a 
general potential form 


As a basis for a phenomenological discussion of 
the interaction potential of atoms, we shall as- 
sume that they are regularly arranged in a lattice 
at points 


Ir: +mrt2+nT3+ Pa (1) 


where /, m and nm are integers, 7, 72, T3 are the 
primitive translations of the lattice, and p, is a 
vector independent of /mn which gives the 
position of the ath atom in a unit cell. For a 
diatomic face-centered lattice of the cubic type, 
T), Ts, T3 May be taken in the form 


a a 0 


a}, 0}, a}, (2) 


0 a a 


respectively,*® while p, may be chosen to be zero 
for one atom (type a, say) and 


a 
0 (3) 
0) 


for the other (type 0), so that a is the closest dis- 
tance between unlike atoms. Thus when the 
atoms are in their equilibrium positions, the 
lattice will be invariant under one of the 230 
space groups and we must expect the equilibrium 
potential function to be likewise invariant (i.e., 
to belong to the unit representation of the space 
group). When the atoms are displaced from their 
equilibrium positions, so that the new position of 
an atom of type a, say, in the /mn-th cell relative 
to its old is given by the vector of component 
xa‘(lmn)(t=1, 2, 3), the potential of the lattice 
will be altered, but we shall assume that it may 
be expressed as a power series expansion in these 
displacement variables in such a manner that all 
terms of a given degree are invariant under the 
group. As is well known, the constant term in this 
expansion need not be considered, while the linear 
terms will vanish as the condition of equilibrium. 
We shall regard the quadratic terms as being of 
principal importance and treat the cubic terms 
as small perturbations. For this purpose, we shall 


8° For a presentation of the theory of crystal symmetry in 
terms of the notation to be presented here see F. Seitz, 
Zeits. f. Krist. 88, 433 (1934), et seq. 
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begin by neglecting the latter and choose a set 
of normal coordinates for which the quadratic is 
in normal form, which must be taken to be real 
since we will want to quantize the entire system. 
We shall consider the crystal to be fixed, so that 
the process of choosing normal coordinates is es- 
sentially one of resolving all of the x’s into 
Fourier-like components and treating the ampli- 
tude-coefficients as a set of variables. Since one 
of these will correspond to pure translations, it is 
important that the potential be invariant under 
infinitesimal translations, and we must express it 
as a function of differences 

xa'(l'm'n’) —x,'(lmn)=&,-‘(l'm'n', lmn) (4) 


of displacements, which will also allow us to ex- 
press the potential change going with a homo- 
geneous deformation of the lattice. It is implied 
that the coefficients in the power series expansion 
are so defined that the entire expression is in- 
variant under infinitesimal rotations because we 
are dealing with a scalar form, so that we shall be 
concerned with all quadratic and cubic forms in 
the differences that belong to the unit representa- 
tion of the space group. Of these a very large 
number may be neglected since they are con- 
sidered to be physically unimportant, and we will 
actually consider but a very small fraction of the 
possible terms. 

For simple crystals in which each atom is at a 
center of symmetry of the type for which the 
corresponding point group is that of the class to 
which the space group belongs, the problem of 
finding the desired form is comparatively easy. 
In this case the translational symmetry may be 
handled completely apart from the rotational, 
and we need only take the sum of the equivalent 
potentials of equivalent cells, these potentials 
being determined by mere rotation group con- 
siderations (i.e., screw operations and glide planes 
do not play a role). 

In the lattice of the NaCl type (Fig. 7) with 
which we shall be entirely concerned, each atom 
possesses the symmetry O” for which the ten 
irreducible representations have been given by 
Bethe.*' The character system going with these is 
listed in Table II, in which [’* and I~ are written 
together with the corresponding signs for the 
elements of determinant 1 and —1. 


"31H, Bethe, Ann. d. Physik 3, 133 (1929). 
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Fic. 7. Face-centered cubic lattice of the NaCl type. 


Although the translation group is of the type 
r.’, it is convenient to designate each atom by the 
integers /mn which specify its position vector of 
components /a, ma, na, since these lie at all 
mesh points of a simple cubic lattice. Of the pos- 
sible interaction terms of each atom with its 
neighbors, we shall consider only those between a 
given one and its twenty-six neighbors. Moreover, 
of these we shall be concerned only with quad- 
ratic terms of the type 


E'(l'm'n', lmn)£*(l'm'n', lmn) (5) 
and with cubic terms of the type 
E*(l'm'n’, lmn)é*(l'm'n', lmn)&*(l'm'n'’, lmn) (Sa) 
and discard those in which the indices of the 


TABLE IT. Character system for crystallographic group O*. 











E 6C2’ 8C3s 3C2 6C4 I 6Ra 8S6 3R, 654 
r 1 1 1 1 1 +1 +1 +1 +1 +1 
re 1 1 —-1 —1 1 +1 +1 Fl +1 +1 
rs 2 2 0 0 -1 +2 +2 0 0 +1 
im 3 —1 1 —1 0 +3 ¥1 +1 ¥1 0 
rs 3 -1 —1 1 0 +3 F1 ¥1 +1 0 








where: E is the unit element; 


6C? is the class of 6 two-fold rotations about axes of the (110) 


ty 


rpe; 
8Cs3 is the class of 8 three-fold rotations about axes of the (111) 


type; 


3C2 is the class of 3 two-fold rotations about axes of the (1Q0) 


type; 
6C; is the class of 6 four-fold rotations about axes of the (100) 


type; 
J is the inversion; 


6Rgq is the class of 6 reflection planes orthogonal to axes of (110) 


type; 
8Seis the class of 6 six-fold rotary reflections about axes of (111) 


type; 


3R, is the class of 3 planes orthogonal to axes of the (100) type; 
6S< is the class of 6 rotary reflections about axes of the (100) 


type. 
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type /'m'n’ in these differ among themselves. 
This simplification has the effect of reducing the 
total number of terms for the twenty-six neigh- 
bors to about a fourth. Of the terms considered, 
those which describe the interaction between un- 
like atoms are the same regardless of which type 
of atom is at the center of neighbors, and need to 
be counted but once for each cell, while those 
which describe the interaction of like particles 
appear with different coefficients depending upon 
the type of atom concerned. 

The problem of analyzing the character system 
of the quadratic and cubic terms will not be gone 
into here, and it will suffice to tabulate the 
results. Regarding the former, it is found that 
there are ten forms to be considered as these 
may be classed as follows. In all cases the vari- 
ables &, n, ¢ will be considered in place of £', &, &, 
and because of the restriction (5) we may 
write &(l'm'n’, lmn) in place of &(l'm'n’, lmn) 
-n(l'm'n', lmn), etc. The summation indicated will 
always imply that two terms in which a cyclic 
permutation of the variables that appear are to 
be added. The character system of the various 
cases is listed in Table ITI. 

I Interaction of given atom with 6 unlike neighbors 
1. Sfe2(d+1mn)+#2(/—1mn)}; 


2. Lse(mn+1)+e(/mn—1) 
+£2(/m+i1n)+2(dm—1n)}. 


Il Interaction of given atom with 12 nearest like neighbors 


3. T{E(m+1n4+1)+e(lm+i1n—1) 
+2(lm—1n—1)+22(lm—1n+1)}; 


4. Lf{e(+1m+i1n)+8(/+1m—1n) 
#2(/+1mn+1)+22(1+1mn—1)}; 


5. Sfem(lt+im+in)+e¢(1+1mn+1) 
—in(l+im—i1n)+é¢(1+1mn—1)}; 


6, 7 and 8 duplicate 3, 4, 5 for atoms of other type. 
III Next nearest unlike neighbors 
9. S{Le(d+41, m+1, n+1); 
+1 
10. © {(n+OeE¢+1m+1n4+1)—(n—H)E(l+1m—1n+1) 
—(n+Oe(l+1m—1n—1) 
+(n+He(1+1m+i1n—1); 
—(nt+o)i(l—1m+1m+1)+(—HEl—1m—1n—1) 
+(n+e)e(l—1m—1n—1) 
—(n—S)E(l—1m+i1in+1)}. 


In the following we shall designate the coefficient of 
each of the ten quadratic terms by a(1), ---, a(10), 
respectively. 


Since we shall employ the cubic terms only 


TABLE III. Character system for quadratic terms considered. 








TYPE 
OF ConpDI- 
NEIGH- TIONS 
BOR (5) E 6C2’ 8C3 3C2 6C4 I GRa 8Se 3R, OSs No. 
I i=j 18 0 0 6 2 0 2 0 12 0 2 
I ix) 18 0 0-2 -—2 0 2 0 -4 0 0 
II i=j 36 2 0 0 0 0 2 0 12 0 2 
II ij 36 2 0 0 0 0 2 0 -4 0 1 
III [=j 24 0 0 0 0 0 4 0 0 0 1 
Ill ixj 24 0 0 0 0 0 4 0 0 0 1 











insofar as they are able to break down the selec- 
tion rules for optical transitions, and we shall as- 
sume that they do not affect the energy spectrum 
given by the quadratic terms, it will not be neces- 
sary to go farther than to write down the set of 
terms arising from (5a) when only the six nearest 
neighbors are considered for reasons that will 
become apparent in $5. The character system is 
given in Table IV. 

The corresponding cubic forms are 


1. >> fe8(/-+1mn) — (l—1mn)} 
2. >} £*n(lm+1n) — E2n(lm—1n) 


&2¢(/mn+1)—#(/mn—1)}. 


TABLE IV. Character system for cubic terms considered. 


RELATIONS IN 


(Sa) FE 6C2’ 8Cs 3C2 6Cy I GRa &8Se 3Ra OSs 


i=j=k 18 —2 2 0 2 0 4 0 1 
ixj=k 36 0 0 —4 0 0 0 0 & 0 1 
ixj#k 6 0 0 2-2 0 2 0 4 0 0 


$3. Normal coordinates: the energy spectrum 


The process of reducing the quadratic expres- 
sion, corresponding to the sum of the possible 
invariant terms introduced in §2, to normal form 
is essentially one of reducing a matrix to diagonal 
form, and since this matrix is invariant under an 
entire space group, we know that the normal 
coordinate system will be one in which the space 
group is reduced. When the Born-von Karman 
boundary conditions are employed, it may be 
shown™ that for all crystals these coordinates 
may be taken to be of the Fourier type 

Xo exp (ior) (6) 
where x. is invariant under the primitive trans- 
lations of the lattice and o is a vector such that 

*F. Bloch, Zeits. f. Physik 52, 555 (1929); F. Seitz, 


On the Reduction of Space Groups to appear in Annals of 
Mathematics, January (1936). 
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o-T,=2rn,. (7) 
In this, T; are the three vectors defining the 
polyhedron at the surface of which the boundary 
conditions are satisfied. For a simple cubic lat- 
tice of type T’., for example, this polyhedron may 
be taken to be a cube so that these vectors lie 
along each of the three axes. In the case at hand 
in which we are dealing with lattice vibrations, 
x. will be periodic in the primitive translations 
and will be a vector function which takes values 
only at the positions of the atoms. In other 
words, for a lattice of the type considered here, 
the functions (6) will represent waves in which 
the atom of type 1 at Ja, ma, na has the dis- 
placement 


xo(1, /mn)fi(o, mn) exp [io-r(lmn)] (8) 


where x, is the scalar amplitude of the wave, f; is 
the unit vector in the direction of displacement 
and r(/mn) is the position vector of the atom, 
while atoms of type 2 at the other mesh points 
will have displacements of scalar amplitude x(2) 
and direction f,, both of which will generally 
differ from x(1) and f;. For an arbitrary, allow- 
able o, the three values of f; and f. which reduce 
the terms of type o to normal form may be deter- 
mined by solving an appropriate secular equation 
to be described later. 

For a face-centered lattice in which the primi- 
tive translations are given by (2), the allowable 
form of a is 


2 Ni +N2—N3 
T 
—n+N2+ns3 (9) 
2Na 
Ni —N2+N3 


corresponding to the boundary conditions being 
satisfied at the surface of the rhombic dodeca- 
hedron of height 2Na (Fig. 8), where m, m2, m3 
and N are integers. The lattice defined by allow- 
ing , m2 and m; to take on all integer values is 
the inverse lattice to the face-centered one and is 
seen to be a body-centered one having a unit cell 
of volume 47/N*a*. Not all values of the 7's cor- 
respond to independent modes of vibration, 
however, and it is readily seen that the only ones 
which need to be considered are those contained 
in a cube whose corners are at the eight points 
corresponding to m, m2, m3; equal to +N/2. 
Within this there are 2N* unit cells so that when 
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Fic. 8. The rhombic dodecahedron at the surface of 
which the boundary conditions are satisfied. It possesses 
one-fourth of the volume of the enclosing cube. 


the three directions of polarization are taken into 
account, it corresponds to 6N* degrees of free- 
dom, which is in agreement with the fact that 
there are 2.V* atoms in the rhombic dedocahedron 
containing N* unit cells of the crystals, and each 
atom possesses three degrees of freedom. 

From the standpoint of the energy spectrum, 
we may expect three energies going with each 
point in the cube described above corresponding 
to the generally different energies going with the 
three directions of polarization of a wave of given 
a. It is to be expected that these energy functions 
E*(c), where 7 ranges over the three polarization 
vectors, will generally vary smoothly from point 
to point except at the surfaces of definite poly- 
hedra, or zones, where they will have discontinui- 
ties analogous to those met with* in describing 
electron waves in solids. The space between two 
zones will generally contain as many points of o 
space as there are unit cells in the polyhedron at 
the surface of which the boundary conditions are 
to be satisfied—N* in this case, so that we may 
expect that there will be two zones in the prob- 
lem at hand, and, in fact, for the lattice I.’, 
these correspond to the cube and truncated octa- 
hedron shown in Fig. 9. Since we do not consider 
regions outside of this cube to be different from 
that inside, the energy surfaces will end at the 


%8 L. Brillouin, Quantenstatistik (Berlin, 1931), Chap. VIII. 

The basis for examining the connection between these 
zones and crystal symmetry is laid in the second of the 
references 32, and will be extended for presentation in this 
journal in the near future. 
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Fic. 9. First two zones for a face-centered diatomic 
crystal of the NaCl tvpe (translation group I’). The inner 
zone is a truncated octahedron and the second is a cube. 
The difference in volume between the first and second is 
equal to that of the first. The numbers indicate particular 
points to be discussed in the text. (Note no. 2 is the center.) 


cube, but may be expected to possess a vanishing 
gradient normal to the surface,*4 while at the 
surface of the truncated octahedron they will 
possess definite discontinuities, being continuous 
on each side of the surface, and will have a 
vanishing normal gradient there. Fig. 6 represents 
the characteristic behavior of E(c) for one direc- 
tion of polarization as ¢ varies along a line passing 
through the origin. The discontinuity occurs 
when o reaches the truncated octahedron, and 
the curve ends at the cube. 

It is a general characteristic of this zone struc- 
ture that each of the points in one zone may be 
viewed as being associated with a definite point 
in another zone in the sense that the waves going 
with the two have exactly the same symmetry 
properties (belong to the same irreducible repre- 
sentation), and in the present case these two 
points are connected by a vector of the same 
length as one of the vectors 


+1 
2n 
s'=—{ +1 (10) 
2a 
+1 


A proof of this may be derived from the work of H. 
Jones and C. Zener, Proc. Roy. Soc. A144, 101 (1934) by 
means of group-theoretical extensions that will be dis- 
cussed in this journal in the near future. 


joining the origin with a cube corner. These 
vectors are the only members of the 2N® con- 
sidered (aside from the origin), that possess the 
property that their scalar product with each of 
the primitive translations is 27, or 0. It is readily 
seen that to each point between the cube and the 
truncated octahedron, the vector ¢+s' lies either 
outside of the cube (and is not to be considered) 
or inside of the truncated octahedron, the only 
ambiguous cases being those in which @ lies on 
the surface of the latter polyhedron, and this 
difficulty, which is never of physical importance, 
may be removed by assuming that one of any 
two such points that are connected by a vector s 
belongs to the inner zone and the other to the 
outer. Because of this fact and the fact that 
a+s' possesses the same properties as o under the 
translation group, we need only consider points 
belonging to the truncated octahedron, and it will 
turn out that the correct multiplicity of solutions 
corresponding to the double-valuedness of each 
point inside will automatically follow. These two 
solutions are, in fact, the modes of vibrations in 
which neighboring atoms are out of phase by 
more or less than 90 degrees, respectively. Later 
on, we shall return to the characterization by use 
of all points in the cube, and consider the three 
energies going with the wave in which unlike 
neighboring atoms are most out of phase to be 
outside of the truncated octahedron. 

Finally, it must be borne in mind that we want 
all waves to be real for the purposes of quantiza- 
tion, and for this reason, we shall replace 
exp (i0-r) in the preceding by both sin o-r and 
cos ¢-r, and have for the Fourier resolution of 
the displacement vector x‘(/mn) of atoms of type 
a and b: 

(1) Atoms of type 1 


. Xa(o, ¢) 
x‘(lmn)= >> 4- —fa'(o, c) cos o-r(lmn) 
" S VM) 


Xa(o, S) 


——f,'(o, s) sin er (11) 
Vm 


where xa(¢, c) is the amplitude of the cosine wave 
polarized in the direction f,‘ (a= 1, 2, 3), etc. The 
allowed vectors f' are to be determined by an 
appropriate set of equations. 


oes eee 
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(2) Atoms of other type 


ei , , , ia . Valo, c) ’ ‘a . a , , , 
xi(l'm'n’)=>> ga'(o,¢) cos o-r(l'm'n’) 
- V Me 
Yala, S) | 
+ g'(o, Ss) sin o-r(l'm'n’) (lla) 
\/ Mle | 


where y and g have the same significance as x 
and f in (1) and ¥, accompanies Xa. 

The introduction of the factor 1/\/m, and 
1/\/me has the effect of reducing the kinetic 
energy form to one in which x2 and ¥ have the 
same coefficient and it is only necessary to con- 
sider the potential form in determining the 
normal coordinates. 

It may be readily shown that in all cases in 
which the atoms occupy points of inversional 
symmetry: 

—1 
—1 
—1/, 


the sine and cosine waves separate completely 
in all of the sums that will be taken so that we 
need only consider the latter and designate 
x(o,c) and y(o,c) by x(c) and y(c) for the 
present. -In considering the cubic terms, it will 
be necessary to return to the previous notation. 
From this it is readily seen that for a given o 
the six quantities xa(c), ¥3(o) (a, 8=1, 2, 3) will 
give rise to a secular equation of the sixth degree 
unless one or more of the polarization pairs is 
known beforehand, in which case the secular 
determinant may be reduced to a two- and a 
four-, or to three two-dimensional ones, depend- 
ing upon whether one or two pairs is known. If 
they are not known, the solution of the secular 
equation will determine these directions in terms 
of the set of f, and g3; that have been assumed at 
the start in (11) and (11a). It is clear from the 
symmetry of the quadratic form that the f’s 
will be mutually orthogonal, as will be the g. 


$4. Wave quantization and selection rules 


The process of quantizing the waves once the 
reduction to normal form has taken place is 
simply one of writing the Schrédinger equation 
corresponding to the Hamiltonian function 
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LE (Plo, ©) +4x2(0)G?(o, €)) 


o a 


+ (q7(o, 5) +K,?(o)q?(o, s))} (12) 


where it assumed that the g,’s are the amplitudes 
of the normal variables of given o so that the 
total energy is the sum over o and X of terms of 
the type (12), the former taking on all values 
corresponding to the cube. The x°(¢) are the 
roots of the secular equation discussed in the last 
section. The Schrédinger equation for this prob- 
lem is simply 

h? d*W(qr(o), m(o))  Ky7*(o) 


_ +—— qr*(a) ¥(qrlo), my(o)) 
8x" dgx*(a) 2 


— E(ny(o))¥qa(o), m(o))=0 (13) 


for which the energy of the stationary state 
¥(gr(o), m(0))(m(c) =0, 1, 2, ---) is given by 
h 


—(n+1 2)x,(o), 
2r 


(13a) 


k,x(7) being the classical oscillator frequency. 

The probability of optical transitions from a 
state y, to y, will be proportional to the square 
of the dipole moment matrix component con- 
necting these two states, and this matrix will be 
zero for all ¢ except o=s' in which case the three 
waves of equal energy corresponding to a normal 
mode in which unlike neighbors are 180 degrees 
out of phase (i.e., the mode going with ¢=s' of 
the last section which we shall designate simply 
by s in the future) possess a nonvanishing ma- 
trix. This matrix p,(s) (7 =1, 2, 3) will be of the 
form 


p,(s)=e(s)q-(s)/V/ VW (14) 


where J =m ,me2/(m,+ me) and e(s‘) is the effec- 
tive ionic charge to be determined from precise 
theoretical considerations. If the ions were rigid 
spherical structures, which the evidence of $1 
indicates to be only roughly true, this would be 
the absolute value of the electronic charge for the 
alkali halides and twice this for MgO. In the 
present paper, however, we are interested only in 
the possibility of absorption and not in the in- 
tensity of the absorption lines so that this point 


% C.J. Brester, Kristallsymmetrie und Reststrahlen, Diss. 
(Utrecht, 1923). 
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is not a critical one. Since p,(S) is proportional to 
q-(s) it will have nonvanishing matrix compo- 
nents only if z’=n+1, so that only light quanta 
of frequency »(s) may be emitted or absorbed in 
the quadratic approximation. This is of course 
just the classical result as was to be expected. 


$5. Alterations induced by the presence of the 
cubic terms 

As we have stated previously, we shall assume 
that the results of the presence of the cubic terms 
will be such as to leave the energy spectrum un- 
altered, but will affect the wave functions enough 
so that the selection rules of §4 are no longer 
strictly valid. These cubic terms will be of the 
form 


Ho LD ater EW)}, (15) 
l d 


ijk wire 


where /, /’, /’’, l'’ designate triples /mn, etc. and 
the al!’ are the coefficients upon which the 
proper symmetry restrictions have been placed. 
If the values of the x's are substituted, it is 
readily seen that this reduces to the form 
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Dd Basy(o10203)galoi)ga(o2)g,(o3), (15a) 


apy 
in which 0), o2, o; in a given term are subject to 
the condition 
o:toeto;=s’, (16) 
where s’ is any vector such that 
s’-r'=2nn;, (16a) 


7; being one of the primitive translations (2) and 
n; an integer. There are rather important restric- 
tions upon the §’s as we shall see later, but we 
shall not be interested in these at the present 
time. 

If we now add the sum of terms (15) to the 
Hamiltonian function, we may obtain the new 
stationary states by use of an appropriate pertur- 
bation scheme such as that of Schrédinger or 
Brillouin-Wigner.* The unperturbed eigenfunc- 
tions of the system will be of the form 


W.= Il,, aW(Gala), Na(c)) (17) 


corresponding to a state in which the coordinate 
Ja(o) is in the m.(o)th quantum state. The corre- 
sponding perturbed function will be of the form 


W(--+, gilo), -e+3 +++ mi(o)) = MaeW(Galc), Mala)) +  ¥ y(n'a(o1), n'3(o2)n',(03)) 


XVv(galor), Na (o1))¥(ga(ee), n3'(o 


in which o, o2, o3 satisfy the rule (16), II’ does 
not include the wave functions of the variables 
Ya(o1), ga(o2), gy(o3) and 


Na’ (o1)= Na(o1) +1, 
ng'(o2)=Ng(o2)+1, (19) 
n,'(o3)=n,(o3) +1. 


The coefficients y are proportional to matrix 
components of the cubic perturbation term (15), 
which we need only assume to be small. The 
conditions (19) follow from the harmonic oscilla- 
tor nature of the unperturbed problem and it is 
to be understood that n’ cannot be negative. In 
order to find the transition probability two states 
of the type (18), it is necessary to examine the 
matrix components of p(s) given by (14) since, 
as was pointed out there, this represents the 
dipole moment for the entire crystal, and from 


919203 apy 


2) )¥(gy(e3), n,' (a3) ) . in a¥(ga(a); n,'(c)), (18) 


the remarks made in §4 it follows that v(s) will 
still be the principal frequency because VW» is the 
prominent term in V. Because of the presence of 
the higher order terms in (18), this may possess 
nonvanishing dipole moment components with 
other states WV’ arising from nonvanishing com- 
ponents between the first-order term of ¥ and 
the second-order terms of ¥’ and vice versa. 

There are several prominent cases that we 
shall handle separately: 

Case I. Let us assume that in WV the various 
waves of ga(a) are in states given by the quantum 
numbers n,(¢). Then among the higher order 
terms, which will be a linear combination of 
other states, there will be characteristic terms in 
which the states corresponding to ga(o1), gs(o2), 


% The Schriédinger perturbation scheme may be found in 
any textbook on quantum mechanics. L. Brillouin, J. de 
phys. et rad. (III) 1, 373 (1932); E. P. Wigner Bull., 
Hung. Acad. (to appear shortly). 
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Yy(os), Ys(S) are given by 
Mal(oi)t1, Nglo2)+1, my(o3)+1, mals), (20) 


where o;, a2, 3; must satisfy the condition (16). 
There will clearly be 2* such states, unless some 
of the n’s are zero in which case only the positive 
sign may appear. These higher order terms will 
possess important nonvanishing components only 
with that state W’ in which the zero-order term 
has quantum numbers ,(¢), excluding those 
going with ga(o:), gs(o2), gy(o3), gs(S) in which 
cases the quantum numbers must be 


Naloi)+1, mal(oo) +1, 


Na(o3)+1, m,(S)+1, (21) 


in which the signs of the first three terms must 
agree with those of (20), so that the energy dif- 
ference between V and W’ is 


h(+va(o1)+v3(o2)+v,(03)+(S)), (22) 


that is, all energies corresponding to the sum (22) 
may be absorbed or emitted, the first case being 
that in which the resultant quantity is positive 
and the second that in which it is negative. If the 
matrix components between the zero-order term 
of W and the higher order terms of WV’ are con- 
sidered, it is readily seen that they introduce no 
new possibilities, while the terms arising from 
interconnection of the higher order term of both 
W and W’ will be cast aside as negligible. 

Case II. If o; (or either of the other two oa’s) 
is equal to s, the argument of case I is practically 
repeatable, the only difference being that but 
three waves qs(e2), ¢g,(¢3), gs(S) need be con- 
sidered. The energies that may be absorbed or 
emitted are now given by 


(23) 
(23a) 


+ vg(o2)+v,(a3) 
or 
+ v(o2)+v,(03)+2v(s), 


which are just the results one would obtain from 
(22) by setting v.(o1) =v(s). 

Case ITI. If either of the o’s were equal to zero, 
difficulty would arise from the fact that the cor- 
responding states belong to a continuous spec- 
trum, and a special treatment would be required, 
but we may safely assume that the absorption 
involved in any single transition of this type is 
negligible and that the cases in which @ is prac- 
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tically, but not quite zero (see below) are the 
really important ones. For these, the discussion 
given above is valid, and the frequencies absorbed 
or emitted are practically 


+va(o2)+vg(o3)+r(S). (23b) 


The additional cases in which all of the o's 
are equal to s or zero do not occur as will be 
seen later. 

It is perhaps wise to repeat that the negative 
signs in (22), (23), (23a) are valid only if the 
corresponding states are excited in VW, since this 
fact will give rise to much of the temperature 
dependent effects; that is, at sufficiently low 
temperatures only the positive signs will be 
important. 

The absorption (or emission) coefficient for any 
given frequency of the allowed type will neces- 
sarily be small if the assumptions made here are 
essentially correct, and in order to obtain a region 
of strong absorption, it is necessary that there be 
a large number of allowable combinations of ¢ in 
the equations of Case I and II for which the 
energy absorbed is practically constant and the 
relation (16) is maintained.*? This implies that 
the energy spectrum in the neighborhood of each 
of the combining o's must be such that its gradient 
is practically zero, which can only occur when 
the o’s are in the vicinity of those regions of the 
zone for which the gradient actually does vanish. 
This implies the need for precise knowledge of the 
zone structure of each lattice dealt with, and an 
examination of the zones described previously 
fer the face-centered lattice shows that there are 
seven values of o for which the gradient vanishes. 
(o’s which are sent into one another by symmetry 
operations are not counted as different). As may 
be expected, these points correspond to the inter- 
action of the zones with symmetry axes, so that 
the polarization direction of the waves may be 
determined at once. Thus in the (100) direction, 
it is evident that there will be one wave polarized 
in this direction (longitudinal wave), and two 
waves of equal energy (differing from that of 
the longitudinal wave) for which the polarization 
directions are orthogonal to this. These points 
are listed below along with the number of differ- 


ent energies going with the three direction of 


8? The importance of this was realized by Born and 
Blackman in their work. 
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polarization. The actual energies in terms of the 
parameters a(1), -, a(10) of the potential 
energy expression will not be given at this time 
since correlations have not yet been obtained. 


1. Intercept of the (111) axes with the cube; there are 
three equal energies. 

2. The origin corresponds to three equal energies. 

3. Intercept of (100) axis with cube. The two transverse 
waves have equal energies differing from that of the longi- 
tudinal wave. 

4. Intercept of the (110) type of axis with the cube. 
There are two different energies corresponding to a wave 
polarized in the (100) direction and two polarized in the 
plane orthogonal to this, the latter two being degenerate. 

5. The center of the hexagonal faces of the truncated 
octahedron. The longitudinal wave is polarized perpen- 
dicular to the face, and has a different energy from that 
of the two degenerate transverse waves. 

6. The center of the edge of a square face of the trun- 
cated octahedron (say the (100) face). There are three 
waves of different energy, one of which is polarized in a 
direction orthogonal to the plane passing through the 
(100) axis and the point. The other two polarization 
directions lie in the (011) and the (011) directions. 

7. The intercept of the (011) axis with the truncated 
octahedron, bisecting one of the edges of a hexagonal face; 
the three waves are of different energies and their polariza- 
tion axes lie along this axis and in the two symmetry planes 
which intersect it, respectively. 


The points corresponding to these cases are 
indicated by the appropriate number in Fig. 9. 

Of these seven cases, it is to be noted that 1 
and 2, 3 and 4, and 6 and 7, respectively, possess 
o's which differ by a vector of the type s‘, so that 
of these three pairs, one member may be assumed 
to have the three, possibly degenerate, energies 
going with the lower part of the discontinuity and 
the other the energies going with the upper part. 
In case 5, for which 20=s"', the discontinuity is 
such that both upper and lower branches have 
values at the same point; that is, the discon- 
tinuity arises from an actual double-valuedness 
of the energy surfaces at this point—a property 
that is intimately tied up with the symmetry 
properties. 

Were the §’s in (15) all nonvanishing, it would 
be safe to say that the spectrum would possess as 
its most prominent secondary lines, those transi- 
tions of the type (23), (23a) and (23b) in which oe 
belongs to the lower side of the discontinuity 
and o; to the upper, or vice versa, and in which the 
signs before both of these are positive. For in 
these cases v(o,) + (02) would vary very slowly as 
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a2 and o; take on values permitted by the relation 
o2+0;=8* since one of these decreases as the 
other increases, and there would be a large accu- 
mulation of transitions of nearly constant energy. 
It is readily seen that there would be 17 possible 
combinations of this type, namely, those listed 
below: 


3 and 4—four lines (two different polarization energies) ; 
5 and 5—ivur lines (two different polarization energies) ; 
6 and 7—nine lines (three different polarization energies). 


Designating these 17 sums, vs(e2)+v,(03), as 
multiplet A, the 17 lines corresponding to (23) 
would appear on the large frequency side of the 
fundamental, v(s), since they are all found to be 
greater than it. Similarly, this same multiplet will 
appear at an equal interval on the large frequency 
side of 2v(s) and 3v(s) corresponding to (23b) and 
(23a), respectively. 

The 17 difference frequencies | v3(o2) —v,(c3) | 
which we shall designate as multiplet B, would 
also be prominent at high temperatures and, 
corresponding to the cases (23), (23b) and (23a) 
we should have the cases 


B, v(s)+B, 2v(s)+B. 

Similarly, if we designate —A by A and —B by 
B, the cases 

2v(s) +B, 2v(s)+A 


—v(s)+A, v(s)+B, 


would be of importance, and indicate that the 
multiplet A would appear between 0 and p(s), 
that the inverted multiplet B would appear be- 
































A 
(23) | _B_ wom 
A inn 
(36)| ~*~ 8B] 8. “ 
A - 
(23a) sites ~ Se 2 <&. 
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Fic. 10a. A schematic representation of the disposition 
of multiplets A, B and their inverses. The dotted lines 
represent lines that will disappear at very low tempera- 
tures. The three rows correspond to results arising from 
the three equations indicated. 
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tween 0 and v(s) and v(s) and 2p(s), respectively, 
and that the inverted multiplet A would appear 
between 0 and v(s) (see Fig. 10a). 

The restrictions on the 6's, however, do not 
allow all of the 617 possible sums and differ- 
ences of (23), (23a) and (23b) to appear, and in 
fact, the cases excluded are those corresponding 
to transitions in which states involving 1, 2, 3 
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course, excludes the 4X3 transitions involving 3 
and 4 mentioned above and also some others not 
yet discussed. This restriction may be brought 
out only by a detailed investigation of the type 
of trinomial expressions in the g’s that may occur, 
the argument being essentially as follows. 

In terms of the normal coordinates ga(c, c) 
and q3(¢, s) the displacement vector for the atom 





of type 1 or 2 situated at /mn is given by 


and 4 are excited in any combination. This of 


x'(Imn)= >> halo, Ofa'(e) palo, lmn) cos o-r(lmn)+ alo, S)fa‘(o) pala, lmn) sinao-r(lmn)}, (24) 


where ga(o) are the final normal coordinates, and o extends throughout the entire cube in o-space, 
fa‘(o), is the ath polarization vector going with the oth wave, which will be identical for like atoms, 
while p.(o, /mn)/p3(o, l'm'n’) gives the relative amplitudes of the atoms at /mn and /'m’n’, assuming 
they have the same phase (this will be 1 if the atoms are identical). For the prominent values of ¢ 
in which we are interested, the polarization directions are identical for unlike atoms, so that we need 
not carry the indices /mn in these quantities. 

Using (24), the quantity £'(/;m jn, lmn) is 


YY tale, c)fa'(a) (plo, limin;) cos o-r(lymjn;) — po, lmn) cos o-r(/mn)) 


+ alo, S)fa'(o) (plo, limn;) sin o-r(limym;) — p(lmn) sin o-r(lmn))} (25) 
and if we set r(/jmjn;) =r(/mn)+p(limjn;), this becomes 
=P I [galo, c)fai(a)A(o, limin;) —galo, 8)fa‘(e) Blo, limn;) | cos ¢-r(lmn) 
+[qalo, S)fa'(a)A(o, limin,) +qalo, 5)fa‘(o)Blo;l.min;) |] sin o-r(ljmjn;)| (26) 
where 
A (a, limin;)= —(p(o, lmn) — plo, limjn;) cos a: p(limn;)) 
(27) 
B(o, limn;)= — po, limn;) sin o- p(ljymjn;) 


so that A is invariant under an inversion in /mn, while B goes into its negative. If terms of this type 
are substituted in one of the invariant forms of §2 and the result is summed over all values of /mn 
corresponding to like atoms, the result may be written in the form 


> ymin) & f{a(os, limn;)a(oe, limin;)a(o3, limjn;) —a(o1, limn;)b(o2, Lymjn;)b(o3, limin;) 
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—b(o;, limn;)a(oe, limjn;)b(a3, lym n;) — b(04, Lymn;)b(o2, limjn;)a(o3, Limjn;)\ (28) 


where o}, o2, 7; must satisfy the relation (16), /;m jn; are summed over the values occurring in a given 
form with the proper coefficients y, and 


a(o, limin;=)> { dale, ¢)fa'A(o, lLmini) —qalo, 5) fa'B(o, limn;)}, 


(29) 
b(o, lamn)=>d {galo, s)faiA(o, limn:)+ gale, c)fa'B(o, limn;)}. 


a 


Now since £'&£*(/+Aiym+uyu, n+v) goes into — £'¢/E*(1—A, m—p, n—v) under inversion in the point 
lmn, the expression 


EEK +N, m+y, n+v)— Ee E(I—d, m—p, N—v) (30) 
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INFRARED ABSORPTION 


must occur in each cubic invariant form if the 
first term does. As a result of this, it follows at 
once that in the result (28), when (29) is substi- 
tuted, all terms not containing B(o,/;mjn,) an 
odd number of times will vanish. Closer inspec- 
tion shows quite readily that this fact does not 
depend upon the assumption that all three prod- 
ucts in a given trinomial possess the same /mn 
and /m jn;. 

Now in the cases 1, 2, 3 and 4 described 
previously, sin o-p(/;mjn;)=0 for all values of p 
so that all intercombinations of these can give 
rise to no secondary maxima since the necessary 
combination of g’s do not appear in (15). 

For the combination of 5 and 1, and 6, 7, and 1, 
in the manner described above, this does not oc- 
cur and all possible transitions are allowed, so 
that there are 13 frequencies v3(¢2) +¥v,(¢3) arising 
from these two cases, and the multiplets A and B 
as well as their inverses will contain this number 
of terms. 

From the restriction discussed in the previous 
paragraph, it follows that the frequencies 2(s), 
3v(s) and 4y(s) may not occur as a combination 
of cases 1 and 2, as far as the cubic terms are con- 
cerned; moreover, it is easily seen that the case 
in which o;+s cannot lead to anything new if all 
restrictions are taken account of properly. 

There is one additional possible source of 
secondary maxima, namely, those cases in which 
o2= —o;3 (and o,=s), so that only transitions in 
which waves on the lower side or on the upper 
side play a role along with s are concerned. Such 
transitions undoubtedly have probabilities com- 
parable with those considered above, but will 
probably not lead to maxima of the same in- 
tensity since the v(e2)+v(¢3) will not be constant 
for as large a range of o; as in the preceding cases. 
It follows as before that only cases 5, 6 and 7 
may occur along with 1. 

At low temperatures, which means tempera- 
tures well below the characteristic temperature 
of the crystal, the peaks corresponding to cases 
in which negative signs occur in (23), (23a), 
(23b) (dotted lines in Figs. 10a and b) will dis- 
appear and only those in which all are positive 
will remain. The number of lines corresponding 
to these is 93 as compared with 260 in the high 
temperature case. 
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ParT C. ON THE CORRELATION BETWEEN 
EXPERIMENTAL AND THEORETICAL RESULTS 


If an attempt is made to find a satisfactory 
correlation between the preceding experimental 
and theoretical results, it is readily seen that the 
former is not complete enough either to confirm 
or deny the correctness of the latter. In the first 
place, it is essential that the fundamental fre- 
quency be unmistakably identified, and this is 
seen to be impossible at present because of the 
limitations on the long wave-length side of the 
spectrometer; that is, the absorption peak at 
15.34 would be the fundamental if it were 
definitely known that there were no additional 
high absorption regions beyond this, but that 
cannot be definitely stated. If the work of Strong 
can be accepted as correct, it would indicate that 
another region of absorption exists in the vicinity 
of 23u. This difficulty may be removed by using 
another prism, such as one of KCl, or by using a 
grating spectrometer. 

Second, the resolving power in the short wave- 
length region (5 to 10u) of the present instrument 
is not sufficiently great to guarantee that the 
secondary structure has been determined to the 
same degree of accuracy as that in the region from 
10 to 154; and, as a result, the position and num- 
ber of peaks in the former region is not sufficient 
to test the symmetry about v and 2», assuming 
the 15.3 line is the fundamental. 

Third, the broadening of all of the lines, be- 
cause of the influence of high temperature on the 
coupling between waves, makes the precise posi- 
tion of the maxima very difficult to determine. 
This difficulty could be considerably reduced by 
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Fic. 10b. The plan is the same as that of Fig. 10a, the 
difference being that the multiplets arising from transitions 
involving waves on but one side of a discontinuity are 
shown. 
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cooling the observed sample to sufficiently low 
temperatures. 

Concerning the precise position of the funda- 
mental the following theoretical remarks may be 
made. If we make the very rough assumption 
that a(1) and a(2) are predominantly larger than 
any of the other 10 coefficients in B §2, and that 
these are related to the components of the stress- 
strain tensor, the fundamental frequency is found 


to be 
v=(1/v2r)(A/M)?, 
where 
M=myme/(m,+me), A=aleutcs), 


a is the nearest distance between unlike atoms, 
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and ¢,; and cy; are the usual elastic coefficients. 
If the corresponding wave-lengths are deter- 
mined using the results for the elastic coefficients 
discussed above, it is found that for NaCl and 
MgO the fundamental is at 49.54 and 17.5u 
respectively. Since the former is actually at 61.1, 
this suggests that the theoretical results are 
definitely low and indicates that the true value 
for MgO may be somewhat farther in the infrared. 

In order to obtain experimental data that are 
capable of being compared with the preceding 
theoretical considerations, we are undertaking a 
thorough study of the absorption of the MgO 
crystal at low temperatures by means of an 
echelette grating spectrometer. 
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The Spectra of Strontium Hydride 


W. R. FREDERICKSON AND M. E. HoGan, JR., Syracuse University AND WILLIAM W. Watson, Yale University 
(Received July 25, 1935) 


New SrH band systems photographed at high dispersion 
include: (1) an abbreviated C*XZ—N?2Z system with (0,0) 
band head at 3808A, (2) a many-lined D?2—N 2Z system 
between 5000A and 6700A, (3) an E 2*II—N 2 system con- 
sisting of a Av=0 sequence only with the (0,0) band at 
5323A. The (0,0) and (0,1) C bands have three violent 
perturbations and a sharp cut-off due to predissociation 
at A’=19. The form and magnitude of the perturbations 


INTRODUCTION 


T has previously been observed! that strontium 
hydride exhibits two band systems in the 
red region of the spectrum, a ?2—*L system with 
principal head at 7020A and a *II-—> system 
whose (0,0) band has heads at 7347A and 7508A. 
These bands lie slightly towards longer wave- 
lengths from the analogous systems of the CaH 
molecule. Since the lowest electronic levels of Sr 
and Ca are so very similar one would expect a 
fairly complete correspondence between the 
hydride spectra of these two elements. 
In addition to the two CaH systems in the red 
there are three other systems? in the visible and 


1W. W. Watson and W. R. Fredrickson, Phys. Rev. 39, 
765 (1932). 

*>Cf. W. Jevons, Report on Band Spectra of Diatomic 
Molecules for references. 


are computed. Quantum analyses give By=3.930 and 
Do= —9.31 10-4 for the C state, B, =3.56 for the normal 
state. Quantum assignments of the lines of the six principal 
branches of the (0,0) band of the E system are presented. 
For the E state A =117, B*o, _5=3.6388, B*o, .; =3.8687 
and perturbations at small J values distort the A-doubling 
curve for the *II; state. 


near ultraviolet. An interesting *"2—?S system 
consisting normally of but two bands, (0,0) at 
3534A and (0,1) at 3696A, exhibits a sharp cut- 
off at K=9 of the v’=0 level due to predissocia- 
tion. There is also in the range 4000A—5800A a 
weak *>—*Y system whose bands have no 
marked heads and whose upper vibrational levels 
all display violent perturbations. Finally, the 
existence of a *II—*r band at 4900A with Q 
branches having very closely spaced rotational 
lines has been noted. SrH band systems corre- 
sponding to all of these have been located by 
the writers. In each instance the system is 
displaced slightly towards the red from the 
position of the similar CaH system. We shall 
discuss in this paper the abbreviated *2—*r 
bands in the near ultraviolet and the *II—-?> 
system in the green region only. 
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Fic. 1. Microphotometer trace of the (0, 0) C band of SrH at 3808A. Note the multiple 
perturbations and the sharp cut-off due to predissociation at K’=19. 


For the sake of brevity and to emphasize the 
similarity between the spectra of these two 
molecules, we shall adopt the same alphabetical 
designation of states as used for CaH. The 
normal state, common to all the band systems, 
is N?=*+. The upper states of the red bands are 
A *II and BS, that of the near ultraviolet system 
is C?, for the many-lined spectrum from 5000A 
to 6700A the upper state is D?Z, and for the 
green bands F II. 


TABLE I. Assignment of frequencies in the C*X—>N?*Z(0, 0) 
band of SrH (cm™ units). 














Int. v (obs) R(K"’)| Int. v (obs) P(K"’) 
3 26548.88 18 1 6291.72 1 

3 47.04 18 2 84.88 2 

1 39.49 17 3 78.68 3 

1 38.03 17 1 77.36 19 

4 11.76 16 2 75.31 19, 13 
4 10.68 16 4 73.07 4, 20 
4 497.16 15 1 71.85 20 

4 95.22 15 0 70.09 14 

2 76.32 14 3 68.30 5 

2 74.40 14 0 65.33 8 

1 72.20 12 6 63.38 6, 18 
1 63.69 12 6 62.14 17, 18 
5 61.79 13 6 59.96 8,17 
5 60.85 13 N 58.63 7 

5 46.81 12 4 56.96 9, 10 
6 45.58 12 5 55.04) 8, 10 
6 31.88 il 53.56; 15, 16 
0 25.70 10 7 53.56 | 11, 12. 12 
4 18.19 10 51.40} 13, 13, 14, 16 
4 17.09 10 0 35.44 9 

1 08.32 9 

1 06.23 9 

4 03.85 9 

4 392.18 8 

4 91.85 x 

3 81.44 7 

3 79.78 7 

0 74.94 6 

1 68.32 6 

1 62.42 6 

0 58.19 7 

3 53.22 5 

4 43.12 4 

3 33.21 3 

3 23.91 2 

2 15.15 1 

1 06.94 0 








Preliminary spectrograms were made with a 
small 120-cm grating, while the spectrograms 
from which measurements were taken were 
obtained in the second order of the 21-foot 
grating in a stigmatic mounting, the average 
dispersion being 2.4A per mm. The light source 
was a d.c. arc carrying 2 amperes between a 
water-cooled copper cathode and a hollow copper 
anode filled with metallic strontium in an 
atmosphere of hydrogen at about 7 cm of Hg 
pressure. Eastman 40 plates were used for the 
C system, Eastman I-C plates for the E bands. 
The exposure time for each group was two hours. 


THE C SYSTEM 


The C system of SrH consists merely of a (0,0) 
band with its head at 3808A, exhibiting multiple 
perturbations and a sharp cut-off at K’=19, 
together with a considerably weaker (0,1) band 
at 3984A, both bands degrading towards the 
violet. The sudden cessation of lines in these 
bands is undoubtedly caused by predissociation 
into the continuum of the less stable D*> state. 
In Table I are listed the intensities, wave 
numbers and quantum assignments of the lines 
of the (0,0) band. These assignments are made 
with the aid of lower state combination differ- 
ences known from the earlier work on the red 
band systems. The nature of the perturbations 
may be seen by inspection of Fig. 1 which is a 
microphotometer trace of this band. The complex 
structure of this band is in decided contrast 
with the simplicity of the corresponding CaH C 
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Fic. 2. Graph showing departures of observed R branch 
lines of the (0, 0) C band from positions predicted by fourth 
degree least-squares equation for the unperturbed lines. + 
indicates observed line has a larger wave number than the 
predicted value. 


band. Whereas the lines in the CaH band are 
single, most of those in the SrH band are close 
doublets due to the larger coefficient y= 0.122 
of the spin doubling in the normal state. The 
first perturbation appears at K’=7 although the 
lines belonging to K’= 6 are broad and somewhat 
displaced. Therefore in evaluating the constants 
of the upper state it is possible to use only the 
first five of the appropriate combination differ- 
ences. This somewhat lessens the accuracy of 
this calculation. From these five differences we 
compute Bo’ = 3.930 and Dy’ = —9.31 10-4. 

If one substitutes the By constants in the usual 
expression for R-branch frequencies (neglecting 
the terms involving the D constants since K is 
small), the wave numbers of these first five lines 
may be calculated. In order to get agreement 
with the observed frequencies, it is evident that 
an extra linear term —0.21K’ must be included. 
This term is considerably smaller than the 
corresponding extra term —2.22K’ in the CaH C 
band.* 

The nature and magnitude of the perturba- 
tions have been determined in the following 
manner. We observe that the R branch lines for 
K”=9 and 11 are intense and single while 
R(10) has the splitting required by the spin 


3B. Grundstrém, Zeits. f. Physik 69, 235 (1931). 
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doubling in the K’’=10 level. Assuming, then, 
that these three lines (taking the mean of the 
two R(10) lines) are practically in their normal 
position, and by using them in conjunction with 
the first six lines of the R branch, the following 
fourth degree equation representing these fre- 
quencies has been determined by the method of 
least squares: 


v= 26300.00+ 6.719.117 +0.3958 \/? 
— 0.00329.M*+ 7.32 X 10-8 M4. 


The variation of the position of each of the other 
R branch lines from that calculated with this 
equation is shown in Fig. 2. These variations 
which yield curves of the ‘‘anomalous dispersion”’ 
type indicate the presence of three perturbations. 

In the construction of Fig. 2 the assignment 
of the perturbed lines to either the R; or R2 
branch has been made so as to give the best 
agreement with the known AoF,"’(J) and Ae F2”’(J) 
values, using P branch lines of about the same 
intensity. This assignment accounts for every 
line on the low wave number side of the origin, 
including even a weak line lying some 15 cm™! 
beyond the head of the band. In the region of 
R branch lines but one line at 26,413.82 remains 
unassigned. The characteristic assignment of the 
weaker, greatly perturbed lines to transitions 
already represented by stronger, less perturbed 
lines is to be noted. 

Since the corresponding band of CaD has so 
many perturbations‘ in contrast to the simplicity 
of the CaH band, the study of the hydrogen 


TABLE II. Assignment of frequencies in the C?Z—>N?*Z(0, 1) 
band of SrH (cm™ units). 











x” R P —" R a 
3 25165.11 25110.20 11 269.86 90.00 
4 75.10 04.93 12 285.27 91.32 
86.48 92.59 
5 Masked 01.13 
13 302.82 94.57 
6 096.42 04.02 
7 211.50 Masked 14 319.51 97.31 
13.20 21.29 98.65 
8 225.62 087.56 15 25341.95 25099.95 
24.48 92.59 43.86 101.13 
9 238.07 91.32 16 59.70 02.19 
92.59 61.19 
10 253.13 90.26 17 89.88 12.07 
54.05 91.67 14.26 











4W. W. Watson, Phys. Rev. 47, 27 (1935). 
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TABLE III. Assignment of frequencies in the E *11—~N?*Z(0, 0) band of SrH (cm units). 














JI” +4 Pi Q1 Ri Pe Q2z Re 
1 18972.59 
2 18868.95 18883.22 18958.56 77.53 
3 18850.98 72.11 95.04 18936.64 56.35 82.71 
4 47.20 75.17 905.04 26.98 53.98 89.30 
5 43.74 78.30 15.77 17.17 52.80 96.25 
6 39.49 81:47 26.86 09.48 52.00 19002.97 
7 36.11 84.81 38.23 02.45 52.00 10.28 
& 32.57 88.35 49.10 894.27 52.50 
9 29.58 91.49 59.81 87.50 53.12 26.38 
10 26.16 95.04 69.18 81.47 53.98 34.98 
11 22.41 98.99 81.56 75.15 55.38 43.79 
12 17.56 902.45 91.99 69.56 56.94 53.00 
13 15.84 05.88 19002.97 64.42 58.86 62.48 
14 12.76 09.48 13.84 59.56 61.28 72.11 
15 09.23 13.22 24.88 55.12 63.69 82.09 
16 06.23 17.17 36.38 50.98 66.38 92.23 
17 03.43 20.50 46.96 47.20 69.18 102.48 
18 00.10 25.41 58.20 43.74 72.59 12.88 
19 798.27 29.43 69.48 40.56 76.27 23.62 
20 96.20 33.71 80.78 37.54 79.55 

21 93.86 38.23 92.23 34.94 83.58 45.22 
22 91.96 42.62 103.60 32.57 87.67 56.24 
23 90.09 47.07 15.07 30.41 91.99 67.62 
24 88.67 52.00 26.45 28.60 96.25 78.39 
25 87.11 56.77 38.03 27.20 19000.76 89.70 
26 85.91 61.57 49.42 26.16 05.22 200.96 
27 84.80 67.15 60.97 25.35 09.43 12.27 
28 83.81 71.68 72.62 24.23 14.84 23.44 
29 83.13 84.05 23.50 34.72 
30 82.58 95.16 45.80 
31 82.00 206.70 57.07 
32 81.40 17.96 68.05 

$1.00 29.02 

40.06 

51.17 

61.58 

73.77 

84.45 





isotope effect in this SrH band should prove 
interesting. This investigation will be reported 
in a forthcoming paper. 

In Table II are listed the frequencies and 
quantum assignments of the (0,1) C band. The 
analysis of this band is somewhat hindered by 
the low intensity of its lines and by the presence 
of several strong atomic lines. Perturbations 
identical with those of the (0,0) band are of 
course present. The few final state combination 
differences yield the values B,’’=3.56 and 
D,"” = —1.7X10-. 


THe E SYSTEM 


The Av=0 sequence of the E system of SrH 
has a principal head at 5323A, the bands de- 
grading to the violet. Because of the near 
equality of the rotational B constants in the two 
electronic states involved, other sequences are 
apparently of negligible intensity. Analysis shows 
the bands to be of the *II-—*S type. Since the 
lines in the P and R branches are very dastinct 
at a distance from the origin, quantum assign- 
ments are easily made. for these lines with the 
help of the combination differences for the final 


N?*> state. These assignments show the £?Il 
state to be regular. 

Below J” =8 the assignment of lines becomes 
very uncertain. This uncertainty is caused by 
the presence of the numerous lines of the (1,1) 
and (2,2) bands as well as lines of the D system 
which extend into this region. Also there are 
undoubtedly perturbations at the low J values. 
It is quite evident that in the *II, sub-state with 
v=0 perturbations occur at J’=16} and 10}. 
In Table III are listed the assignments of 
frequencies for the Pi, Ps, Qi, Qe, Ri and R: 
branches of the (0,0) band. 

The assignment of quantum numbers to the 
QQ; and Q: branches is based on A-doubling 
considerations. Since A/B=31 and since the 
nature of the doubling indicates that the p 
constant must be positive, it is assumed that the 
A-doubling of the FE ?II state should lie somewhat 
between that exhibited by SiH and ZnH.® Fig. 3 
shows the A-doubling for the assignments chosen, 
Due to the uncertainty of the P;- and Ry-branch 
assignments for J values below /=8, the A- 
doubling curve in this range for the Il, levels 








*R.S. Mulliken and A. Christy, Phys. Rev. 38, 105 
(1931). 
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Fic. 3. A-doubling curves for the /°II state of SrH. 


may not be accurate. If it is correct as given, 
the failure of the curve to return to zero at J=0 
must be attributed to another perturbation. 
Because of this possibility of error in the P; and 
R, branches at low J values, however, we do 
not evaluate the usual fo and go quantities 
involved in the A-doubling. The other molecular 
constants, determined from the A »F values in the 
customary manner,’ are presented in Table IV. 

Satellite branches should be observed in this 
band. We have located the “Q2; branch at large 
J values but the other satellite branches cannot 
be definitely established because of the many 
weaker lines of the other bands present in this 
region. 


DIEKE 


TABLE IV. Rotational constants of the E*Ml—-~N2?Z(0, 0) 
SrHI band at 5323A. 








N 2Y state: Bo =3.6198 Do = —1.287 X 10-4 
El state: Bo,* —} =3.6388 Do = —1.03 X 10-4 
Bo,* +4 =3.8687 Do = —1.95 X 1074 


A =117 


ELECTRONIC STATES OF SrH 


In our earlier work! on the red bands of SrH 
we noted that the A *II and B °S states stand to 
each other in the relation of ‘‘pure precession,” 
being undoubtedly 5po*5pzx7Il and 5po*5po 2X 
respectively. The NV **+ normal state is probably 
5so0*4do *X. The E state may be either 5so*4dz ?II 
or 5so5po5pr 711. The D state might then well be 
S5sa5po"*X, and the lack of strong interaction 
with the E state could be due to the large 
difference in the B constant of the D state from 
that of all other states. For our analysis of one 
of the bands of the D system, with v’ as yet 
unknown, indicates B’ to be only 1.91. The lack 
of strong interaction may also be due to the 
presence of the C*S state above the E state, 
thus neutralizing in part the effect of the D state. 

The D system of SrH as well as the C and B 
systems of SrD will be discussed in a further 


paper. 
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The 3p*> —2s*> Bands of HD and D, 


G. H. Dieke, The Johns Hopkins University 
(Received July 18, 1935) 


The 3p*2—>2s° bands of HD and D, which form the most conspicious regularities in the 
near infrared are given and the band constants calculated. 


N a previous communication! the so-called 
Fulcher bands (3p*II-—+2s*2) of the isotopes 
HD and D; of the hydrogen molecule were given. 
There is a band system 3p*2—>2s*> chiefly lying 
in the extreme red and near infrared closely con- 
nected with the Fulcher bands. It has the same 
final state, and the initial states belong to the 
same complex, i.e., they differ only by the value 
of the projection A of the orbital momentum on 
the internuclear axis. 
~ tG.H. Dieke and R. W. Blue, Phys. Rev. 47, 261 (1935). 


The present paper contains chiefly the empir- 
ical data of the 3p°2—>2s*= system. A theoretical 
discussion of the 3p complex is given in the fol- 
lowing paper, which will also contain an account 
of the perturbations in both systems. The general 
remarks about the production of the spectrum, 
the wave-length measurements, etc., made in the 
previous paper! apply also to the present in- 
vestigation. 

In the far red and infrared the sensitivity of 
the photographic plates varies rapidly with the 








BANDS OF 


HD 


AND D, 


TABLE I. The 3p'=—>2s*d bands of HD. 

















P BRANCH R BRANCH P BRANCH R BRANCH P BRANCH R BRANCH 
v I v J v I v I P v I v I 
_ 4 666.08 9 985.87 5 4-3 
0 _ — 11 664.67 2 5 574.12 9 961.41 4 0 _ _ 11950.77 4 
1 11574.34 1 694.49 4 6 472.42 7 925.10 3 1 1187462 3 970.71 5 
2 514.36 3 713.90 4 7 361.56 5 2 819.05 1 978.51 5 
3 444.97 3 722.78 4 & 242.21 3 3 752.12 2 974.16 4 
4 366.59 2 721.22 2 9 114.83 1 q 674.32 1 957.75 3 
5 279.73 1 709.06 O 5 586.21 0 929.41 2 
3-1 6 889.33 0 
1—0 0 - — 14 559.84 7 
0 — — 13 465.77 2 1 14477.34 0O 580.96 8 5-2 
1 13377.31 4 491.64 10 2 415.97 5 590.42 6b 0 — —_ 15 330.13 4 
2 315.37 5c 505.20 10 3 342.07 6 569.03 5 1 15254.16 4 343.93 5 
3 242.09 5 506.34 10 4 257.98 4 551.11 6 2 192.39 4 341.43 5 
4+ 157.88 6c 495.14 10 5 144.89 5 §16.13 5r 3 115.32 4 322.68 4 
5 063.27 5 471.79 0 6 037.72 4 4 023.30 3 287.72 4 
6 12958.84 5 436.34 1 7 13916.10 2c 5 coincides with Dg 
7 845.16 4 388.39 0 6 169.93 2 
& 722.88 2c 3-42 
9 592.63 3 +De 0 — —_ 12 459.88 7 5-3 
1 12379.33 7 483.05 7 0 — — 13 330.16 2 
1—1 2 322.13 8 496.62 9 +D2/1 13256.17 3 345.97 2 
0 11 262.98 0 3 254.43 9 481.43 6 2 198.48 3 347.51 3 
1 290.98 1 4 178.46 6 471.69 6 3 127.41 4 334.78 2 
2 308.77 0 5 075.59 5 446.90 4 4 043.38 4 307.80 2 
3 316.41 0 6 980.54 4 5 12946.87 5 266.77 1 
4 313.49 00 7 6 838.43 3 
5 
4—1 6-3 
2—0 0 — _ 16 050.74 2 0 — -—— 14 586.93 4 
0 _ —_ 15 165.34 6 1 15970.62 2 066.61 3 1 14615.05 3 598.37 5c 
1 15078.98 6 187.32 6 2 906.81 2 066.24 3 2 455.29 4 593.49 6 
2 015.02 7 195.00 & 3 827.66 2 049.70 2 3 379.80 4 571.55 7 
3 1493779 6 188.37 6 4 733.67 1 017.18 3bv 4 289.31 4 534.70 5 
4 847.70 7 167.47 ir 5 625.53 0O 5 183.63 2 
5 745.23 7 132.56 6 6 503.70 00 6 065.32 1 
6 631.14 5 083.85 5 
7 505.88 5 021.64 4 4-2 6 +4 
8 370.28 3 14946.35 3 0 oo = 13950.71 5 0 — _ 12684.92 4 
225.81 1 1 13872.65 6 968.64 9 1 698.35 3 
2 812.95 7 972.41 9 } 2 12559.18 4 +De 697.34 3 
2-1 3 740.00 6 962.05 9 an,Qs | 3 489.62 5 681.37 3 
0 _ — 12962.60 7 4 654.22 7 937.64 4 4 406.90 5 652.33 3d 
1 12878.32 6 986.71 8 5 556.18 4 899.34 4 +D:/5 310.38 3 
2 818.64 8 998.65 5 +D:/ 6 446.43 2 847.43 1 6 
3 747.76 9 998.15 4r 7 325.66 1 781.62 2 
8 194.85 1 +Ds 





wave-length. The photographic intensities de- 
pend very much on the plate used and the sensiti- 
zation employed. Relative intensities appear 
sometimes completely distorted even over a com- 
paratively narrow region. This accounts for al- 
most all the apparent intensity anomalies found 
in the tables. 

The frequencies and intensities of the bands 
are given in the Tables I for HD and II for Ds. 
The bands contain a large part of the strongest 
lines in the infrared and are the most prominent 
feature in that part of the spectrum. Only those 
bands which are fairly complete are listed. There 
are several weaker bands, the position of which 
can be calculated and of which a number of lines 
can be found. They, as well as bands with higher 
values of V’, are left for a later occasion as un- 
doubtedly their identification will be much more 
certain when the other strong lines in this region 
have been classified. 

The rotational constants were calculated in the 


usual way,! and they are listed in Table III. The 
spaces left open are the vibrational levels which 
are so strongly perturbed that no reliable values 
can be obtained. 

As shown more in detail in the following paper, 
the B, values calculated directly from the obser- 
vations cannot be used for the calculation of the 
moment of inertia or other constants which have 
to do with the vibration and rotationless state of 
the molecule. It is necessary to take into account 
the influence of the interaction with the 3p*II 
state. If we correct for this in the way indicated 
in the following paper the new values for B, are 


HD D, 
Bo 20.257 13.553 
B, 19.259 13.045 
B; 18.236 12.504. 


This includes the B, of all those levels which are 
far enough from any 3)II level, so that only the 
regular A doubling has to be taken into account. 

B,. can then be calculated from B,=8B, 





608 G. H. DIEKE 


TABLE II. The 3p*X—2s*S bands of D>. 





P BRANCH R BRANCH P BRANCH R BRANCH P BRANCH R BRANCH 
J v I v } J v v J ¥ rf v 
0-0 4+-+1 5 465.60 4 685.42 4 
0 _ - 11676.15 0 0 15 402.94 2 6 395.27 6 652.89 2 
1 696.39 0 | 1 15347.75 1 415.83 2 7 317.431 
2 11575.39 1 709.82 3 | 2 305.68 3 418.84 4 8 233.36 4 
3 528.90 0 716.39 1 Ss 254.16 2 411.95 2 9 140.57 0 
4 476.25 1 716.14 2 4 193.39 3 395.20 4 7-33 
5 417.65 00 709.29 1 . 123.61 1 368.66 1 0 _ 15 385.52 2 
6 353.38 0 695.27 0 6 045.09 0 332.44 2 | 1 15 336.14 1 393.59 1 
7 958.12 0 286.68 0 | 3 9501 2 300.78 3 
1—0 8 862.88 1 231.56 Od 3 243.15 2 375.96 2 
0 - —_— 13 163.44 3 9g 759.97 0 | 4 181 02 > 350.49 3 
1 13 103.88 2 181.55 x 10 649.58 00 | § 107.72 0 313.84 lv 
2 062.67 4 191.80 4 P 024.820 O 66.09 1 
| » at. £00. 
3 014.06 3c 194.21 4d 4-2 : 93181 0 07.16 Ode 
4 12958.25 6 188.66 4 0 . 13657.08 4 _ - ia 
5 895.44 5 175.35 1 1 13603.08 3 671.16 4 ea 
6 825.92 x 154.26 2 2 563.27 9 676.43 6 | O . = 13 773.41 2 
7 749.94 3 125.61 0 3 515.18 10 672.95 3 1 13725.11 1 782.60 1 
8 667.80 3 089.43 4 458.94 § 660.71 6 2 686.18 4 781.96 2 
9 579.79 5 5 304.74 2 639.80 3 | 3 637.59 3 770.44 2 
10 486.28 2 6 322.80 6 610.28 | 4 579.80 3 749.27 2 
7 243.74 2 572.32 ¢ } 5 51193 5 717.84 3 
2-0 s 157.50 1d 526.04 1 |} 6 435.44 01 670.83 3 
0 - . 14.583.13 3c 9 064.32 0 7-5 
1 14524.61 3 599.13 4 10 964.79 1 | Oo 7 12 226.45 2 
2 482.36 6 606.22 5 | 1 12179.25 1 236.70 1 
3 431.64 4 604.40 4 4-3 2 238.22 2 142.43 3 
4 372.66 4 593.66 Sr 0 11978.87 3 3 097.06 1 229.91 1 
5 305.66 3 574.07 2 1 11925.94 2 994.00 3 4 043.55 3 213.05 2 
6 230.89 4 545.71 4 2 888.39 5 12001.52 Se | 5 1198099 2 187.11 1 
7 148.93 (10) +HD 508.76 1 3 843.56 1 001.35 3c | 6 910.89 2 152.15 1 
8 059.22. 1 463.43 (4) ass,Q: | 4 791.73 2 11993.54 4 a3 
9 5 733.11 O 978.14 1 . <a 
mee = 0 - _ 16 354.17 00 
10 $ 85 0 aa = 1 16306.03 00 359.57 00 
seit . = _ a = 2 263.600 352.46 0 
@ _ = 12:768.25 ¢ ‘ bobs 3 209.25 1D 332.93 00 
: "sma & . | 4 142.76 2 301.22 00 
1 12710.91 4 785.42 5 5-2 | = 06467 00 
2 670.95 7 794.83 8 0 . - 1469.00 4 |” sso7ssa 0 
3 623.74 7 796.49 6 1 14815.99 3 $81.02 5 aoi, Ro | sila 
4 569.38 9 790.37 7 2 775.22 5 883.28 6 | 8-4 , 
5 508.12 7 776.54 3 3 725.04 4 876.25 5 0 - — 14 742.04 3 
6 440.23 6 755.09 4. 4 665.83 § 861.98 5 | 1 14695.06 748.53 3 
7 365.94 4 726.14 3 5 598.10 3 $23.48 1 2 654.81 3 743.63 5 
8 285.57 4 + 689.73 2 6 524.23 1 3 603.54 3 727.35 4 
9 199.432 646.08 1d 7 427.32 3 +HD 4 541.44 4 700.00 5 
10 107.74 2 5 468.86 2de 
S83 6 386.18 1 
3—0 0 som : 13 190.79 2 85 
0 — -_ 15934.99 0 1 13138.91 0 203.86 4r 0 _ “ai 13195.18 id +HD 
1 15877.57 0 948.84 0 2 100.29 4 208.44 4c 1 13149.22 0 202.66 0 
2 $34.23 3 952.75 1 3 053.46 3 204.71 3 2 111.09 1 199.97 1 
3 781.35 1 946.68 O 4 1299865 5 +HD 194.85 1 3 063.27 (5) +HD 186.84 0 
4 719.16 2 930.64 3 5 936.46 3 161.74 0 4 005.24 2 1603.78 0 
5 647.91 0 904.74 1 6 869.16 3 5 12937.96 1 
6 567.86 1 869.05 1 7 779.89 1 6 861.63 2 
7 479.29 00 9-4 
8 382.55 00 6-2 0 as _— 15607.40 1 
0 — , 16.006.34 0 1 15562.27 0 610.26 1 
3-1 1 15954.63 0 015.70 O ) 520.17 1 5090.88 2 
0 — — 14.120.06 4 2 912.52 4 014.00 3 3 46527 1 376.22 1 
1 14063.78 3 135.10 5 3 859.70 2 001.26 0 4 30772 1 530.19 2 
2 022.79 10 141.33 8 4 796.48 3 15977.62 (4) ax, Ps} = 31769 0 488.72 0 
3 13973.42 5 138.73 5 5 723.15 0 942.88 00 6 25.36 (0 
4 915.86 10 127.31 5 6 639.74 1 897.36 2 ei. 9-5 
5 850.34 5 107.16 3 a 
6 777.16 3 078.35 3 6-3 0 ~~ ~ 14060.33 1 
7 696.562 040.90 34 +HD| 0 - — 14328.03 4 4 | | 1401633 0 —s 
8 608.85 3 13995.33 1 1 14277.47 2 338.60 (10) M0: | 2 1397639 2 056.09 3 
9 514.39 2c 2 237.50 10 +HD 33005 4 THD) 3 924.65 4d 035.64 1 
10 413.50 Or 3 188.08 4 329.65 4 4 861.43 3 002.98 3d 
4 129.28 5 310.34 5 5 786.74 0 957.84 0 
3-2 5 061.36 3 281.20 2 6 
0 _— -- 12 374.23 6 6 984.63 4 242.23 2 9-6 
1 12318.99 8 390.40 5 7 899.34 4 +HD = 193.58 1 0 — — 12577.78 5 
2 280.39 7 398.91 8 1 1253467 0O 582.55 46-4,P3 
3 234.45 6 399.71 6 6-4 2 496.87 3 576.58 3 
4 181.41 7 392.83 7 0 _— — 12715.98 3 3 448.35 1 559.18 4 
5 121.54 5 378.27. 5 1 1266644 3 727.54 3 4 389.29 4 530.79 2 
6 055.06 7 356.20 5 2 628.74 7b 730.27. 5 5 319.86 1 490.82 Ob 
7 11982.25 3 326.65 2d 3 582.55 4 9-6,R: 724.14 3 6 
8 903.36 3 289.75 3 4 528.10 7 709.18 5 7 
8 
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TABLE III. Rotational constants of the 3p*D level. 











HD D2 
By D By D 

0 20.077 0.0089 13.475 

1 19.100 87 12.974 0.0042 
2 18.095 77 12.440 39 
3 + 11.910 39 
4 16.044 71 11.367 37 
5 15.011 77 + 

6 + 10.201 34 
7 + 

& 8.920 27 
1) 8.024 39 





! 
| 
| 





—a(l’+1,2). In evaluating these constants the 
By value of De which is less reliable than the 
others was disregarded. 

We find then 


HD D. 
B, 20.766 13.856 
a 1.010 0.541 
I, 1.3319-10-% 1.9961 -10~* 


le 1.09662 -1078 1.09649 - 1075 
average ro = 1.09655- 1078. 


The ratio of the B, values is 0.6672 which com- 
pares well with the theoretical value 0.6671. The 
ratio 0.534 of the a is lower than the theoretical 
value 0.544, but the accuracy with which the a 
can be determined in this case is hardly better 
than this agreement. 

If the constants would be calculated in the 
same way from the observed B, without the 
correction, the ratio of the B, values would be 
().6692 which is considerably different from the 
theoretical value 0.6671. 

With the rotational energy of initial and final 
states known, the origins of all the bands can be 
calculated and they are given in the Tables IV 
and V. The values marked with a + are the 
origins of those bands which have strong pertur- 
bations. They are obtained only from the P(1) 
line as the J’=0 level cannot be perturbed.’ As 
they are calculated from one line only they are 
less accurate than the other values. 

The constants in the formula for the vibra- 
tional energy 


wt — xv? + yo? —svt+ +++ (v= V+1/2) 
2G. H. Dieke, Phys. Rev. 47, 870 (1935). 
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TABLE IV. Zero lines of the 3p*X—2s*d bands of HD. 








yr 

Vv’ 0 1 2 3 

0 11 624.55 

1 13 427.56 11 224.80 

2 15 129.24 12 926.46 

3 14 525.41* 12 525.41* 

4 16 018.66 13 918.64 11 918.67 
5 17 400.16 15 300.16 13 300.18 
6 14 




















vr" 
J ~ 0 1 2 3 4 5 
O 11649.14 
1 13 137.49 
2 14558.24 12 743.35 
3 15911.16 14096.23 12 350.41 
4 15 380.22 13634.40 11956.15 
5 14 847.31* 13 169.12* 
6 15985.94 14307.64 12695.57 
7 15 366.35* 13 754.23* 12 207.29° 
& 16 336.30 14724.19 13 177.25 
9 15591.37 14044.37 
are 
w s y z Ve 
Hy 2196.13 65.80 —0.433 11 838.91 
HD 1905.17 51.70 0.522 —0.091 11 825.68 
D» 1556.64 34.51 0.287 0.04 11 813.40 
D::HD 0.81706 0.6675 0.550 —0.44 . 
p" 81675 0.6671 0.544 0445 — 


The values for He are obtained from the formula 
given by Richardson and Das.* They cannot be 
compared directly with the constants for HD 
and D, as they were calculated in a different way 
from less complete data. 

The ratio of the constants for D, and HD, 
which is given in the fourth row, is again ap- 
proximately equal to the theoretical p" values 
calculated from the mass ratios. The remarks 
about the agreement in the Fulcher bands! apply 
also here. In another paper’ it was shown that the 
agreement cannot be expected to be better, and 
that if a higher accuracy is desired, various 
corrections must be applied. In that paper also 
the shifts between the origins vy, of the band sys- 
tems were discussed. 

I wish to thank Dr. C. R. Jeppesen for his help 
with the measurements of the lines. 


sO. W. Richardson and K. Das, Proc. Roy. Soc. A122, 
688 (1929). 
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Many peculiarities first observed in the Fulcher bands of the hydrogen molecule can be 
explained by the interaction of the 3p°II with the 3p°= level. Such peculiarities are: irregular 
A-doublings, perturbations in the 3p*II and 3p*E state, and the absence of P and R branches 
for large values of the vibrational quantum numbers. This last phenomenon is due to the pre- 


dissociation of the 3p*II state. 





$1 


HE present investigation is concerned with 
the interaction of the 3° and 3p'II levels 
of the hydrogen molecule. The 39*II—-2s*> 
bands, the so-called Fulcher bands were the first 
known regularities in the He spectrum. Richard- 
son and Das! to whom the complete analysis of 
the Fulcher bands is due, noted that the Pand R 
branches of the bands showed marked irregulari- 
ties and were completely absent if the vibrational 
quantum number of the initial state exceeded the 
value three, whereas the Q branches had a per- 
fectly regular structure. 

The study? of the analogous bands of HD and 
D, led to an explanation of these phenomena as 
being due to the interaction of the 3° and the 
3p*II levels and the analysis of the 3p*S-—2s*d 
bands’ confirmed this view completely. The pres- 
ent paper is concerned with a closer investigation 
of this interaction. As examples of interaction 
between molecular levels which give rise to 
perturbations in which both levels are completely 
known are extremely rare, the present case is of 
considerable general interest, especially as the 
same kind of interaction can be studied for the 
three different isotopic molecules. 

The empirical material for HD and Ds is given 
in two previous papers”: * 4 and the data for He 
are taken from Richardson and Das.' 

The theoretical background for the treatment 
of such interaction has been summarized in a 
previous paper® to which frequent reference will 
have to be made as “‘theoretical part.” 

QO. W. Richardson and K. Das, Proc. Roy. Soc. A122, 
688 (1929). 

2G. H. Dieke and R. W. Blue, Phys. Rev. 47, 261 (1935). 

5 See the preceding article. 

‘ As more extended measurements have become available 
the analysis of the Fulcher bands of HD and D,» could be 


considerably amplified. Also the 3p°Z—>2s°E bands of He 
could be extended. 


§2. EmprricaL Facts 


It is perhaps most convenient to begin with a 
summary of the observed empirical facts which 
require an explanation. They are: 


1. All the Q branches are quite regular. 

2. The A-doubling of the 3p°II state depends in a very 
erratic way on the vibrational quantum number, and is 
quite abnormal for some vibrational levels. 

3. Whereas the dependence of the A-doubling on the 
rotational quantum number J is only slightly abnormal in 
most cases, there are pronounced perturbations for certain 
levels. 

4. The Hs, HD and Dz levels show the same general 
behavior, but the irregularities are found at different 
places for the three molecules. 

5. For V’>3 for Hz and HD and V’>4 for D2 the P 
and R branches are completely absent in the Fulcher bands. 

6. The 3p%S levels show also irregularities, but only the 
big perturbations are easily recognized, as there is no 
A-doubling which registers abnormalities with much 
greater precision. 


§3 


If we consider Fig. 1 which is a plot of the 
vibrational levels of the two electronic states, the 
explanation of many of the points mentioned 
above becomes apparent. In the first place we 
note that there are certain pairs of levels which 
lie very close together. Such pairs are 3p°=, V=4 
and 3p*II, V= 1 for He, 3p°=, V= 3 and 3p°11, V=0 
for HD, and 3p*2, V=5 and 3p*II, V= 1 for De. All 
these levels show big perturbations which makes 
it practically certain that the abnormalities 
listed in §2 are due to the interaction of the 3p*Z 
and 3p*II levels, and a closer investigation con- 
firms this. 

For the notation, etc., reference must be made 
to the theoretical part.® As the interaction is be- 
tween a II and & level it must be of type A, i.e., 


5G. H. Dieke, Phys. Rev. 47, 870 (1935). 
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Fic. 1. The vibrational levels of 3p°S and 3p*Il. The 
relative positions of the He, HD and D, levels are adjusted 
so that the vibrationless (electronic) level of 3p°= has the 
same value for all three molecules. 








the coupling of the two states is due to the rota- 
tion of the molecule. 

The Il state is always split up into a pll* 
and a PllI~ state® due to the rotation of the 
molecule (A-doubling) and of these only the pII* 
state can interact with the pX* levels. The pllI- 
level is the initial state for the Q branches and is 
completely unaffected by the interaction. This 
explains point 1 of §2, that all the Q branches are 
completely regular. 

If there were no interaction the pIl- level would 
coincide exactly with the pIl* levels, so that the 
amount of the A-doubling is a direct measure for 
the magnitude of the interaction. We have seen 
in the theoretical part that the interaction be- 
tween two levels is by far the greatest if they 


® Usually the pIl~ state is called pill, or pil. and the pit* 
state Pll, or plla. ; 


have the same vibrational quantum number, and 
we call the effect of this kind of interaction the 
regular A-doubling. If the vibrational quantum 
numbers differ, the interaction is much smaller 
and its effect on the energies is of any conse- 
quence only if the levels are close together. This 
gives us the irregularities in the A-doubling. 


§4. THE REGULAR A-DOUBLING 


According to (16) and (24) of the theoretical 
part the regular doubling is in first approximation 
given by 

2B,7J(J+1)/ 4, 


in which 6 is the distance between the p> level 
and the corresponding II level with the same 
rotational quantum number J and the same vi- 
brational quantum number IV. As the equilibrium 
distances are different for the two interacting 
levels, a better approximation is obtained if we 
substitute B,,= /RBR’rdr for B,= { RBRr dr. 
(R is the vibrational wave function.) This follows 
from (17) if the second small and uncertain term 
is omitted. This correction is very small in our 
case. With harmonic binding we get, e.g., 


Boy= (2a ‘1 +a?) le? /2(1+0%) Bo 
in which 
a=(w,'/w,)', b?=42*cw,’u(ro—ro’)?/h 
and the values of the constants are for De 
3p*Il = ro=1.047-10-°,  w, = 1678.22 
3p*S rp’ = 1.096-10-8,  w,’= 1556.64. 


This gives for De Boop=0.941 Bo, and similarly 
small corrections for the two other molecules. 

The regular doubling for the V=0 vibrational 
level is found to be in this way 0.321 J(J+1) for 
He, 0.180 J(J+1) for HD, and 0.078 J(J+1) 
for De. 

In Figs. 2, 3 and 4 these calculated regular 
doublings are compared with some of the ob- 
served values. Table I gives all the observed 
doublings.’ It is seen that the order of magnitude 


7 Asa pil* and a pli” level with the same J never can 
combine with the same lower state, it is not possible to 
obtain the values of the doubling directly from the observa- 
tions by using only the combination principle. However, 
as the structure of the lower (2s*S) state is very regular, 
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TABLE I. A-doubling of the 3p*I1-level. 

Hy HD D» 
J Vv=0 1 2 3 0 1 2 3 0 1 2 3 4 
I 0.46 —0.31 0.92 1.89 —0.04 O48 0.29 0.67 0.13 0.09 O18 O13 0.23 
2 125 —-163 2.68 4.43 —0.43 1.34 0.84 1.89 0.35 0.20 0.62 0.37 0.74 
3 2.54 — 6.50 4.76 6.46 — 3.26 2.45 1.46 3.34 0.50 0.21 1.23 0.70 1.36 
4 3.90 6.87 8.28 +9.06 3.83 1.95 4.88 1.00 —-O17 1.97 0.98 2.04 
5 5.19 8.36 7.67 5.33 1.75 6.10 149 —3.08 2.71 0.75 2.62 
6 6.06 7.49 6.59 6.55 1.86 +8.60 3.54 —0.87 2.96 
7 2.18 6.22 


of the observed doubling is well represented by 
the theoretical curves,* but that the individual 
values show considerable systematic deviations 
from the curves, positive in some cases, negative 
in others. 


$5. IRREGULAR A-DOUBLING 

These deviations must be due, for the greatest 
part,’ at least, to the interaction of the fll level 
with that pY level which is its closest neighbor. 
We take first the case where the interacting levels 
are not so close that we have typical perturba- 
tions. The interaction of two levels with different 
vibrational quantum numbers is given in first 
approximation by (17) of the theoretical part in 
which we can omit again the uncertain second 
part. 

The quantity B,,.. can be calculated easily if we 
make use of the well-known properties of the Her- 
mitian polynomials.'? We get, e.g., Bos=0.02By 
for De. This is much too small to account for 
the empirical deviations of the V=0 level of 39II 
from the theoretical curve (Fig. 2) which would 
require Bo, to be of the order of magnitude 0.1 Bo 
to 0.2 Bo. 


the differences between successive rotational levels of it 
can be calculated with any desired accuracy. In fact it is 
necessary only to calculate the difference between the 
J=0 and J=1 levels of 2s°E with the help of the term 
formula, and all the rest including the A-doublings can be 
calculated with the help of the combination principle only. 
The values in the table were obtained in this way. 

‘It may be good to emphasize that no adjustable con- 
stant was used for the calculation of the theoretical curves. 

® There are, of course, other reasons why there should be 
deviations of the observed A-doublings from the theoretical 
curves. The formula used is only a first approximation and 
can be expected to hold only for small J. The figures show 
that for all levels the observed values for the larger values 
of J seem to have a tendency to be too low, and this kind 
of deviation must be due to the insufficiency of the formula 
used. 

10 The calculations necessary for this are of the same kind 
as those used for the calculations of band intensities in 
electronic band systems. As the calculations do not give 
the right magnitude anyway they are omitted here. 


The explanation for this discrepancy must be 
sought in the fact that the calculations are made 
under the assumption that the binding is har- 
monic.'' They are, however, very sensitive for 
deviations from the harmonic law. The calcula- 
tions could be carried out also with a general an- 
harmonic law of force but would be much more 
cumbersome, and the present state of affairs 
would hardly justify going into them. Possibly 
also the neglected term in (17) may be appre- 
ciable. We hope to come back to this ina later 
paper but restrict ourselves now to some quali- 
tative considerations. 

Whether the A-doubling of a given Il level is 
diminished or increased depends on whether the 

Y level which is its nearest neighbor lies just 
above or below the PII level. There is no simple 
law which governs the distance of a given Pll 
level from the closest pL level, therefore the 
relative position of such pairs can be considered 
to be regulated more or less by chance. This ac- 
counts for the erratic way in which the irregulari- 
ties of the A-doubling depend on the vibrational 
quantum number (point 2 of §2). Also there is 
absolutely no connection between the relative 
arrangement of such pairs in He, HD and Dp 
which can be verified by a glance at Fig. 1. This 
explains point 4 of §2. 

If the 3p> level lies just above the neighboring 
3pll level, the effect of this will be that the 3fII 
level will be depressed, i.e., the A-doubling 
diminished. We have such a state of affairs, e.g., 
for V=0 in De, V=2 in HD and V=0 in Haz. 
These are the levels represented by an x in the 
Figs. 2-4, and the A-doubling lies in all these 
cases below the theoretical value, as it should. 

'! Some of the effects of the deviations from the anhar- 
monic binding have been taken into account; e.g., if we 
use B,=B,.—a(V+}), the empirical value of @ is used. 


If the binding is regarded purely harmonic, @ would have 
the same order of magnitude but opposite sign. 
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Fic. 2. A-doubling for 3p°II of De. Curve: calculated 
regular A-doubling for 1’=0. x Observed values for =0, 
o for V=2. 
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Fic. 3. A-doubling for 3p*Il of HD. Curve: calculated 
regular A-doubling for 1 =0. x Observed values for |'=2, 
o for V=1. 
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Fic. 4. A-doubling for 3p*II of He. Curve: calculated 
regular A-doubling for /=0. x Observed values for )=0, 0 
for V=2. 


It may even occur that the deviation is larger 
than the regular doubling in which case the total 
doubling becomes negative. 

The opposite state of affairs exists if the 3p= 
level lies just below the neighboring 3/II level as 


is the case for V=2 in De, V=1 in HD, and 
V=2 in He. Then the A-doubling should be 
larger than the regular doubling and that is 
found to be the case (see Figs. 2-4). In judging 
the agreement it must be borne in mind that the 
drawn out curve is the regular A-doubling for 
V=0, whereas the regular A-doubling for V=1 
and |’=2 is lower (about 25 percent lower for 
the V=2 level of Hz). The figures show that the 
deviations from the regular curve become rela- 
tively much larger for larger J when the doubling 
is below its regular value and the opposite is 
true when it is above its regular value. 

The reason for this is that when the Il level 
lies below the p> level the distance 6 between 
rotational levels with the same J decreases with 
increasing J, as the spacing of the pIl rotational 
levels is larger than the spacing of the p> levels. 
This means that the higher levels are relatively 
more affected because of their closer proximity 
than the lower ones. If the pI state is above the 

> state just the opposite takes place. All this is 
in agreement with the experimental facts. 


$6. TRUE PERTURBATIONS 


We come now to the cases where the p> levels 
lie so closely above the PII levels that the inter- 
action becomes considerable and we have true 
perturbations. Fig. 1 shows that this is the case 
for V=1in Hzgand Dz. and V=0in HD. The data 
are most complete for the pair = 1 of 3p*II and 
V=5 of 3p*> of De and therefore we shall con- 
sider this case more in detail. The two other 
cases mentioned are of the same kind. 

The magnitude of the perturbation can be 
obtained with great accuracy in the following 
way. The total A-doubling is the sum of the 
regular A-doubling which can be calculated in the 
way indicated in §4 and the perturbation due to 
the V=5 level of 3p . By subtracting therefore 
the calculated regular A-doubling from the total 
doubling we obtain the amount of the perturba- 
tion. According to theory the perturbation of the 
V=5 state of 3p> is of the same magnitude but 
opposite sign. By subtracting resp. adding the 
amount of the perturbation to the pair of the 
perturbed levels we obtain the unperturbed 
states.” The difference between pairs of such 








'2 Strictly speaking the interactions of the 3pII V = 1 level 
with the 3p2V'=1 and the 3pXV=5 levels are not inde- 
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Fic. 5. The solid curve represents the calculated per- 
turbations for the 3p*Il,  =1 level of D2. o are the observed 
perturbations. The dotted curve gives the distance between 
corresponding rotational levels of 3p°=, V=5 and 3/*II, 
V=1. (The scale for this curve is ten times smaller than 
that indicated by the figures on the left.) 
levels with the same J is represented by the 
dotted curve in Fig. 5. The o mark the observed 
perturbations (on a tenfold magnified scale). 
It is seen that the perturbations become rapidly 
larger as the distance between the levels de- 
creases and that they change their sign when the 
two levels cross over. The solid curve represents 
the perturbations calculated according to (13) 
and (18) of the theoretical part, the constant a 
being chosen in such a way as to give the best 
agreement for small values of J. It is seen that 
the agreement is very good for small J but that 
after the break there are discrepancies much 
larger than the errors of measurement. They are 
due to the approximate character of our calcula- 
tions and undoubtedly the agreement could be 
made much better by going to higher approxima- 
tions. 

The value of the constant a is 4.9 from which 
it follows according to (17) that B,;= 1.5. That is 
the same order of magnitude as was found for 
Bos in §5. 

In all cases the perturbations of the 3p state 
pendent as they are caused by the same interaction force. 
It can be easily shown, however, that they can be treated 
as if they represented separate interactions, if first the 
regular A-doubling is taken into account and the new 
levels, shifted from their original position by the amount of 


the regular A-doubling are taken as the unperturbed states 
for the interaction between the close pairs. 
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are of the same magnitude and opposite sign as 
those for the 3pII state with which they interact. 
Only they cannot be ascertained with the same 
precision as those of the pIl states as the unper- 
turbed states are unknown whereas for the PII 
levels they are exactly given by the unaffected 
pll~ states. 

There is a perturbation of smaller magnitude 
than those mentioned above in the 3p2V=7 
level of De. As Fig. 1 shows, this perturbation 
cannot be due to the interaction with a 39ll 
state. It is probable that it is caused by the 
interaction with the 4p* level which, however, 
has not yet been established empirically with 
sufficient certainty to verify this unambiguously. 


$7. PREDISSOCIATION 


The vibrational levels of the 3p state are 
discreet only up to a certain limit. Beyond this 
limit which is indicated in Fig. 1 the molecule 
is dissociated into a normal and an excited atom 
in the second quantum state, and the energy 
levels are continuous. The nature of the wave 
functions is not changed in such a way as to 
affect materially the interaction with the 39II 
state. This means that the continuous levels inter- 
act with the neighboring 3II levels in essentially 
the same way as the discreet levels. A character- 
istic feature of such interaction is that the mole- 
cule oscillates back and forth between the two 
interacting states. As the 3p2 state is dissociated 
the 3pIlI state will also be dissociated by the inter- 
action (predissociated) i.e., the 3pII* levels a- 
bove the dissociation limit of the 3p state will 
be destroyed as discreet levels. The bands 
originating from these levels will therefore be 
absent.'® The Fig. 1 shows that this must be the 
case for V>3 in Hz and HD and for V >4 in Dz, 
and this agrees exactly with the experimental 
observations (point 5 of §2). 

All this applies only to the 3pII* level which 
gives the P and R branches, as the 3II~ level 
which gives the Q branches is not affected at all 
by the interaction. This explains why only the 
P and R branches are absent, whereas the Q 
branches show no changes. 

In a previous paper’? this phenomenon was 
used to determine the heat of dissociation of He. 


18 The continuous spectrum which should come in its 
place, is too weak to be noticed. 
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The Magnetic Anisotropy of Crystals of Sn, and of Sn with Added Sb, Cd, or Ga* 


HAROLD J. HoGe, Sloane Physics Laboratory, Yale University 
(Received June 27, 1935) 


The principal magnetic susceptibilities of single crystals 
of white (tetragonal) tin have been determined. Crystals 
were grown by the traveling furnace method. Perpendicular 
to the tetragonal axis the susceptibility (per unit volume) 
is x, =0.197 X10~; parallel to the axis it is ky =0.176 107°. 
The ratio xi /x, is strongly influenced by the addition of 


small amounts of certain elements. For tin containing no 
added element this ratio is 1.120. Adding 0.60 percent Sb 
the ratio is 1.200; 0.76 percent Cd, 0.923; 0.89 percent Ga, 
0.897. The relation between the change in this ratio and 
the number of valence electrons of the added element is 
compared with available data for bismuth. 





HAT the crystal structure of tin controls 
its magnetic properties is evident from the 
fact that white tin is paramagnetic, grey tin 
diamagnetic. It therefore seemed important to 
investigate the magnetic anisotropy of the non- 
cubic form (white tin) and to modify this aniso- 
tropy in any feasible way. The method recently 
developed by Goetz and Focke' and applied by 
them to bismuth (diamagnetic) were obviously 
suitable with such minor changes as were made 
necessary by the lower absolute values of the 
principal susceptibilities of tin. The investigators 
just referred to grew single crystals suitable for 
magnetic measurements by the Gouy method 
and found how the two principal susceptibilities 
were affected by impurities. 

Impurities which can be entrapped in crystals 
are of greater interest than those which cannot. 
Whether there is true solution in all such cases is 
not yet clear, but it appeared reasonable to try 
first impurities known to dissolve in tin. Of these 
it has been possible to study only antimony and 
cadmium. Gallium was added to the list because 
of its low melting point, peculiar crystal struc- 
ture, and favorable position in the periodic table. 


EXPERIMENTAL PROCEDURE 


The tin selected for definitive experiments was 
a particular lot of Baker's “‘C. P. Analyzed”’ con- 
taining As 0.0000 percent; Pb 0.003 percent; 
Fe 0.002 percent; Zn 0.000 percent; Cu 0.000 
percent. Since not over one percent of Sb, Cd, 
or Ga was added, ordinary commercial grades of 
these elements could be used. Alloys were pre- 
pared by weighing and were melted under hy- 

* Part of a dissertation presented for the degree of Doctor 
of Philosophy in Yale University. 

1 A. Goetz and A. B. Focke, Phys. Rev. 45, 170 (1934), 
and preceding papers. 


drogen in Pyrex. The low range of solid solubili- 
ties made it inconvenient to obtain a series of 
the desired concentrations by diluting a richer 
alloy. After thorough mixing a stick 7 inches 
(18 cm) long and 0.1 inch (0.254 cm) square was 
cast in a two-piece slate mold. 

This polycrystalline casting was placed in a 
snugly fitting graphite box (Fig. 1). One end of 


























Fic. 1. Graphite box in which crystals were grown. 


the box was closed by a steel plug. At the other 
end a seed crystal could be brought into contact 
with the casting. 

This assembly was then placed in a horizontal 
position inside a crystal growing apparatus simi- 
lar to that described by Hasler.2 The tubular 
furnace which traverses the casting was made 
quite short (4 cm) and the temperature gradient 
on the exit side was made steeper by directing 
air jets against the glass wall of the tube contain- 
ing the graphite box. Cooling of the furnace by 
this jet system was prevented by an asbestos 
screen. Preliminary tests showed that a steep 
temperature gradient is essential to success, and 
that 2 millimeters per minute is a suitable furnace 
speed. The use of a short furnace ensures against 
too much change in concentration of soluble im- 
purities. Hydrogen was used to prevent oxida- 
tion. The temperature gradient was further in- 
creased by causing the hydrogen to flow past the 
crystal in the direction in which the furnace 
moved. That the temperature gradient requires 


2M. F. Hasler, Rev. Sci. Inst. 4, 656 (1933). 
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special control was shown by the fact that it was 
possible to initiate the growth of crystals with 
the tetragonal axis at a large angle to the axis of 
the casting by the use of a steel plug with its end 
beveled at a sharp angle. This procedure was 
suggested by the experimental results of Palibin 
and Froiman,’ and was a convenient way to get 
seed crystals of the desired orientation, which is 
with the tetragonal axis perpendicular to the 
axis of the specimen. 

The castings thus prepared were etched in 
hydrochloric acid solutions properly modified to 
prevent a too rapid attack. The etching showed 
at a glance whether a single crystal had been ob- 
tained, and if so, permitted its orientation to be 
determined by optical means. The most impor- 
tant etch planes were found to be {101}! and 
|301},4 these indices referring to the smaller of 
the two commonly chosen unit cells. 

Susceptibilities were measured by the Gouy 
method, in which the force on a long specimen of 
uniform cross section is given by F=3(k—ko) 
A (II? —IJ1,?). Here x and ko are the susceptibilities 
(per unit volume) of the specimen and of the 
surrounding medium (air), 77 and //) are the field 
strengths at the two ends of the specimen, and A 
is its cross section. Since x) is known and //;? is 
negligible for such long specimens, the quantities 
to be measured are F, A, and //. 

The force F was measured by counterpoising 
the specimen on a sensitive microbalance, and 
increasing the counterpoise as necessary when the 
field was applied. Readings were taken to the 
nearest 0.01 milligram. High sensitivity was re- 
quired because of the small magnitude of F, 
which was of the order of 3 milligrams. 

The fields 7/7 and //) were measured with a 
search coil attached to a pendulum and connected 
to a low resistance (16 ohms) Leeds and Northrup 
type HS-2285e galvanometer which was cali- 
brated with a standard mutual inductance. Field 
strengths // of about 24,000 gauss were employed, 
in which case //) was of the order of 650 gauss. 

The cross section of the specimen was found 
by lowering it a known distance into a tube of 
ethyl alcohol and measuring the rise in the liquid 


8’ P. A. Palibin and A. I. Froiman, Zeits. f. Krist. 85, 322 
(1933). 

' Details of the optical apparatus used and of the identi- 
fication of the etch planes will be given elsewhere. 
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level with a traveling microscope. This rise was 
most easily measured in a connecting side tube of 
small diameter. The apparatus was calibrated by 
immersing pieces of drill rod of known cross 
section. 

Sections 3.5 inches (9 cm) long were cut from 
the good crystals, selected to give the most uni- 
form cross section in the non-homogeneous part 
of the field. The bar to be cut was laid in the 
jaws of a pair of tin snips and allowed to fall into 
the hand when cut, so that all strains would be 
localized at the ends of the specimen. 

The crystal hanger contained a device which, 
when the balance was arrested, allowed the 
crystal to be lifted, rotated about its axis, and 
reseated without lateral displacement. A crystal 
of the proper orientation could thus be set with 
its tetragonal axis either parallel to or perpendicu- 
lar to the field. Since all directions perpendicular 
to the tetragonal axis are equivalent, the suscep- 
tibility in any arbitrary direction in a tetragonal 
crystal is given by x=x«,, cos? a+x«, sin? a, where 
a is the angle between the applied field and the 
tetragonal axis. Measurement of «,, and «, there- 
fore gives complete information about the mag- 
netic properties of any crystal. Since the poly- 
crystalline susceptibility of tin is temperature 
independent® it seemed unlikely that small 
temperature fluctuations would introduce any 
error, and all measurements were made at room 
temperature. 


RESULTS 


For pure tin: x, = 0.197 X 10~, «,, =0.176 X 10°*, 
k,/k,,=1.120. Taking the density of tin as 7.3 
the mass susceptibilities are: x, =0.0270 10~°, 
X;, = 0.0241 X 10~®. The Voigt-Thomson® relation 
gives for an ideal polycrystalline aggregate: 
x =3(2x,.+x,,) =0.0259 X 10°, which is in good 
agreement with previously published results.? 

The experimental data are given in Table I. 
Figs. 2 to 4 show the variation of the two princi- 
pal susceptibilities for each impurity as a function 
of its concentration. Fig. 5 shows the ratio of these 
susceptibilities (x,/x,,) for each impurity as a 
function of its concentration. The effect of Sb is 
to decrease both «x, and x, «,/x,, increasing. 


°H. Endé, Sci. Rep. Tohoku Imp. Univ. 16, 201 (1927). 
®°W. Voigt, Lehrbuch der Kristallphysik (1928), p. 960. 
? International Critical Tables, Vol. 6, p. 355. 
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TABLE I. Susceptibility data on alloys if single crystals of tin. 








CRYSTAL ALLOY WT. ORIEN- Cross FIELD 1? —Hy 
NO. % TATION SECTION Forcre.* Forcru* (gauss) x10°% 
19 0.09 Sb 89.9° 0.06493 cm?® 3.371 mg 2.920 mg 24,700 6.0967 
3.375 2.925 24,700 6.0967 

20 .26 Sb 86° 06278 3.135 2.660 24,400 5.9494 
3.138 2.663 24,400 5.9494 

18 .09 Sb 89.9° .06278 3.143 2.705 24,304 5.9026 
3.130 2.703 24,378 5.9386 

17 Stock B 90° .06329 3.198 2.795 24,353 5.9265 
3.193 2.788 24,304 5.9026 

13+ Stock B 89.9° -.06046 3.265 2.865 24,304 5.9026 
& Stock B 87° .06193 3.103 2.745 24,378 5.9386 
21 0.76 Cd 85.5° .06244 2.823 3.105 24,378 5.9386 
2.830 3.108 24,378 5.9386 

23 .50 Cd 87° .06391 3.033 3.104 24,333T 5.9167 
16 Stock B 89° .06123 3.130 2.725 24,324t 5.9123 
3.135 2.733 24,324f 5.9123 

27 0.25 Cd 89° 05631 2.775 2.683 24,324 5.9123 
2.775 2.675 24,294 5.8978 

28 .12 Cd 88° 06391 3.153 2.888 24,294f 5.8978 
3.143 2.863 24,294 5.8978 

31 .60 Sb 88° -06262 2.963 2.383 24,259 5.8808 
2.945 2.360 24,2597 5.8808 

32 42 Ga 87° .06193 3.075 3.035 24,288 5.8948 
3.058 3.020 24,288t 5.8948 

33 89 Ga 88° -06262 3.133 3.555 24,298 5.8997 
3.128 3.553 24,298T 5.8997 
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* After subtracting suspension correction of 0.1 mg. 





t After adding air correction of 0.029 x 10~, 


+ Fluxmeter calibrated at time of taking this reading. 
+ Not averaged. 


1.12 
1.12 
1.029 
1.032 


1.077 
1,082 


1.198 
1.201 


1.011 
1.011 


8971 
8961 








When Cd is added «x, decreases and x,, increases; 
k,/«, therefore decreases. At 0.4 percent cad- 
mium «,=x,, and the crystal is magnetically iso- 
tropic. When Ga is added x, changes very little 
while x,, increases rapidly; x,/«,, decreases as it 
does when Cd is added. 

It may be seen that Cd and Ga, the two ele- 
ments having fewer valence electrons than tin, 
cause a decrease in «,/x,,. On the other hand Sb, 
which has one more valence electron than tin, 
causes an increase in this ratio. The significance 
of this finding is increased by the fact that x, /x,, 
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FG. 2. Susceptibilities of tin containing antimony. 
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Fic. 3. Susceptibilities of tin containing cadmium. 














| Wh 6a | | | 





GCdt#éiewsts tits se «és & 


Fic. 4. Susceptibilities of tin containing gallium. 
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is not subject to errors made in measuring either 
field strength or crystal cross section, since these 
variables cancel out. Therefore, although the 
probable error in «x may be as high as 2 percent, 
the probable error in x,/x,, is considerably less. 
The precision in measurement of x was checked 
by determining the mass susceptibilities of dis- 
tilled water and fused quartz. Results were, for 
water, x= — 0.738 10-* compared with the ac- 
cepted value of —0.720X10-*; and for fused 
quartz x= —0.420 10-* compared with a value 
of —0.42X10-® obtained by Dr. C. T. Lane in 
these laboratories but with other apparatus. 


DISCUSSION 


It is interesting to compare Fig. 5 with the cor- 
responding data given by Goetz and Focke for 
bismuth. With tin, elements toward the alkali 
side of the periodic table decrease «,/x,,, while 
those toward the halogen side increase this ratio. 
The exact opposite was found for bismuth. How- 
ever, when account is taken of the fact that the 
susceptibilities of tin are both positive, while 
those of bismuth are both negative, this difference 
is perhaps to be expected. Assume for definiteness 
that an added impurity contributes diamagnet- 
ism parallel to the principal axis and nothing at 
all perpendicular to it. Then in the case of tin we 
are adding a negative quantity to the smaller of 
two positive quantities, while in the case of bis- 
muth we are adding it to the smaller of two nega- 
tive quantities. The ratio for tin will be increased 
and for bismuth will be decreased. Something 
more complicated than this actually occurs, for 
the magnitudes of x, and x,, both change. 

Goetz and Focke explained their results with 
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F1G. 5. «,/«, ratios for tin containing impurities. 
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the aid of the Ehrenfest-Raman*: *: '":" large 
electron orbit hypothesis, which had previously 
proved useful in explaining the abnormally high 
diamagnetism of bismuth, antimony, and graph- 
ite crystals, and the diminution in these values 
when the crystals were made submicroscopic or 
melted. Assuming the impurity to collect in cer- 
tain crystallographic planes, Goetz and Focke 
pointed out that these planes might act as reflect- 
ing walls to electrons moving in large orbits. A 
reorientation of these orbits would result, and 
since they are presumably large enough to sur- 
round many atomic nuclei, the diamagnetic sus- 
ceptibility in certain crystallographic directions 
would be considerably increased at the expense 
of that in other directions. 

Processes of this kind are adequate to explain 
some of the effects of impurities on tin. The five 
planes having the greatest atom densities are, in 
order: {100} {110} {112} {101} {210}. By as- 
suming impurities to be concentrated in {110} 
planes, the changes in «x, x, may be explained. 
Planes of the form {110} are suggested as more 
probable than {100} by the fact that several of 
the crystals were covered after etching by fine 
lines which corresponded to traces of {110} 
planes. These traces were most regular and con- 
spicuous on the lateral faces of the crystal, and 
were best developed in the direction most nearly 
parallel to the sloping isothermal surfaces in the 
freezing melt. That is, one of the planes (110), 
(110) was conspicuous; the other was not. It is 
difficult to explain these lines without attributing 
them to impurities, particularly since they were 
most pronounced in crystals containing the 
largest amounts of impurities. Although the lines 
were so widely spaced as to be visible without 
magnification, they certainly do not make less 
plausible the existence of planes of impurity 
spaced more or less uniformly at much smaller 
intervals throughout the crystal. A considerable 
amount of work indicates the presence of some 
sort of secondary structure in metals, so that 
there is no longer much question of its existence. 
We are now more interested in the sizes, causes, 
and regularity of these secondary patterns. 


8 P, Ehrenfest, Physica 5, 388 (1925). 

® P. Ehrenfest, Zeits. f. Physik 58, 719 (1929). 
10C, V. Raman, Nature 123, 945 (1929). 

1 C, V. Raman, Nature 124, 412 (1929). 
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In the present work the magnitude of the ob- 
served changes (20 percent in x,/x,, for less than 
1 percent added impurity) makes it possible to as- 
sociate the changes not with the mere presence 
of the impurity but with the effect of its presence 
on the rest of the crystal. Goetz and Focke point 
out that a distribution of impurity in widely 
separated planes could be expected to produce a 
larger effect than a uniform volume distribution. 
Furthermore, a distribution in planes accounts 
much more readily for anisotropic changes in 
susceptibility. 

To explain why an element having more 
valence electrons than tin should have an oppo- 
site effect to one having less, we may assume 
either that different impurities go into different 
planes, or that they go into the same planes but 
affect the crystal differently. 

One limitation of the large orbit hypothesis 
should be pointed out. A mere rearrangement of 
existing orbits cannot affect the calculated ideal 
polycrystalline susceptibility, for any increase 
in x, would result in a corresponding decrease in 
k,, and vice versa. As can be seen from Figs. 2 to 4, 
the corresponding polycrystalline susceptibility 
(which would be given by a line 1/3 the distance 
from x, to x,,) does not remain constant. 

If all added impurities form reflecting walls, 
then to explain changes in «x,/«,, we must assume 
that Cd or Ga collects in planes parallel to the 
axis while Sb collects in planes perpendicular to it. 
For example, the addition of Cd decreases x,/,,. 
If Cd collects in any set of planes parallel to the 
tetragonal axis, we would expect an increase in 
the area of large orbits projected on these planes, 
and an increased diamagnetism perpendicular to 
the axis. The net positive value of x, would there- 
fore be lowered relative to x, and «,/x,, would 
decrease, as is actually observed. 

One limitation of the large orbit hypothesis 
has already been mentioned. Another appears 
when we assume that impurity collects in planes 
parallel to the tetragonal axis, for unless condi- 
tions of growth favor one plane more than an- 
other there is no reason to expect impurity to 
collect in, say, (110) and not to collect in (110). 
If it collects in both the planes of form {110} the 
crystal is divided by walls of impurities into nar- 


row rectangular prisms, and no large orbits can 
exist at all. We may avoid this situation only by 
assuming that all impurities collect in (001) 
planes, which divide the crystal not into prisms 
but into layers. But if we do so there is still the 
difficulty of explaining the opposite effects of 
different impurities. 

Honda and Shimizu” have suggested that in- 
crease in lattice constants may produce diamag- 
netism. In this connection the results of Bowen 
and Morris Jones" have some bearing on the 
present work. By x-ray investigation they found 
that when Sb is added to tin there is a gradual 
expansion of the tin lattice up to the limit of 
solid solubility, with no change, however, in the 
ratio of the lattice parameters. Fig. 2 shows a 
decrease in paramagnetism, as this line of argu- 
ment would predict. It would be interesting to see 
if Ga, which increases the paramagnetism of tin, 
simultaneously decreases its lattice constants. 

In comparing the effects of impurities on tin 
with the effects on bismuth it should be pointed 
out that there is little evidence in the present 
work for a critical concentration of impurity cor- 
responding to the complete filling of a set of uni- 
formly spaced planes with a one- or two-atom 
layer of impurity. The nearest thing to such evi- 
dence is the dip in the Cd curves (Fig. 3) near the 
origin. Since all the critical concentrations re- 
ported by Goetz and Focke lay below 0.1 percent, 
it is not impossible that minima in these curves 
occur at lower concentrations than any yet used. 
The stock tin was not, however, pure enough to 
warrant studies at such low concentrations. It will 
be observed that the changes in «,/x, for tin are 
much smaller and much more nearly linear than 
those for bismuth over the same range in concen- 
tration. The approximate linearity may be due 
merely to the smaller range of «,, «,,, and of their 
ratio, in the present case. 

In conclusion the author wishes to express his 
appreciation to Professor L. W. McKeehan, 
director of the research, for his ready advice and 
suggestions. He is also indebted to Dr. C. T. 
Lane and Dr. A. B. Focke for much valuable aid. 


12K, Honda and Y. Shimizu, Nature 135, 108 (1935). 
18 E. G. Bowen and W. Morris Jones, Phil. Mag. 12, 441 
(1931). 
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The expansion coefficients were measured at 5° intervals from — 115°C to — 160°C and mean 
values found from —160°C to —190°C. The upper temperature limit was set by a recrystal- 
lization process which set in between —115°C and —110°C. a, was found to be larger and 
to increase more rapidly with the temperature than a,. The volume coefficient plotted in 
conjunction with Carpenter and Oakley’s measurements give values of the volume coefficient 
from the melting point to — 160°C where the coefficient apparently begins to decrease rapidly 


with the temperature. 





INTRODUCTION 


HE properties of single crystals and the 

variation of those properties with the tem- 
perature are of considerable importance in theo- 
retical work on the solid state. The only previous 
measurements of the principal expansion coeffi- 
cients of mercury are those of Griineisen and 
Sckell' who measured only the mean values of 
the expansion coefficients over the large tempera- 
ture interval from —190°C to —78°C and were 
unable to find how the coefficients vary with the 
temperature. Even the volume coefficient has 
not been satisfactorily measured. Dewar? found 
the value 8.87x10-° from —188.7°C to the 
melting point; Grunmach,’ 12.310-° from 
—80°C to the melting point; and Bridgman‘ 
obtained the values 12.5X10-° and 16.5x10~° 
just below the melting point quite incidentally 
in the course of an investigation on the change 
of volume on freezing. Carpenter and Oakley°® 
measured the volume expansion coefficient over 
the range —90°C to the melting point, but their 
measurements were not accurate or extended 
enough to show other than a linear variation 
with the temperature. 

The present research was begun in 1932 in an 
attempt to find how the volume coefficient varies 
with the temperature by measuring the linear 
expansion coefficients of polycrystalline rods. 
However, the results were so erratic that it 
became evident that it would be necessary to 


1E. Griineisen and O. Sckell, Ann. d. Physik 19, 387 
(1934). 

2 J. Dewar, Proc. Roy. Soc. A70, 237 (1902). 

3... Grunmach, Physik. Zeits. 3, 134 (1902). 

4P. W. Bridgman, Proc. Am. Acad. 47, 343 (1911). 

5L. G. Carpenter and F. H. Oakley, Phil. Mag. 12, 511 
(1931). 


use single crystals and measure the principal 
expansion coefficients. 
CRYSTAL STRUCTURE 

Mercury crystallizes in the trigonal system® 
so that it is necessary to measure only the 
coefficients parallel and perpendicular to the 
trigonal axis.’ The coefficient in a direction 
making an angle @ with the trigonal axis is given 
by the relation: 


ag= ay, COS” +a, sin® 6 
=a,,+(a,—a,,) sin? 6, (1) 


so that ag is linear with sin? 6. The angle @ which 
the trigonal axis makes with the axis of the 
single crystal rod will hereafter be called the 
orientation. If the coefficients of a number of 
crystals with different known orientations are 
measured and plotted against sin® 6, the inter- 
cepts of the straight line through them at 
sin? 62=0 and sin? @=1 will be a, and ay, 
respectively. 

The specific resistance in a direction making an 
angle @ with the trigonal axis is given by the 
relation :* 

poe= py, Cos? 6+ p, sin® @ 
sin? 6= (pe— py)/ (Ps — Pu)- (2) 


This offers a means of determining sin*® @ since 
Sckell® has measured p,, and p, at certain known 
temperatures, the values being later corrected by 
Griineisen and Sckell.'’ Sckell’s curves show that 


5M. C. Neuburger, Zeits. f. anorg. allgem. Chemie 212, 
40 (1933). 

7W. Voigt, Lehrbuch der Kristallphysik, p. 290. 

SW. Voigt, Lehrbuch der Kristallphysik, p. 344. 

®O. Sckell, Ann. d. Physik 6, 932 (1930). 

1 E. Griineisen and O. Sckell, Ann. d. Physik 19, 387 
(1934). 
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p,, and p, are sufficiently linear in the liquid air 
region to allow interpolation from the two 
temperatures at which he made measurements in 
that region to any other temperatures in the 
region. Thus, measurements of pg and the tem- 
perature allow sin? @ to be calculated from 
Eq. (2). 

Sckell found that p,, was smaller than p, and 
Eq. (2) shows that p,, should therefore be the 
smallest possible value of p. In the present work 
three crystals were found which had smaller 
values of p than that given by Griineisen and 
Sckell by 0.2, 0.3 and 0.5 percent. Although only 
the last is outside the limit of experimental 
error it appears quite probable that the value of 
p,, given by Griineisen and Sckell is still high, 
especially since their crystal No. 13 was also 
lower by 0.4 percent than the value they gave 
as most probable. 


EXPERIMENTAL 


Preliminary work with single crystals showed 
that they frequently recrystallized so that it was 
necessary to use a method whereby the specific 
resistance and hence the orientation could be 
determined both before and after the expansion 
measurement. This could be done if the crystals 
were grown in the form of cylindrical rods of 
uniform cross section. A method was adopted 
which was quite similar to that employed by 
Griineisen and Sckell. 

Plaster molds, 30 cm long by 0.8 cm inside 
diameter, were made by clamping a _ thinly 
paraffined drill rod on the axis of a glass tube 
which was thickly paraffined on the inside. 
Patching plaster was poured between the tube 
and the rod and the air bubbles pumped out 
with an aspirator. After the plaster set, the mold 
was removed from the glass tube by melting the 
paraffin and the rod removed from the center of 
the mold in the same manner. A mold made in 
this way was straight, circular, of uniform 
diameter, and the inside had an almost polished 
appearance. Deep grooves were made length- 
wise of the mold so that it could be broken free 
from the crystal. Pyrex tubes ending in various 
types of tips were waxed into the bottoms of 
the molds. In some the tubes were drawn down 
to straight capillaries 0.75 mm inside diameter 
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and 4 cm long, in others the tubes were melted 
down to thick-walled capillaries 0.5 mm to 1.5 
mm inside diameter below which bulbs of from 
2 to 4 cc capacity were blown. The third type 
was similar to the second but had one or two 
tungsten wires sealed into the bottom of the 
bulb and extending up to about the center of 
the bulb. The type of tip seemed to have very 
little influence on the percentage of single crystals 
obtained or on the orientation. The orientations 
were apparently random. 

The molds were filled in a vacuum by placing 
a glass reservoir containing the mercury at the 
top end of the mold and sealing both mold and 
reservoir in a larger glass tube which was then 
evacuated to fluorescence for several hours, the 
mercury being heated in the mean time to drive 
out water vapor and gases. The mercury was 
then poured into the mold by tipping the 
apparatus. 

The crystals were grown either by lowering 
the molds at the rate of about 0.2 cm/min. into 
a liquid-air bath by means of a set of reduction 
gears or by suspending the molds in a large 
Dewar flask so that the tips were about 1/3 
immersed in liquid air in the bottom of the flask. 
In the second method the process was hastened 
by lowering the molds from time to time, care 
being taken that the crystallization was always 
well above the liquid-air level. This method 
required about 2? hours to freeze the 28 cm by 
0.8 cm crystals and was employed for most of 
the crystals. 

The mercury was cleaned with dilute po- 
tassium hydroxide, dilute nitric acid, and then 
distilled once in a vacuum and twice in an air 
stream of 2.5 cm pressure by Hulett’s method." 

The mold containing the crystal was removed 
from liquid air and the ends sawed off with a 
band saw to the desired length of 25.5 cm. The 
mold was broken free from the crystal along the 
previously made grooves and the crystal fastened 
to a grooved wooden block to prevent deforma- 
tion in handling. The block containing the crystal 
was then clamped in a previously cooled wooden 
measuring bench which consisted essentially of a 
rod of seasoned hard wood containing three 
steel needle points spaced 10 cm apart and used 


"(G, A. Hulett, Phys. Rev. 33, 307 (1911). 
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as potential probes. Current leads were clamped 
on the ends of the crystal and the whole lowered 
into a Dewar flask. The distance between the 
needle points was determined in liquid air by 
immersing the wooden rod so that the needles 
protruded above the air. The cross-sectional 
area was calculated from the length, mass and 
the density measurements of Dénitz.” 

A Leeds and Northrup type K potentiometer 
was adapted to measuring the resistance by 
adding a second shunt to extend the range by a 
factor of 10. The individual and combined re- 
sistances of the two 10 cm sections of the crystal 
were determined separately, a necessary con- 
dition that the crystal be single being that the 
two sections have equal resistances. Errors from 
thermoelectric effects were eliminated by re- 
versing the current through the crystal and 
through the potentiometer and repeating the 
measurements. Readings of 10-4 volt could be 
duplicated to 10-7 volt. The resistance of a 
10-cm length of a crystal was of the order of 
10°* ohm and the measured resistance of the 
two sections always agreed with the sum of their 
separate resistances to about 0.1 percent. 

The expansion apparatus consisted of a quartz 
tube, 50 cm long and 0.8 cm inside diameter, 
suspended from a rigid support and hanging in a 
cylindrical copper vessel, 40 cm deep and 7.5 cm 
in diameter, containing isopentane. The copper 
cylinder was supported in a large Dewar flask 
48 cm deep and 10 cm in diameter. The crystal 
was placed in the quartz tube and a smaller 
quartz tube, about 26 cm long, set on the crystal 
so that the top of the smaller tube extended out 
of the larger one. On the top end of the smaller 
tube a piece of plane glass was cemented and 
served as a support for the back leg of an optical 
lever, the front legs of which rested on the same 
support used for the larger quartz tube. The 
optical lever served to measure the relative 
expansion of the crystal and an equal length of 
quartz. 

Isopentane (boiling point 20°C to 30°C) 
served admirably for a bath down to about 
— 165°C. It became quite viscous near — 165°C 
and froze a little below —167°C when stirred. 


! E. Griineisen and O. Sckell, Ann. d. Physik 19, 389 
(1934). 
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To aid in the proper stirring of the bath, a 
baffle plate was soldered across one side of the 
copper cylinder and extending from 1 cm from 
the bottom to 4 cm from the top. The bath was 
stirred by passing compressed air down a tube 
behind the baffle almost to the bottom of the 
cylinder, as the bubbles rose they kept the liquid 
circulating. 

In order to be sure that the temperature of 
the crystal was the same as that of the bath, 
the quartz tube containing the crystal was 
ground completely open in slits on opposite 
sides of the tube beginning 1 cm from the 
bottom and extending upward 22 cm so that 
the crystal was in contact with the bath. 

The temperature of the bath was measured 
with a platinum resistance thermometer cali- 
brated at the boiling point of oxygen, the 
sublimation point of carbon dioxide, the freezing 
point of mercury and the freezing point of water. 

The bath was cooled to — 165°C by pouring 
liquid air between the cylinder and the Dewar 
flask, the crystal was removed from liquid air 
and suspended in the bath and then the bath 
was allowed to warm slowly by radiation and 
conduction. The rate of temperature rise varied 
from about 10°/hr. at —160°C to 8°/hr. at 
—100°C. Readings of the expansion and tem- 
perature were made every 5 to 12 minutes. 

The length of the crystal at —160°C was 
found from cathetometer measurements between 
fiducial marks on the quartz tube in which it 
was contained, and on the tube set on the 
crystal, and from the known depth and length 
of those two tubes. 

The change in length of the crystal was plotted 
against the temperature on a large scale graph. 
The slope of this curve was measured with a 
straight-edge at 5° intervals and the expansion 
coefficient calculated from the equation: 


(ae) r= (1/Lr)(AL/AT)r. 


Since the bath could not be used below 
—165°C the mean expansion coefficient from 
— 190°C to —160°C was found. The crystal was 
suspended in a liquid air bath and the necessary 
readings taken, then the liquid air bath was 
quickly replaced with the pentane bath which 
had been previously cooled to —165°C and 
readings taken as the bath warmed. 
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A correction for the expansion of quartz was 
made from the equations of Souder and Hidnert." 
It amounted to less than one percent. 


RECRYSTALLIZATION 


In the preliminary work on polycrystalline 
rods it was found that either plastic deformation, 
due to the weight of the crystal, or crystal 
growth set in at a temperature of about — 90°C 
for, when the crystal was cooled a few degrees 
and then allowed to warm up again, the ex- 
pansion curve fell below the original curve. At 
higher temperatures this effect increased but at 
lower temperatures it was not noticeable. 

In single crystals a recrystallization or a 
breaking up of the single crystals into smaller 
crystals seemed to take place as they were 
warmed. This was probably caused by deforma- 
tion of the extremely soft crystals since Sckell" 
and Reddemann" who did not remove their 
crystals from the glass molds in which they 
were grown were not bothered by the effect. 
Three crystals which were kept in liquid air for 
three days showed no change in orientation 
during that time but Table I shows that (except 
for No. 23 which was badly bent) the three 
crystals which were warmed to the highest 
temperatures experienced the greatest changes in 


TABLE I. Changes in resistance of sections and in orientation 
during the expansion measurements. 








CRYSTAL °% Dirr. BETWEEN sin? @ Asin?@ WaARMED 
No. SECTIONS TO 
before* after*  before* after* 


23** 0.2% 1.3% 0.654 0.432 0.222 —112°C 
6 





1 9.0 991 824 = .167 —78 
20 1 10.0 077 = .238 ~——« 161 —93 . 
19 a 7.5 285 403 «118 —93 
18 J 3.2 787 =.688) = «6099 — 112 
32 2.0 2.0 540 564 .024 —112 
14 0.5 2.4 850 827 023 —I11 
26 a 0.6 975 .994 019 —112 
27 a - 907) =.920_ O13 — 113 
35 I es 150 3=.143) 007) —112 
12 0 a .046 053 007 > —I111 
21 0 6 .000 =.000 =6—.000 )=——112 
30 4 l 000 =.000)=6.000 )3=— —113 
37 0 — -000 — — —143 
36 mf — .072 — — —143 
34 1.0 — .098 — — —143 
33 0.2 — .922 — — — 143 








* Before and after the expansion measurement. 
** No. 23 was badly bent before the second measurement. 
— The dashes indicate that no measurement was made. 


18 W. Souder and P. Hidnert, Bur. Stand. Sci. Pap. No. 
524 (1926). 

4Q,. Sckell, Ann. d. Physik 6, 932 (1930). 

1H. Reddemann, Ann. d. Physik 14, 139 (1932). 
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sin? 6. Moreover, the expansion curves of these 
three crystals showed corresponding changes in 
slope and curvature at a temperature between 
—110°C and —115°C. For this reason the 
expansion measurements were carried to only 
— 115°C. 

Frequently only one section of a crystal would 
change its orientation. The change was usually 
toward a mean value of sin? 6. This behavior 
can also be noted in Fig. 2 of Griineisen and 
Sckell. 

RESULTS 

The measurements of ag are plotted in Figs. 1, 
2 and 3. Figs. 1 and 2 show ay plotted against 
sin’? @ for all the crystals at the indicated tem- 
peratures. In Fig. 3 the mean expansion coeffi- 
cients from —190°C to —160°C for the six 
crystals measured in that interval are plotted 
against sin? @. 

Table II contains the values of a, and a, 


TABLE II. Values of the coefficient of expansion parallel 
and perpendicular to the trigonal axis. 





a,X10° (ay+2a,) X10° 








i. © ay, X 10° 
—115 4.98 3.77 12.52 
—120 4.96 3.75 12.46 
—125 4.88 3.72 12.32 
— 130 4.80 3.69 12.18 
— 135 4.75 3.68 12.11 
— 140 4.66 3.60 11.86 
—145 4.64 3.56 11.76 
—150 4.60 3.55 11.70 
—155 4.55 3.54 11.63 
—160 4.49 3.52 11.53 
—190/—160 34 10.94 





4.26 3. 








read from Figs. 1, 2 and 3, and also the volume 
coefficient, a,+2a,. In Fig. 4, a, and a, are 
plotted as functions of the temperature. The 
mean values from —190°C to —160°C are not 
plotted since the other measurements are at 
definite temperatures, but the mean values indi- 
cate that the coefficient may begin to decrease 
rapidly just below —160°C. The temperature 
range is not sufficient to warrant drawing the 
curves. a, decreases more rapidly with the 
temperature than a, as Griineisen predicted 
from the elastic constants, the mean expansion 
coefficients from —78°C to —190°C, and the 
known behavior of the thermal pressure coeffi- 
cients.'® 


1% E. Griineisen and O. Sckell, Ann d. Physik 19, 407 
(1934). 
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ge es ee Fe Carpenter and Oakley.'? These authors repre- 
a 4 sented their measurements with a straight line, 
¢ 
- /N7 but when they are compared with the present 
% e 1 measurements it is evident that except for the 
ok -/ ¢ point at —89.7°C their data are better repre- 
—_ + + . = - 
| x A Fd 7 sented by the smooth curve drawn in Fig. 5. 
e + ~ . 
' i J “4 The mean value of the volume coefficient from 
mo i . —78°C to —190°C, found by extending the 
“a 8 al smooth curve to —190°C, is 12.2310-*. This 
x atl ' 4 agrees well with the value 12.210-* found by 
- [po “| Griineisen and Sckell. If the lower end of the 
x | | curve should bend down, as the mean value from 
°‘FereCeC rere ee SS a — 160°C to —190°C indicates that it does, the 
190° 140 -90° 40° 


TC 


Fic. 5. The volume coefficient as a function of the 
temperature. The circles are the present measurements, the 
+signs are the measurements of Carpenter and Oakley. 


The volume coefficients as functions of the 
temperature are plotted in Fig. 5. The circles 
represent the values found in the present 
research, the + signs are the measurements by 


agreement would not be much affected. The data 
are insufficient to give the form of the curve 
below — 160°C. 

The writer wishes to express his appreciation 
to Professor George Winchester for suggesting 
the problem, for his many helpful suggestions, 
and for his constant interest throughout the 
course of the work. 


71. 
(1931). 


G. Carpenter and F. H. Oakley, Phil. Mag. 12, 518 
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Quenching of Cadmium Resonance Radiation by Foreign Gases 


Harry C,. Lipson AND ALLAN C. G. MITCHELL, Department of Physics, New York University, University Heights 
(Received July 13, 1935) 


Measurements have been made on the quenching of the 
cadmium resonance line 3261 by He, De, CO, NHs3, Ne 
and CH,. The experiment was carried out by comparing 
the intensity of 3261 emitted from a cell containing Cd 
vapor and foreign gas with that from a cell containing Cd 
vapor only. The cells were illuminated simultaneously 
and the results obtained by the method of photographic 


photometry. The experiments were performed under con- 
ditions (low vapor pressure, short light path through the 
vapor and small foreign gas pressure) which warrant the 
application of the Stern-Volmer formula. The quenching 
cross sections obtained therefrom are: Ho, 0.67; De, 0.19; 
CO, 0.14; NH3, 0.041; Ne, 0.021; CH,, 0.012 *10-" em?. 





INTRODUCTION 


HE term quenching implies a reduction in 

the intensity of resonance radiation due to 
the presence of a foreign gas. Collisions of the 
second kind between foreign gas molecules and 
excited atoms lead to a removal of the latter 
from the excited state before they have time to 
radiate, thus causing a diminution in the in- 
tensity of the emitted radiation. Under certain 
conditions a quenching collision cross section, 
oq’, can be determined, which in turn may be 
correlated with the particular process taking 
place, if the energy levels of the foreign gas 
molecule are known. 

After Wood's first discovery! of the quenching 
of resonance radiation, a great deal of work 
was done on Hg and Na in which the cross 
sections for various gases were fairly well es- 
tablished, especially in the case of mercury. The 
first quantitative work in this field was done by 
Stuart on mercury” whose work was repeated and 
extended by Zemansky* and Bates.‘ With appa- 
ratus similar to that of Stuart, Mannkopf® 
measured the quenching of sodium resonance 
radiation. The work of Bates® and Bender’ on 
the quenching of Cd resonance radiation was 
purely of a qualitative nature and no attempt 
was made to determine collision cross sections. 
Cadmium exhibits a spectrum of the same type 
as that of mercury, the main difference, as far as 
this experiment is concerned, is that the energy 
1 R. W. Wood, Physik. Zeits. 13, 353 (1912); Phil. Mag. 
27, 1018 (1914). 

2H. Stuart, Zeits. f. Physik 32, 262 (1925). 

*M. W. Zemansky, Zeits. f. Physik 36, 919 (1930). 

4]. R. Bates, J. Am. Chem. Soc. 52, 3825 (1930); 54, 
569 (1932). 

®°>R. Mannkopf, Zeits. f. Physik 36, 315 (1926). 


® J. R. Bates, Proc. Nat. Acad. Sci. 14, 849 (1928). 
7 P. Bender, Phys. Rev. 36, 1535 (1930). 


difference between the normal state and the *P 
levels is less in the case of cadmium (see Fig. 1) 
than in the case of mercury (6°P2—6*P,=0.57 v, 
6° P,—6®P,)= 0.20 v, 6°P;—6'S p= 4.86 v). 


EXPERIMENTAL PROCEDURE 


The top view of the full experimental arrange- 
ment is represented in Fig. 2. Light from a 
quartz discharge tube, S, containing cadmium 
vapor, and through which hydrogen was circu- 
lated to carry the discharge, was focused by 
means of a quartz lens, L, on two resonance cells 
placed one above the other. In one arm of each 
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Fic. 1. Energy level diagram of cadmium. 
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Fic. 2. Diagram of apparatus. 


cell pure cadmium was distilled; the other arm 
was connected to the pumping system and to the 
source of a foreign gas. The cells could be isolated 
from the rest of the system by stopcocks and 
magnetic cut-offs to prevent cadmium from 
condensing outside the furnace. The vapor pres- 
sure of cadmium was determined by the tempera- 
ture of the side arms which was maintained at 
210°C throughout the experiment. A separate 
electric heating circuit maintained a considerably 
higher temperature in the front part of the 
furnace to prevent cadmium from condensing on 
the windows of the cells. The resonance radiation 
was taken out at right angles, passed through a 
filter to cut off the 2288 line, and 3261 was 
photographed by means of a camera fitted with a 
quartz lens. 

On the photograph there appeared two images, 
one from each resonance cell. The lower cell 
usually contained cadmium vapor plus foreign 
gas while the upper cell contained cadmium 
vapor alone and served as a standard of com- 
parison. In this manner any fluctuations in the 
light source were compensated. 

Reasonable care was taken to insure purity of 
the gases used. Hydrogen was admitted to the 
system by diffusion through a heated palladium 
tube. Nitrogen was generated by heating sodium 
azide and was stored over metallic sodium. 
Carbon monoxide was generated by the action of 
sulphuric acid on sodium formate, while methane 
was prepared by heating sodium acetate and 
sodium hydroxide. Anhydrous ammonia, ob- 
tained from a tank, was stored over potassium 
hydroxide and was further purified by condensing 
in liquid air and pumping off the residual 
noncondensible gases. Deuterium (100 percent 
pure) was obtained through the courtesy of the 
Frick Chemical Laboratory, Princeton Uni- 
versity. 
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The procedure adopted was first to reflect the 
light from the are by a mirror, ./, into a Hilger 
quartz spectograph, // (see Fig. 2) and place 
intensity-time calibration marks on a _photo- 
graphic plate; then to remove the mirror, allow 
the light to fall on the two resonance cells and 
photograph on the same plate the resonance 
radiation from the two cells, first in the absence 
of the foreign gas and then in the presence of the 
varying amounts of the gas in the lower cell. 
Each pressure change was made slowly and the 
cadmium vapor was allowed to reach equilibrium 
before a photographic exposure was made. 

The plate blackenings were traced on a micro- 
photometer and then translated into light in- 
tensities by means of the calibration marks 
mentioned above. The quenching, Q, was de- 
termined by dividing the intensity of the 
resonance radiation from the cell containing the 
foreign gas by the radiation intensity from the 
cell containing cadmium vapor only and the 
reciprocal of this quenching was plotted against 
the foreign gas pressure in mm. 


EXPERIMENTAL RESULTS 


The experimental results are given in Figs. 3, 4, 
5 and 6 and Table I. As is seen in Fig. 3, hydrogen 
quenches very efficiently, the effect being first 
noticeable at pressures as low as 0.01 mm and the 


TABLE I. Quenching slopes and collison cross sections for the 
quenching of Cd resonance radiation by foreign gases. 
































Gas H. D, CO NH; Ne CH, 
Slope (mm~') {2.45 0.50 0.15 0.063 0.022 0.016 
g 9°( X10" cm*) 0.67 0.19 0.14 0.041 0.021 0.012 
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Fic. 3. Quenching of Cd resonance radiation (3261) by 
H> and Ds. 
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Fic. 4. Quenching of Cd resonance radiation (3261) by CO. 
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Fic. 5. Quenching of Cd resonance radiation (3261) by No, CH, and NHs 


intensity of the resonance radiation being re- 
duced to half of its original value at about 0.4 
mm of He. The quenching by deuterium is of the 
same order of magnitude as that of hydrogen 
although the “‘half-value” occurs at 1.5 mm. 
Carbon monoxide, Fig. 4, is less efficient than 
hydrogen, 0.4 mm decreasing the resonance 


radiation intensity to only about 95 percent of 
its maximum. Ammonia, Fig. 5, quenches moder- 
ately, the curve lying between those of nitrogen 
and carbon monoxide. Nitrogen and methane 
are very ineffective as quenching agents, the 
phenomenon not being observed until a pressure 
of a few millimeters was reached and the intensity 
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F1G. 6. Quenching of Cd resonance radiation (3261) by foreign gases. 


dropping to approximately half-value at 45 mm 
in the case of Ne and 60 mm in the case of CH. 
The relative quenching efficiency of all the gases 
is shown in form of slopes in Fig. 6. 


DISCUSSION AND INTERPRETATION OF RESULTS 


The correct interpretation of quenching ex- 
periments is affected mainly by two factors: 

1. The diffusion of imprisoned resonance radiation, i.e., 
reabsorption of the emitted resonance radiation on its way 
out through the gas to the exit window, giving rise to 
secondary radiation, which absorbed again will produce 
tertiary radiation, and so on. 

2. The alteration of the width and the position of the 
absorption line relative to the exciting line, which is due 
to the collisions of atoms with foreign gas molecules 
(Lorentz broadening). 


Unless the above difficulties are overcome, a 
quenching curve can give no information leading 
to the computation of collision cross sections and 
in some cases*:® it is not even evidence that 
quenching takes place at all. 

The present work was performed under con- 
ditions in which these two effects were largely 
eliminated. Thus, Lorentz broadening was re- 
duced by operating with pressures of foreign gas 
as low as was compatible with the quenching 
process. As to diffusion of the imprisoned reso- 
nance radiation this was eliminated, first by 
keeping the cadmium vapor pressure low (5.6 
10-4 mm corresponding to 210°C at which the 


absorption coefficient at the center of the 3261 
line is about 0.4) and second, by placing both the 
entrance and the exit window of the resonance 
cell near a common edge, so that only a small 
portion of the Cd vapor was exposed to the 
incident radiation, thus permitting only the 
primary resonance radiation to be re-emitted. 

The fulfilment of the above conditions warrants 
an application of the Stern-Volmer formula:* 


Q=1/(1+7Z9) (1) 


where Q= quenching of the resonance radiation, 
t= the lifetime of the excited atom, and Zg= the 
number of quenching collisions per sec., per cc, 
per excited atom. 

Since Zg varies linearly with the foreign gas 
pressure, p, a plot of the reciprocal of Q against p 
should result in a straight line. This was found to 
be true in the present investigation, as attested 
by the above graphs. Above a certain pressure, 
the deviation from straight lines can probably be 
accounted for by Lorentz broadening.’ 

The quenching collision cross sections, ¢@’, 
were determined as follows: 

The Stern-Volmer formula, (1), gives a relation 
between the quenching Q and the total number of 
quenching collisions Zg. The latter can be 
calculated on the basis of the Maxwellian distri- 


*O. Stern and M. Volmer, Physik. Zeits. 20, 183 (1919). 
® Mitchell and Zemansky, Resonance Radiation and 
Excited Atoms (The Macmillan Co., 1934), p. 224. 





61 
he 
ce 
all 
he 
he 





QUENCHING OF RESONANCE RADIATION 629 


bution of velocities, resulting in the equation 
Zq=2Nnog*2rkT(1/Mi+1/M2)}' (2) 


where J; and 1/.=molecular weights of the 
colliding particles, V and n=their molecular 
concentrations, and k= universal gas constant or 
expressed in terms of the pressure in mm of Hg 


Zo= 2666.60 o?[ 27N/kT(1/Mi+1/M2)]}-p. (3) 


From (1) and (3) it follows that the collision 
cross section is given by 
1/(Q-p) 
¢Q°= ——. (4) 
72666.6[27N/RT(1/.1,+1/ Me) }} 

The values 1/(Q-p) are the slopes of the 
quenching curves and could be read for each gas 
from Figs. 3, 4 and 5; they are given in Table I. 
The temperature of the cadmium vapor was 
210°C and the mean life for the 5°P, state of Cd 
was taken as r= 2.5 X 10° sec.'® The values of the 
quenching collision cross sections thus computed 
from (4) are given in the same table. 

The high quenching efficiency of hydrogen was 
interpreted by Bender’ as due to a collision of an 
excited Cd atom in the 5°P, state with a normal 
He molecule, resulting in a formation of an 
unexcited cadmium hydride molecule and an 
atom of hydrogen 


Cd (5°P;)+ HsoCd H+H 
3.78v 444v 0.67v _ 


(5) 


This process seems to be highly probable, as 
indicated by the energy balance of the equation, 
in which the energy required for the above 
reaction differs from the energy available by less 
than half of one percent. This was further 
substantiated in Bender’s experiment by the 
appearance of CdH bands in the resonance cell 
and the formation of atomic hydrogen found as 
water frozen out by the liquid air trap, due, 
probably, to a reduction by atomic hydrogen of 
some oxide present in the Cd supply. 

However, simple calculations made on the 
basis of the data given by Jevons" indicate that 


1” Reference 9, p. 147. 
"\W. Jevons, Report on Band-Spectra of Diatomic Mole- 
cules (The Physical Society, London, 1932), p. 269. 


Hz has a vibrational level of the normal electronic 
state at 3.73 volts, which, from purely energetic 
considerations makes the following process also 
possible: upon collision, the excited Cd atom 
gives up 3.78 volts to raise the Hz molecule to a 
3.73-volt vibrational level: 


Cd (5*P1) +H2Cd (5'So) + (He2)- 
. (6) 


3.78 v 3.73 v 


Although in this reaction the energy discrep- 
ancy is larger (about 1.5 percent), still the above 
process seems to be very plausible. 

Similar considerations can probably apply to 
deuterium (energy of dissociation of Ds» is the 
same as that of He; there is also a vibrational 
level at 3.72 volts) although no attempt was 
made to detect the existence of the CdD bands. 

The quenching by CO, Ne and CH, cannot be 
accounted for by a dissociation process, since the 
energies of dissociation of these molecules are 
far in excess of the energy available (10.0 volts 
for CO, 9.1 for Ne, 16.3 for CH,4). However, the 
carbon monoxide molecule has a vibrational level 
at 3.83 volts, which being near the excitational 
level of the 3261 Cd line makes the following 
reaction possible: 


Cd (5°P,;)+CO3Cd (51S) +(CO), 
: i 
3.78 v 3.83 Vv 


In this case the energy discrepancy is 1.5 percent. 
The nearest vibrational level of nitrogen is at 
3.48 volts, indicating an energy discrepancy of 8 
percent; concerning the CH, molecule there are 
not sufficient band spectra data known to con- 
sider the quenching by methane from above 
point of view. However, at the temperature used 
(210°C), the average kinetic energy of translation 
of a molecule (3k7) is about 0.07 volt, which 
corresponds to the energy difference between the 
excited 5°P; and the metastable 5*P, state of Cd. 
Thus, perhaps, the quenching by Ne and CH, can 
be explained in terms of a transfer of the excited 
Cd atom to the metastable state, resulting in a 
loss of 0.07 volt to the molecules in a form of a 
translational kinetic energy, e.g., 

Cd (5°P:)+N20Cd (5°Po)+(Ne) K.E. 
. (8) 

3.78 Vv 3.71 Vv 0.07 v 
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Of course, the reverse process is plausible and 
undoubtedly takes place also, in which a collision 
between a metastable atom and the gas molecule 
results in a transfer of the Cd atom back to the 
58P, state, but as has been shown by Klein and 
Rosseland® this type of a collision (collision of 
the first kind) is less probable than the former, 
since a fast molecule will remain in the neighbor- 
hood of an atom during a shorter interval of time 
than a slow one. This point of view seems to be 
confirmed by Bender's observation of an increase 
in the relative intensity of the 3404 line over the 
remainder of the spectrum in the presence of Ne 
and CO but not in the presence of He. This line 
results (see Fig. 1) in the absorption of an energy 
quantum from the 5*P, state to the 5°D, state 
and a subsequent re-emission, and thus a heavy 
population of the metastable state will favor an 
increase in its intensity. 

Because of the meager data existing about the 
energy levels of the ammonia molecule it is not 
possible to ascribe any definite process to the 
quenching by ammonia. No evidence was ob- 


tained for its decomposition under the influence’ 


of excited cadmium atoms, in agreement with 
Bates and Taylor. The test is not conclusive, 
however, since the intensity of the source may 
have been too small to produce a measurable 
chemical effect, although it was large enough for 
the measurement of quenching. In view of the 


12Q,. Klein and S. Rosseland, Zeits. f. Physik 4, 46 


(1921). 
18 J. R. Bates and H. S. Taylor, J. Am. Chem. Soc. 50, 
771 (1928). 
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fact that some decomposition might have oc- 
curred and the resulting products might have had 
an effect on the intensity of the resonance 
radiation, exposures were made at the same 
ammonia pressure, first with a fresh sample of gas 
and next with a sample which has been irradiated 
previously for 15 minutes. There was no de- 
tectable difference in the quenching between the 
two Cases. 

Lastly, the collision between any one of the 
above molecules and the Cd atom cannot result in 
raising them either to an excited electronic level 
or to some rotational level, in the first case, the 
energies required being too high (the first exci- 
tational level of He is 11.1 v, CO—4.8 v, Ne 
—8.5 v) in the second case the energies being 
only of the order of a fraction of a volt. 

In conclusion, it ought to be emphasized that 
it is difficult to interpret quenching in terms of 
one process alone, the reactions taking place 
within the resonance tube being of a very 
complicated nature. It may be that quenching by 
a particular gas can be explained by all the 
processes mentioned above, there being only a 
degree of effectiveness of one reaction over the 
other. Thus, perhaps one can venture an opinion 
that the process, in which the discrepancy be- 
tween the energy required and the energy 
available is small, is more responsible for the 
quenching phenomenon than the one in which a 
corresponding discrepancy is large. 

The authors wish to acknowledge their in- 
debtedness to Dr. R. L. Garman of the Chemistry 
Department of the Washington Square College 
for making the microphotometer traces. 
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Intensity of Diffraction of Electrons by ZnO 


H. J. YEARIAN, Purdue University 
(Received July 18, 1935) 


Electron diffraction patterns of ZnO powder have been 
obtained in the range from 10 to 80 kv by using two types 
of camera. The intensity distribution has been determined 
photographically and shows marked deviations from the 
intensity distribution of the corresponding x-ray pattern. 
It has been shown that these deviations are not due to 
impurities, difference in structure due to methods of prepa- 
ration, coincidences with another form of ZnO, or due to 
preferred orientation. The scattering as a function of 
sin (@\) is discussed on the basis of additional surface 


reflection, dynamic reflection and refraction of the elec- 
trons, and it is found that also these effects cannot ac- 
count for the anomalous intensity distribution. Complete 
agreement between experimental and calculated curves 
for the scattering factor is obtained by assuming a dis- 
tortion of the electron cloud (M shell of the zinc and L 
shell of oxygen) resulting in a different parameter for the 
distance Zn nucleus—oxygen nucleus, and Zn electrons- 
oxygen electrons. The resulting structure factor has the 
form 


S=Z(dn)— VF q, (Zn)e2* in, +-Z(O)e2*ial— Fy, (O)e2*i(atbn, od! 


[' is known that while the scattering of x-rays 
is not influenced by the nucleus of the 
scattering atom, electron diffraction is affected 
by the nucleus as well as by the electronic cloud ; 
electron diffraction is therefore able, at least in 
principle, to locate the position of both nucleus 
and electrons in a crystal lattice. 

In view of these considerations experimental 
investigations of the intensity of diffraction of 
high speed electrons in a powder, using photo- 
graphic recording, have been undertaken, with 
the object of determining the atom factor for the 
diffraction of electrons.' 

The wave-mechanical theory of the electron 
diffraction atom factor FE, has been put in an 
especially simple form by Mott,’ 


E=(Z—F), (sin @,d)?. (1) 


@ is the half-angle of scattering, Z the atomic 
number of the scatterer, and F its x-ray atom 
factor. This expression is derived for a free atom 
and is based on the assumption that the atom is 
spherically symmetric about its nucleus, and 
that the plane wave of the electron beam is 
undistorted by the atom. The first of these 
assumptions, as Mott points out, may not be 
valid for an atom bound in a crystal lattice, but 


'H. J. Yearian and K. Lark-Horovitz, Phys. Rev. 42, 
905 (1932); K. Lark-Horovitz, H. J. Yearian and E. M. 
Purcell, Phys. Rev. 45, 123 (1934); K. Lark-Horovitz, 
I. J. Yearian and J. D. Howe, Phys. Rev. 47, 331 (1935). 

*>N. F. Mott, Proc. Roy. Soc. Al27, 658 (1930); G. P. 
Thomson, Proc. Roy. Soc. A125, 352 (1929); H. Mark and 
R. Wierl, Zeits. f. Physik 60, 741 (1930). 


the second should be acceptable for electrons of 
sufficiently high velocity. 

This expression has been found to agree quite 
well with experiments carried out on gold, silver 
and aluminum.? 


EXPERIMENTAL PROCEDURE 


Diffraction photographs were obtained in two 
cameras illustrated in Figs. 1a and b. 


The first camera, Fig. la, consists of a brazed and sol- 
dered brass box B 25 cm long closed with a lid and rubber 
gasket. Electrons from a very small conical filament of 
tungsten, F, are accelerated to the anode by a potential 
of the order of 20-30 kv, and a fine beam is formed by 
the slits S; and .S, drilled with 0.15 mm and 0.1 mm holes, 
respectively. Ss is a shielding slit on which the specimen 
is mounted. The diffraction pattern is recorded on a plate 
(Eastman lantern slide) at P. A two-stage Hg diffusion 
pump backed by a Cenco Megavac produced a vacuum of 
10-* mm Hg. The potential for this camera is supplied by 
a 30 kv full wave filtered rectifier. 

In order to obtain the multiple exposures necessary for 
establishment of the intensity-density characteristic of 
the plate, a large diaphragm / with two 30° segments 
cut away is placed in front of the plate. This diaphragm 
is rotated about an axis parallel to the beam by a weight 
(not shown) and stopped in five adjacent positions by 
a catch operated by a solenoid. The diaphragm is con- 
trolled automatically by a motor driven switch which 
also automatically times the five exposures in the ratio of 
1:2:3:4:5 by means of a pair of small Helmholtz 
coils. All the necessary exposures are easily made in 

‘The apparent discrepancy in the case of Al disappears 
when a temperature correction is applied. The comparison 
with the analysis of a number of other substances will be 
presented in later papers. 
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approximately one second more than the total exposure 
time (of the order of 20 seconds), so that current-voltage 
conditions are relatively easy to maintain. 

The second camera, Fig. 1b, has been used up to 80 kv. 
It is assembled of brass parts with the removable joints 
sealed with sulphur free rubber gaskets. Electrons from a 
small filament F, mounted adjustably on a flexible bellows, 
B, form a narrow beam defined by the 0.15 mm hole in the 
anode S;, and the 0.10 mm hole in slit S.. The specimens 
are mounted either on the shielding slit S; or on a universal 
mounting operated by the double cone joints 7. The 
camera is mounted directly Holweck molecular 
pump backed by a Cenco Megavac, giving a very high 


on a 


pumping speed. The potential is supplied by a 120 kv 
full wave filtered rectifier 

The pattern is recorded on plates held in the magazine P. 
This magazine holds seven plates each of which may be 
pulled over the camera opening, exposed and returned, 
by means of a traveling carriage operated by a ground 
joint J. Also a screen with a rectangular opening may be 
brought into position over the camera so that five expo- 
sures may be made in strips on each plate in order to 
obtain its intensity-density characteristics. The exposures 
are timed by contacts on a motor-driven screw; these are 
adjustable to 20 sec. in steps of 1,20 sec.* 


Specimens of ZnO were formed by catching 
the fumes from a zine are on the last diaphragm 
of the slit system. The deposited crystallites grow 


*This camera has recently been improved to a quite 
universal form. A detailed description will be published 
shortly. 


into the slit opening, and with a sufficiently fine 
deposit excellent diffraction patterns are ob- 
tained, Fig. 2. 

The measurements of the rings in this pattern 
agree very closely with those calculated from 
x-ray data, as is shown in columns 1 and 2 of 
Table I. Columns 3 and 4 give a comparison of 
estimated intensities. It is evident that there are 
large discrepancies in the intensities which must 





Fic. 2, Pattern of ZnO: 71.5 kv, 1.0 ma, 0.75 sec. exposure. 
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TABLE I. 








ZnO Hex. Close packing a@=3.24A c/a=1,599 











Est. Est. Mean Si Est. Est. Mean Si 

sin 6/X int. int. rel. Ave. Obs. | sin 6/r int. int. rel. Ave. Obs. 
hkl obs. calc. elec. x-ray int. dev. f S? mean hkl obs. cale. elec. x-ray int. dev. f So mean 
001 .0965 0 301 .5432 0 
100 .1783 .1782 4 4.5 1.60 1! 6 1 2.14 213 .5580 .5533 4 4.75 .524 .105 24 3 10.6 
002 .1941 .1930 6 3.25 2.09 17 2 4 2.64 302 .5682 3 3.75 .273 .023 12 4 10.2 
101 .2042 .2026 5 5 3.46 -27 12 3 1.87 115 .5727 0 
102 .2637 2627 5 4 1.46 .10 12 3 4.64 | 006 .5791 .5789 ? 1 041 O11 2 4 10.0 
003 .2895 0 205 .6000 .5999 3 3.75 .277 024 12 3 12.6 
110 .3093 .3086 4 5 1.40 15 6 4 5.21 106 .6078 .6059 2 O55 .004 12 1 10.8 
111 .3234 0 303 .6078 0 
103 .3407 3394 6 4.75 2.97 32 12 3 8.19 214 .6082 .6093 1 .048 .002 24 1 7.39 
200 .3594 .3563 2 2 243 014 6 1 6.74 220 6173 1 3.5 .080 .024 6 4 10.0 
112 .3654 .3640 4 5 1.48 .15 12 4 6.17 221 .6248 0 
201 .3718 .3692 3 3.5 528 .065 12 3 4.46 310 6424 .066 014 12 1 14.4 
004 .3872 .3860 1 1.75 236 .018 2 4 7.10 222 .6505 .6468 .155 .026 12 4 11.4 
202 4074 4053 2 2.25 .266 054 12 1 7.26 311 6497 131 021 24 3 8.63 
113 .4232 0 312 6586? .6706 .062 .004 24 1 11.2 
104 4269 -4252 3 1,25 197 .040 12 1 7.08 215 .6776 .6743 .264 .033 24 3 13.5 
203 4614 -4592 + 4 571 .079 3 8.91 007 6754 0 
120 4752 4715 1 2.9 226 .053 12 1 10.5 206 6798 .026 003 12 1 10.6 
211 4835 4812 2 3.5 370 .066 2 3 5.89 223 .6741 
005 4826 0 107 .7028 .6986 
114 4973 4943 2 2.5 .320 054 12 4 7.13 313 .7099 .7047 
212 .5165 .5094 5 2.5 .216 .115 24 1 9.26 224 .7280 
105 .5144 4 .574 -115 12 3 12.45 216 .7466 .7467 
204 .5297 .5254 1.25 065 010 12 1 7.78 
300 5386 .534060 2 3 271 .036 6 4 11.9 








| 
| 
| 
} 
| 
| 
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be accounted for. The discrepancy is also ap- 
parent in the comparison of electron and x-ray 
microphotometer records in Fig. 3. 

Inspection shows that planes for which the 
third index is small are relatively too weak in the 
electron diffraction pattern. For example, in the 
x-ray pattern 110 is stronger than 103, but this 
relation is reversed in the electron pattern. It is 
interesting to note that these same differences 
are apparent in the electron diffraction patterns 


of ZnO reproduced by Ponte‘ and Dauvillier® 
and by Finch and Quarrell.* According to formula 
(1) electron intensities should differ from the 
X-ray intensities, but only as a monotonous 
function of sin (@/A) and not irregularly from 
peak to peak. 


H. Ponte, Comptes rendus 188, 244 (1929), 
5M. A. Dauvillier, J. de physique 3, 229 (1932). 
®°G. I. Finch and A. G. Quarrell, Proc. Roy. Soc. A141, 
398 (1933). 








Fic. 3. Microphotometer curves of ZnO diffraction patterns. A, x-ray pattern; B, electron 
pattern; C, ‘‘most probable" trace; D, trace C converted to optical density; £, curve D 
reduced to flat background and partially resolved. 
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DISCUSSION 


Intensity anomalies of this kind might be 
caused by impurities, preferred orientation, 
method of depositing powder and influence of its 
support, or by an erroneous application of the 
kinematical theory of diffraction neglecting ab- 
sorption, extinction, refraction and the possible 
influence of cross gratings. In the search for the 
possible causes of the unexpected intensity 
distribution, over 200 plates have been taken 
under different conditions at potentials ranging 
from 10 to 80 kv. 

No difference has been found between the 
diffraction patterns from C.P. and_ spectro- 
scopically puref zinc. Specimens have been pre- 
pared on brass, copper, gold, platinum and 
molybdenum slits without change. To check the 
effect of the form of powder deposit, pictures 
have been taken with both commercial ZnQO, 
and oxide from the arc collected on glass and 
pressed on the slit. There was no difference in 
the diffraction patterns from deposits of very 
large or very small particle size. Zinc oxide from 
a zinc arc sometimes forms long threads of 
crystallites. Patterns have been obtained from 
such single threads, from powder deposited on 
silk fibers and on thin films of celluloid. All 
these patterns are identical with the one illus- 
trated. 

There is no indication of preferred orientation, 
but to test for an axis of orientation parallel to 
the electron beam, a deposit was made on a slit 
which could be rotated. No change was produced 
when the slit was rotated in the range from 90° 
to 30° with respect to the beam. X-ray photo- 
graphs of one of the specimens, moreover, gave 
the regular x-ray intensities. 

All x-ray data on the structure of ZnO are in 
agreement with a single space group C,,‘*. A 
different surface structure, undetectable by 
x-rays, might affect the electron intensities. A 
comparison, however, with the intensities of the 
other possible space groups’ shows that no 
combination of these can account for the ob- 
served electron intensities. 

Coincidence of lines of some contamination 


+ The pure zinc was obtained from Adam Hilger, London, 
England. 
*Wryeoff, The Structure of Crystals, first edition, p. 275. 
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could falsify the intensites. An actual impurity 
is improbable since no extra lines were observed. 
A comparison of the ZnO pattern with the cubic 
pattern reported by Bragg and Darbyshire,* the 
pseudomorphic form of the oxide reported by 
Finch and Quarrell,® and with zinc, shows several 
coincidences but in all cases there are equally 
strong or stronger lines in these patterns which 
do not coincide and would be observable if 
these substances were present. 

From these considerations it would 
impossible that the anomalous intensities ob- 
served can be due to contamination or to 
preferred orientation. Further it appears very 
improbable that they can be explained by 
absorption or extinction effects. In primary 
extinction one might expect that intense re- 
flections should be weakened more than weaker 
ones, but there is no apparent connection 
between the deviations in the intensities and the 
expected strength of reflection. 

As will be shown in more detail later, the 


seem 


nearer a plane is to being perpendicular to the 
c axis, the more its intensity deviates from 
normal. To account for this by extinction would 
require a corresponding regularly increasing 
degree of imperfection of the crystal for such 
planes. This seems very artificial. The dynamical 
theory of electron diffraction developed by 
Bethe® in considering the effects of extinction 
and interaction of diffracted beam with the 
incident beam, predicts the appearance of re- 
flections not possible in the usual kinematical 
theories. Faint “inner rings,’’ discussed below, 
are sometimes found in the ZnO pattern which 
might be attributed to such dynamical effects. 
Their intensities, however, are less than one 
percent of the stronger regular reflections, and 
would indicate that the deviations of fifty 
percent or more observed in the regular pattern 
are not explainable by dynamical considerations 
of extinction, etc. 

One could possibly explain the dependence of 
the intensity deviation on the direction of the 
diffracted rays in the crystal, as being due to 
absorption. The absorption in crystallites much 
longer parallel to the c axis than perpendicular 


’W. A. Bragg and J. A. Darbyshire, Trans. Faraday 
Soc. 28, 522 (1932). 
*°H. Bethe, Ann. d. Physik 87, 55 (1928). 








DIFFRACTION OF 


to it would weaken reflections from planes 
parallel to the axis more than it would reflections 
from planes perpendicular to the axis. Kirchner'? 
has found that the total absorption coefficient 
(scattering plus true absorption) varies, in mica, 
approximately inversely as the electron velocity. 
Such a rapid variation of absorption should 
produce considerably different intensity devia- 
tions at different velocities. No noticeable differ- 
ence in the intensity patterns has, however, been 
obtained over the whole range from 7 kv to 
80 kv, and one must conclude that absorption 
does not account for the anomalities observed. 
There is still the possibility that refraction 
produces the effect. On a few plates one or more 
very faint but sharp rings smaller than the first 
permissible reflection (‘‘fractional orders’) have 
been noted; these usually occurred when the 
specimen contained large crystallites. These have 
been explained by K. Lark-Horovitz' as a 
refraction effect since he has shown that although 
the refraction is negligible when high velocity 
electrons are incident normally on a diffracting 
layer and are diffracted from an inner plane, 
when the diffraction occurs from an exposed 
surface, refraction is not negligible. Indeed a 
lattice potential of the order of only 4 volts is 
necessary to produce a ‘“‘half order’ diffraction 
from 002 of ZnQO.* This effect, however, can 
hardly be the cause of the anomalous intensity 
distribution since calculations have shown that 
the refraction is appreciable only when the 
incident or diffracted beam, or both, makes a 
very small angle with the crystal surface, and 
can practically enter, therefore, only in the case 
of planes which are developed faces of the 
crystallites. This necessity of very small angles 
also requires that the surface be exceedingly 
flat, ‘electron optically smooth,’’ and of con- 
siderable extent, of the order of 10°A. From these 
considerations it is obvious that the whole pattern 
could hardly be affected by refraction in the way 


10 F, Kirchner, Ann. d. Physik 13, 45 (1932). 

"KK, Lark-Horovitz and H. J. Yearian, Phys. Rev. 43, 
276 (1933). 

* This explanation of “fractional orders’’ does not 
require that they be rational fractions of some plane; 
they can be any fraction depending upon the lattice po- 
tential and the plane spacing. Recently also G. P. Thomson 
has advanced this explanation (G. P. Thomson, Phil. Mag. 
18, 640 (1934)). 
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observed. For example, the 100 and 110 in- 
tensities are reduced to approximately half as 
compared with reflections from 001. This would 
require that at least half the crystallites have 
these planes developed as large electron optically 
smooth faces and at the same time have no 001 
faces so developed. But in a pseudolayer-lattice 
like ZnO one would expect the latter type of 
development to be at least as frequent as the 
former, and probably much more frequent.** 
Furthermore, the products of the many necessary 
refractions should appear either as_ visible 
“inner” rings, or as a broadening of the ab- 
normally weak rings on the inner side;f this 
does not occur. Any broadening observed affects 
all rings alike without influencing their relative 
intensities. 

The intensities would also be influenced by the 
diffraction from two-dimensional lattices super- 
imposed upon that from the space grating. The 
occurrence of cross lattice diffraction has been 
frequently observed, ' and has been theoretically 
treated by v. Laue.'* Transmission through two- 
dimensional networks distributed at random 
results in a ring pattern. This has the appearance 
of the usual space lattice pattern but because of 
the variable angles of incidence on the various 
crystals, these rings, although sharp on their 
inner boundaries, fall off less rapidly on their 
outer edges. In this type of diffraction from the 
hexagonal basal planes of ZnO, the rings formed 
coincide with the space grating reflections of the 
type hkO, which were shown above to be too 
weak. Although cross grating reflections from 
prism planes coincide with some of the observed 
reflections which are too strong, it is almost 
certain that the basal plane effect should pre- 
dominate, so that the observed intensities cannot 
be explained by cross grating effects. Further- 
more, the asymmetry of the lines associated with 
this type of diffraction is not observed in the 
pattern. 





**See for instance the development of zinc crystals, 
Handbuch der Physik (1934), Vol. 24, No. 2, pp. 811-821. 

t An effect similar to the one observed by L. H. Germer 
in another connection. Phys. Rev. 44, 1012 (1933). 

2 Kikuchi, Jap. J. Phys. 5, 83 (1928); W. L. Bragg, and 
F. Kirchner, Nature 127, 138 (1931); Raether, Zeits. f. 
Physik 78, 527 (1932); Steinheil, Zeits. f. Physik 89, 50 
(1934). 

Sy, Laue, Zeits. f. Krist. 82, 127 (1932). 
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ATOM FACTOR AND DEVIATION FROM 
SPHERICAL SYMMETRY 


Since none of the possibilities discussed above 
seemed to explain the observed intensity distri- 
bution, an atom factor calculation has been 
made in order to investigate the atomic scattering 
in the crystal lattice. 

The intensity calculations were made in the 
following way. Each of the five exposures on a 
plate was microphotometered on a recording 
instrument and the optical density of several 
peaks was computed and plotted against expo- 
sure. The resulting density-intensity curves were 
found to be linear below densities of 0.9 to 1.0. 
An exposure which nowhere exceeded this limit 
was then photometered two or three times at 
different positions. The superposition of these 
traces gave a ‘‘most probable”’ trace practically 
free of uncertainties due to grain of the plate. 
This trace was then converted to optical density 
by means of an improved design of the “loga- 
rithmic protractor” described previously.'* 

After drawing in a background curve, the 
converted trace was reduced to a flat background 
and the positions of all peaks marked. The 
overlapping peaks were then resolved under the 
condition of consistent summation of ordinates 
and areas, and their areas taken as a measure of 
the corresponding intensities. The various steps 
are indicated in Fig. 3. 

ZnO is a hexagonal close packed lattice Zn 
(000, 1/3 2/3 1/2) and O (009, 1/3 2/3 1/2+ ); 
p= 0.375 and c/a= 1.599." This gives a structure 
factor which may be written 


S= [Ezt+Ece?*'?' J[1 fe? til Bh+2/3k+1/20) | 
= S,’ x So. 


S,’ is the structure factor of the molecule and S, 
takes into account the two molecules in the 
basis. By (1), (on the assumption of spherical 
atom symmetry, etc.) 


Z-F 


sin? 6/d? 


4H, J. Yearian, Rev. Sci. Inst. 4, 407 (1933). 
1 Structureberichte, 1913-1918, P. P. Ewald and C. 
Herman, page 79. 
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so that we may define 
Si = S;' Xsin® 0/xN= (Zz — Fz) +(Zo— Fo)e?*™”'. 
In terms of the observed intensities 


f\S/? 





=, , 
sin? 0 

where f is the frequency factor, and 5S; can be 
calculated 


S1= C sin* O\/I/f Se) - 





since S: is a purely geometrical term. 

The mean intensities of fourteen of the best 
ZnO plates in the range 18-20 kv are shown in 
column 6 of Table I.* The mean deviations 
given indicate a quite high reliability, particu- 
larly for the stronger peaks. The experimental 
values of S, deduced from these mean values are 
shown in column 10 and are plotted against 
sin (@/A) in Fig. 4. In Fig. 5 theoretical F factors 
for Zn and O are shown and the theoretical 
values of the scattering functions S; calculated 
from them. Because of the falling-off of the F 
curves with angle, each Z—F term in S; rises 
with angle and the exponential term produces a 
separate curve for each value of index /. It is at 
once obvious how greatly the experimental 
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F1G. 4. Observed scattering curves. 


* Calculations have also been carried out for electron 


velocities up to 80 kv and for temperatures as low as —183 


without appreciably different results. 
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Fic. 5. F curves for zinc and oxygen. FT, theoretical 
(Pauling and Sherman); FX, experimental x-ray; FE, 
experimental electron. S,, ideal scattering curves derived 
from FT. 


intensities differ from the theoretical. Most of 
the experimental scattering curves at small 
angles are much too high; the displacement 
above the expected position is in general greater 
the smaller the angle and the larger the value of /. 

This behavior gives a hint of possible explana- 
tion, for if g is the angle between a reflecting 
plane and the c axis, then 


sin g= (A/sin 6) (1/2c) 


so that the larger / or the smaller sin (@/) the 
more nearly is the plane perpendicular to the c 
axis, and according to the experimental curves, 
the higher is its reflecting power above normal. 
This could be accounted for if the x-ray F curves 
should be decreased below their normal values 
at the larger values of y. This would correspond 
to atomic electron clouds drawn out from a 
spherical symmetry into an ellipsoidal form with 
the major axis in the c direction. 

X-ray Debye Sherrer photographs of ZnO 
analyzed in the same way as the electron diffrac- 
tion photographs, but using the appropriate x-ray 
intensity formula, gave scattering curves too 
consistent among themselves to justify the large 
deviations in F necessary to fit the observed 
electron curves under the assumption mentioned 
above. From these x-ray data the F curves 
(FX) for Zn and O were found. (The parameter 
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p was evaluated to be p=0.378+0.003.)+ The 
greater slope of these curves, as compared with 
the theoretical, (FT), may be caused by anoma- 
lous dispersion, since Cu K radiation was used. 
An absorption correction after the method of 
Rusterholz'® or Blake’? would result in still 
steeper curves. These differences in slope are 
quite consistent with a temperature correction 
of the form e~8 si"? 4%. Tf the Debye-Waller 
exponent is used as a rough approximation for 
the correction in a composite hexagonal crystal, 
the data are satisfied by a characteristic temper- 
ature of = 355°. 

Another explanation of the rise of the observed 
electron scattering curves with increasing angle 
¢ between c axis and reflecting plane, is suggested 
by the fact that the electron and nuclear scatter- 
ing is out of phase. If the center of the electron 
clouds is not quite coincident with the nucleus 
but is displaced slightly along the ¢ axis, then 
an additional phase shift will enter in the 
structure factor and this will increase the re- 
sultant amplitude. This increase would be 
greatest for large values of / and for F values 
comparable with Z, i.e., at small values of 
sin (@/A). This possibility seems reasonable also 
from considerations of the ZnO structure. In an 
ideal hexagonal close packed lattice with axial 
ratio c/a= 1.633 and parameter p= 3/8= 0.375, 
each atom is at the center of gravity of a tetra- 
hedron formed by atoms of the opposite kind. 
In the actual ZnO lattice these relations are 
disturbed. The axial ratio is c/a=1.599 and 
p=0.375. With these values the Zn —O distance 
along c and the altitudes of the tetrahedrons in 
this direction are decreased about two percent. 
The distance from a Zn atom to O atom 
along c is 1.94A, but to its three other O neigh- 
bors the distance is 2.041A, a difference of 5 
percent. The angles between Zn—O bonds are 
also considerably affected. To maintain equi- 
librium with these changes, a slight shift of the 
electron center away from the nucleus might 
very well arise. 

A first approximation to the effect of such a 
shift was obtained in the following way. Since 


+ Precision measurements are now being carried out by 
K. Lark-Horovitz and collaborators. 

16 Rusterholz, Helv. Phys. Acta. 4, 68 (1931). 

17 F.C. Blake, Rev. Mod. Phys. 5, 169 (1933). 
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Fic. 6. Calculated scattering curves. Model of ZnO 
molecule with electron shells displaced to produce dotted 
curves; (scattering effect of zinc N shell negligible). 


the scattering curves for /=0 and /=1 would be 
least affected by any dissymmetry along the c 
axis, these two experimental curves of Fig. 4 
were used to calculate the experimental F curves 
(FE) shown in Fig. 5. Assuming the parameter p 
for the distance between zinc and oxygen electron 
centers, g for that between zinc and oxygen 
nuclei, and a for that between the zinc nucleus 
and electron center, the structure factor 5S; 


becomes 
Si=Zan- Fz,e77'*! + Zoe?* 4! — Foe? tier at. 


With the experimental F values and assumed 
values of p, g and a, this complex summation 
was carried out. All combinations within the 
range 0.37<p<0.39; 0.36<q¢<0.43; —0.05<a 
<+0.05 were tried. The only combination which 
gave curves approaching the experimental ones 
was with p= 0.38, g=0.37, and a= —0.015. 

In these trials it was soon found that a very 
satisfactory set of curves could be constructed if 
it were allowed to let parameter a increase with 
angle yg. This would correspond to saying that 
the effective electron diffracting center of the 
atom, as approximated by a spherical distribu- 
tion, moves farther from the nucleus for planes 
approaching perpendicularity with the c axis. 
This would seem to be at least possible for an 
electron distribution somewhat pear shaped in 


YEARIAN 


form. The best set of curves obtained under this 
assumption, with p=q=0.38, a=+0.01, for 
¢ <55° and a= + 0.02 for ¢ greater than 55°, is 
shown as full lines, in Fig. 6. Comparison with 
the observed curves, Fig. 4, shows a rather 
remarkable agreement. 

The satisfactory nature of this agreement was 
thought to warrant an attempt to justify the 
rather arbitrary assumption made of a variable 
nucleus-electron parameter. Since there is at 
present no theoretical basis for a rigorous calcu- 
lation of F factor for a nonspherical electron 
distribution, an approximation was made by 
assuming each electron shell displaced a slight 
but fixed distance along the c¢ axis, each shell 
remaining spherical. 

Using F factors of the various shells from the 
data of Pauling and Sherman,'® computations 
show that satisfactory scattering curves are 
obtainable with only the J/ shells displaced. 
From Bohr-Coster diagrams and the size-screen- 
ing constants of Pauling and Sherman, the 
energies and average radii of the 3,0; 3,1 and 
3,2 sub-shells of zinc are approximately : E= 136, 
90 and 12 volts per electron; and R=0.25, 0.28 
and 0.39A, respectively. Keeping the relative 
displacements of the shells in rough conformation 
with these data, curves were constructed for 
various displacements. Denoting the displace- 
ments of the electron shells of zinc and oxygen 
by dn: and b,), respectively; the set of values, 
a,;=a.=),=0; 43,0= @3,,= 0.01 ; a3,2= 0.04= be 
corresponding to the model shown gave the 
dotted curves in Fig. 6. Thus distortion of the 
M shells alone,* if sufficiently great, may account 
for the observed intensities. 

The rather large changes produced by these 
displacements are, of course, due to the 180° 
phase difference which exists between the nuclear 
and electronic scattering. At small angles where 
F approaches Z, an additional small phase shift 
in the summation has a large effect in the sum. 
These electronic displacements have little effect 
on the x-ray intensities. X-ray scattering curves 
computed with the same parameter values do not 
deviate by more than 5 percent from those 
calculated without displacements. 


18 L. Pauling and J. Sherman, Zeits. f. Krist. 81, 1 (1932). 
* We are indebted to Dr. Guido Beck for a valuable 
discussion of this point. 
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It has been suggested that calculation of the 
scattering produced by atoms slightly extended 
from the spherical form along the tetrahedral 
bond direction, as would be consistent with the 
peculiar tetrahedral binding in ZnO, might lead 
to curves in agreement with the experiments.* 
Rough calculations of the effect of displacing 
part of the zinc and oxygen scattering power 
along the four bond directions, instead of along 
only the ¢ axis as formerly, have been made by 
the writer. Although such considerations are 
probably poor approximations of the actual 
tetrahedral effect, it seems quite conclusive that 
as long as the scattering centers are equivalent 
and displaced equally, the general features of 
the observed curves cannot be reproduced. A 
dissymmetry among the four directions approxi- 
mately equivalent to those found in the earlier 
trials seems to be essential. 

Considerations of the possible explanations of 
the anomalous intensities present in all electron 
diffraction photographs of zinc oxide, lead to the 


* We are indebted to Dr. H. A. Bethe for a discussion of 
this possibility. Theoretical calculation of the electron 
distribution in tetrahedral binding is now in progress in 


. these laboratories under the direction of Dr. L. Nordheim. 


conclusion that distortions of the atomic elec- 
tron clouds, which are plausible in principle, 
must be assumed. The degree of distortion 
necessary for explanation of the experiments is 
greater** than would be expected, unless the 
detailed quantum-mechanical investigation of 
electron distribution gives an unexpectedly large 
dissymmetry. 

Effects similar to those observed in ZnO are 
apparently present in the diffraction from BeO, 
which is of the same crystal type, although 
sufficiently good photographs for a trustworthy 
analysis have not yet been obtained. Indications 
of the effect have been found in diffraction from 
ZnS and probably also exists in the case of MoS» 
as noted by K. Lark-Horovitz. 

My thanks are due Professor K. Lark-Horovitz 
for the suggestion of this problem, and for his 
constant aid throughout its course. My thanks 
are also due Mr. A. I. May and Mr. W. C. 
Holden, of the Physics Shop, for the great care 
with which they have carried out the design of 
the apparatus used in the investigation. 


** This has been pointed out to us in discussion by Dr. 
G. Beck. 
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Further Evidence for the Existence of K“° 


Recently Nier' has obtained a peak at the 40 position in 
the mass-spectrographic analyses of the positive ions 
generated in potassium vapor, which he has been able to 
assign to K**, Mass spectra of positive thermions emitted 
from various potassium sources have been obtained in 
this laboratory which confirm the existence of this isotope. 

The potassium ion source in the present experiments 
was a platinum disk, 2 mm in diameter, held in position 
6 mm above the immergent slit by a tungsten filament to 
which it was spot-welded. The disk was impregnated 
with potassium by first coating with a potassium-con- 
taining material, basalt and KsPO, being used, and then 
heating to redness in vacuum for about one hour. After the 
heating the disk was washed in distilled water to remove 
any of the coating which might distort the accelerating 
field. During the heating process the platinum dissolves 
sufficient potassium to become an excellent emitter of 
K* ions. 


1A. O, Nier, Phys. Rev. 48, 283 (1935). 
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LETTERS TO THE EDITOR 


Prompt publication of brief reports of important 
discoveries in physics may be secured by addressing 
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twentieth of the preceding month; for the second 
issue, the fifth of the month. The Board of Editors 
does not hold itself responsible for the opinions 
expressed by the correspondents. 


The accelerating potential between the disk and slit 
was maintained at 580 volts. The pressure, as recorded by 
an ionization manometer, was 3 X10~® mm. 

A representative analysis of the region between the K* 
and K* peaks is shown in Fig. 1. 

It will be observed that the peak is flat-topped and that 
the background is comparatively small. 

The K**/K* abundance ratio has been measured for a 
large number of inorganic sources and has been found to 
be close to 14.25 with a probable error of about 0.03. 
The K**/K* abundance ratio estimated from curves of 
the type shown in Fig. 1 lies within the range of 8300+100 
in most every instance. The principal source of error is in 
the estimation of the background. These ratios are close 
to those obtained by Nier who found K**/K* = 13.96+0.1 
and K**/K*° = 8600. 

The possibility of the observed peak being due to Ca*® 
is small since large changes in the calcium content of the 
source had no effect on the observed ratio; the measure- 
ments were made at a filament temperature below which 
a Ca* emission could be detected. The possibility of the 
peak being due to K*® H* is remote. It could not be due 
to a hydride formed about K* ions in transit between the 
disk and slit since the chance of such formation is small 
and since ions so formed give a diffused and not a sharply 
defined peak. Again it is difficult to ascribe the peak to 
K** H* ions emitted from the platinum since fresh filaments 
gave essentially the same ratio as did those that had been 
heated for 12 hours at 700°C; the high dissociation pressure 
of potassium hydride should materially reduce its concen- 
tration during this heating period. The fact that the 
abundance ratio is essentially the same as that obtained 
by Nier from an entirely different type of ion source 
strongly suggests the existence of an isotope of potassium 
of mass 40. 

A. KEITH BREWER 

Fertilizer Investigations, 

Bureau of Chemistry and Soils, 
U. S. Department of Agriculture, 
Washington, D. C., 
August 30, 1935. 
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